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SOVIET PUYSICS: =GCRYSTALLOGRA PHY 


FROM THE EDITORS 


Translated from Kristallografiya, Vol. 
p. 3, January-February, 1956 


This country was one of the leaders in crystal- 
lography even before the revolution. Gadolin made 
a considerable contribufion to the subject with his 
derivation of the 32 symmetry classes, while 
Fedorov is world-famous for his work, including 
the deduction of the 230 symmetry groups for crys- 
tal structures, the invention of the theodolite goni- 
ometer and universal stage for optical examination 
of crystals, and the founding of chemical crystal- 
lography. Wulff made major contributions to physi- 
cal and geometrical crystallography. Later de- 
velopment was even more extensive, however, be- 
cause so many opportunities arose to apply the sub- 
ject in industry and in neighboring disciplines. 
The publication of a journal devoted to crys- 
tallography is a natural consequence of this 
development. 

The journal will cover the following areas. 

1. Theoretical crystallography: symmetry, 
scattering theory, physical theories of solids, theory 
of crystal growth, computational techniques. 

2. Analysis of crystal structures: x-ray,elec- 
tron-diffraction, neutron-diffraction, and other such 
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methods. 


3. Chemical crystallography (crystal chem- 
istry): of organic, inorganic, and mineral com- 
pounds; chemical bonding in crystals; methods of 
describing crystal structures; order—disorder ef- 
fects. 


4, Physical crystallography: optical, elec- 
trical, mechanical, magnetic, and other properties 
of crystals and textures. 


5. Physicochemical crystallography: poly- 
morphism, isomorphism, enantiomorphism, phase 
transitions, liquid crystals. 


6. Crystal growth: nucleation, growth, equi- 
librium, morphology, structure, goniometry, and 
dissolution of real crystals. 


7. Applied crystallography: methods and ap- 
paratus for crystal growing, crystal working, and 
technical uses of crystals. 


8. History of crystallography and the training 
of crystallographers. 

The journal will also carry brief communica- 
tions, reviews, discussions, and conference reports. 
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SOVIET PHYSICS —"CRYSTALLOGRAPHY 


COLOR SYMMETRY GROUPS 


VOSS. NOs 1 JAN.-FEB., 1956 


N. V. Belov and T. N. Tarkhova 


Institute of Crystallography, Academy of Sciences of the USSR; 


and Gor'kii State University 


Translated from Kristallografiya, Vol. 1, No. 1, 


pp. 4-13, January-February, 1956 


Original article submitted November 1, 1955 


A new method of selection from the 230 space groups is given for the 46 black- 
and-white (two-color) groups; this can be extended to 3, 4, and 6 colors; 15 of 
the color groups present in the 230 space groups are illustrated by specimens 


of infinite mosaics. 


Mosaics of low symmetry with 5, 7, or more colors that 


do not fit into the 230 space groups occur because there can be more than 14 
Bravais lattices in a space in which one dimension is qualitatively different from 


the other two. 


The infinite two-dimensional antisymmetry 
groups [1], or, in the usual representation, the 
black-and-white groups, were first derived by se- 
lecting from the 230 space groups those in which a 
figure remaining in a plane gives rise to derived 
figures, some of which may be reversed (turned 
over). This means that, if the figure (€.g., an ir- 
regular triangle) is viewed from the side facing us 
(colored white on that side) in the initial position, 
then in half of the derived positions it shows its 
other (black) side. There are 80 groups that leave 
a figure in the same plane [2-5]; deleting 17 one- 
sided planar symmetry groups (which do not reverse 
the figure), and the same number of groups in which 
each figure consists of two identical figures super- 
imposed (principal plane a mirror plane; the gray 
groups in Shubnikov's terminology), we are leftwith 
46 plane two-color groups in the strict sense of the 
word. 

These groups can be deduced without resort to 
the 230 groups by the method used [6] in deducing 
the 1651 black-and-white (Shubnikov) antisymmetry 
space groups. To the well-known five planar 
Bravais lattices we add the two-color planar 
Bravais lattices, namely: a monoclinic one PE (Oe 
p') in which any edge (or diagonal) of the element- 
ary parallelogram is centered in the colored way 
(i.e., with the other color); three orthorhombic 
(rectangular) ones: Py, (one edge centered in the 
colored way), pd (centered-in the other color with 
respect to the area of the elementary rectangle), and 


c' (both edges of the ordinary centered rectangular 
cell centered in the colored way). The last planar 
color lattice p{ is square and centered in the colored 
way. There are thus five color planar Bravais lat- 
tices in addition to the ordinary five, or ten in all 
(Fig. 1).! 

We give below the international expressions for 
the 17 ordinary plane groups [7] containing not more 
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Fig. 1. The planar two-color Bravais nets (lattices), 


1See the paper cited above for proof that these lattices exhaust the 
possible two-color ones in a plane (also for the 36 in space). 
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Fig. 2. The 46 two-color groups. 


COLOR SYMMETRY GROUPS 


TABLE 1. List of the Planar Two-Color Groups 


itial 
Rational Cochran's fInitial space group rs oe eee ge 
No. : containing related 
symbols symbols with reversing 5 
: objects at 
symmetry elements 
two levels 
1 2 3 4 
4 | pe’ p2’ Cri On Bidar 
ll ae pt C, = PAlb CL=A 114 
' ’ 9) 
3 Py2 pt Coy = Pl > GL = A112 
4 | pm’ pm’ Oh 12M Ca IPH 
Sa ire. Pe’ C2=P 12,4 C2 = Palt 
8 | em! em’ Cr i C2 Ce 
7 | Pom pm 0’ C2, = Pmb C3 = Amit 
Pos pst! C3) = Pb2b C1 = Abit 
Py im pm +m’ Cra Pome Cl = Aim 
10 Py 1g ps +38’ C2 = P2xat C2 = Alas 
41 | pom pm +g’ C2, = Pm2yn C3 = Imit 
12 | Pos pgm’ C8, = Pb2n C4 = Ibt4 
13 | o'm cm -+ m! CP = Cm2a C3 = FmAt 
14 pmm’ pmm OL = pk CA Pine 
45 | Pmim’ pm'm' Di = P222 C2 Pee 
, , 2 : 
16 | pms pm'g C4, = Pi—i | CS,= Pea 
2 
17 | pmg’ pmg’ Cop = PS Cee Bare 
/ 2 one 
18 | pegs" Pee Cy =P Coy = Pbn 
19 | ps's" 88" De P27252 Cl = Pnn 
20 pin'g’ pm’ g’ DE == P2222 OR. = Pen 
2 
24 cmm' cmm Oh. = Ce C32 = Cme 
22 em’m' em'm DS = C222 C32 = Cee 
23 pmm pm,m-—-m’ D3, = Pmmb C34 = Amm2 
24 | pymg pm,e+e' D!} = Pmab CM — Ama2 
25 | pysm pg, m +m’ D3, = Pbmb C13 — Abm2 
26 | Pgs Peg +s D8, = Phab C}? = Aba 
Ail Pomm pm+g’, m+ g’ DP = Pmmn Ces = Imm 
28 Pong pm-+g', g+m’ D3, = Pman OS = lhe 
29 Pose pg+m’, gtm’ D§, = Phan O27 = Iba, 
30 c/mm cm +m’, m+ m‘ D3} = Cmma C25 = Rian 
34 py’ pe’ Ste Pe Ci = Pay 
E { 4 
32 Pod p4t CP ChE 
33 p 4’ mm’ Pp 4 mm’ D5 = P4m2 GE = P4,mce 
34 p 4 m'm Pp 4’ m'm Dig == P42m Css = P4.cm 
35 p4am'm' Pp 4m'm’ Dy = P 422 Opie Aco 
36 p 4 gm’ p 4’ gm’ Die =P4g2 Cx = P 4, gc 
37 p 4 gim p 4 gim D3, = P42.m Om = P4,nm 
38 Pp 4g'm' p 4gim’ Dy P42 Cin Pane 
39 Pc 4mm p 4m+g’, m+ m’ Df, = Pamm C2, =14mm 
40 Po 4gm p 4g+m’, m+m’ DE yn fom on =L4gm 
41 | p 3m'4 p 3m'4 Di = P312 C2 =P 3 
42 | p 31m’ p 31m’ D? = P 324 Of, = PBhie 
2 PB po’ Clie BS C8 = Pb, 
= 2 
44 p 6’ m'm p 6’ m'm Dis Prot = Gr = P 63cm 
45 p 6’ mm’ p 8’mm’ e == JPR & 1 Ce = P63 mc 
46 p 6m'm’ p 6m'm' Di = P 622 C2, = P bec 
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than two principal generating symmetry elements; 
these we make colored (Shown by primes) in turn, 
and then both together. 

The theorem proved in the original derivation 
cited above indicates that any symmetry element in 
the direction of any color translation (in the two- 
color translational groups) either alternates with a 
(similar) uncolored element or coincides with the 
latter, so the complete list of the two-color groups 
is given by replacing in each of the 17 ordinary 
groups the uncolored symbol for the Bravais lattice 
in turn by all the colored ones of the same system. 
The two-color Bravais lattice may not relate in the 
same way to the generating elements in positions 1 
and 2 in the symbol, in which case the number of 
corresponding groups is doubled. Table 1 gives all 
46 two-color groups arranged in the order of deduc- 
tion. The first column gives the representation in 
terms of colorless and two-color Bravais lattices, 
which we consider to be the most logical. The se- 
cond column gives Cochran's rather illogical sym- 
bols from his paper in Acta Crystallographica [8], 
while the third column gives (in Schoenflies and in- 
ternational symbols) the space group (out of the 230) 
that corresponds to that color and which was the 
initial one in the original derivation of the 80 two- 
sided two-color groups. The last column gives a 
further space group from the 230 from which the 
two-color group may be deduced (see below for de- 
tails of this). 

Figure 2 gives schemes for all 46 two-color groups, 
which have been given before: in the first derivation by 
Alexander and Hermann [2], by Weber [3], and then by 
Shubnikov inhis "Atlas" [4] and "Symmetry" (5]. Some 
of the schemes do not agree with those given in these 
sources, firstly because, in these sources, the or- 
thorhombic lattice is centered (in the colored way) 
either on the a—p}, edge (mostly) or on the b— p} 
edge, while we use always pj; secondly because 
preference was sometimes given in the earlier 
works to the color symmetry elements (i.e., the 
corresponding plane or axis was brought into coin- 
cidence with a coordinate one) when color symmetry 
elements alternate with uncolored ones (in groups 
with colored Bravais lattices). We have given pre- 
ference to the uncolored elements, by virtue of the 
basic theorem; for instance, group pag is given as 
pis on account of the alternation of g with m' (Fig. 
3a), whereas, in the earlier works, it is given as 
pgm' (Fig. 3b); group piss (Fig. 4a) is shown as 
pim'm!' in Fig. 4b. 

Antisymmetry was the basic concept in the 
early derivations of the two-color groups, especial- 
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ly by Shubnikov (front and rear in a plane figure, or 
plus and minus sides); this is not represented other 
than by implication if the two sides are given differ- 
ent colors. The groups may be deduced in another 
way, namely by selecting from the 230 space groups 
those that place all figures derived from a givenone 
(an irregular triangle) at two levels only. Let the 
height of the three-dimensional cell be c, and the 
height of the initial figure be z; then half of the fig- 
ures related by the symmetry of the group will re- 
main at level z, while the other half will lie at z 
+c/2. However, z can always be taken as 0 insuch 
groups. The symmetry elements that do this are 
screw axes (2,, 4), 63) and c and n symmetry planes. 
In addition, the A, B, I, and F Bravais lattices can 
act in this way. The figures at z(0) are given one 
color, andthose at(z+c/2) or c/2 are given the other. 

This method has the advantage of completely 
avoiding the one-color plane groups and the gray 
ones; it gives only the 46 two-color groups. Some 
of the aspects of the method are given below. 

It is reasonable to extend this method to figures 
with more than two colors. The groups of multi- 
color symmetry are derived by isolating from the 
230 those that place the figures at three levels only, 
at four levels only, and at six levels only, these 
levels being denoted by different colors. 

The concepts of classical crystallography rule 
out five-color, seven-color, and other such sym- 
metries (see below). The color symmetry elements 
(for numbers exceeding two) are 3,, 35, 44, 43, 64, 
65, 6), and 6, axes together with a d plane and 
Bravais lattice R. 

The tetragonal color groups are P4, (I, I), P4s, 
14, (III), and multiples of these groups, namely 
14, md (IV—V), 14, cd (VI). The trigonal subsystem 
contains P3, (X), P3., R3 (XI), R3m (XII), R3c (XID) 
while the hexagonal one has P6, (XIV), P6;, P6, (XV), 
P6,. The sole group in the lower systems is the 
orthorhombic Fdd (VII, VII, IX), which has figures 
at four levels. 

We give color mosaics consisting of fairly large 
numbers of cells for llof these 15 groups (one of 
each enantiomorphic pair) as a separate table; the 
triangles used in the analogous schemes for the 46 
two-color groups do not give a clear conception of 
these. The symmetry is given under each mosaic, 
while the above lists of groups contain (in parenth- 
eses) the Roman numerals representing the corre- 
sponding mosaics. Many of these mosaics are tobe 
found in inlaid floors, etc., but some of them ap- 
pear novel and have been constructed solely from 
the graphs of the corresponding space groups (14,cd 
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Fig. 3. Group pig in different aspects. 


—VI, R8c—XII). 

It is convenient to place a piece of tracing paper 
over the mosaic and mark in a unit cell, entering 
the symbols for the symmetry elements [7, 9]; the 
color action of these elements can then be followed. 

A 4, axis places the figures at four levels,which 
in our figures are shown as yellow, blue, red, green, 
and then yellow again; operations 4, or 2, take us 
from yellow to red or from blue to green, as does 
the action of a c- or n-glide plane. The six colors 
for a 6; axis are white, red, green, black, yellow, 
blue, and white again. The 2, axis present in a 6, 
axis corresponds to the transitions white —black, 
red—yellow, and green—blue. These latter are al- 
so effected in hexagonal groups by c- and n-glide 
planes. 

We at first thought that only the two enantio- 
morphic groups 6; and 6; were possible for six- 
color symmetry, but an examination of space group 
R38ec gave us mosaic XIII not so far found in inlaid 
floors; this showed that the three levels generated 
by the triad screw axis are further divided in twoby 
ac plane, and mosaic XIII may be split up either in- 
to hexagons of three kinds (each containing two 
colors with a difference of c/2) or into triangles of 
two kinds, one containing three colors with differ- 
ence c/ 3 and the other also three colors with the 
same difference, but the complementary colors dif- 
fering by c/2 from those in the first kind. The hex- 
agons contain triad rotation axes at their centers; 
the triangles, triad screw ones. 

The simpler pattern IV was obtained at first for 
group 14, md, and it seemed that V represented a 
new symmetry group; but the two patterns have the 
same symmetry, the mirror planes in IV being pa- 
rallel to the sides of the square plates, whereas, in 
V, they run along the diagonals of these plates.* 

The orthorhombic mosaics corresponding to 
group Fdd are derived from the tetragonal ones with 
d planes by rhombic deformation, i.e., transforma- 
tion of the square plates into rectangles or into 
rhombic ones, in accordance with whether the d 
planes lie along the sides of the plates (as in IV) or 
along the diagonals (in V). Elongation and compres- 
sion always occur along the d planes, which persist 
in the orthorhombic symmetry and there define the 
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Fig. 4. Group pigg in different aspects. 


coordinate directions, so the initial tetragonal cell 
must be taken in aspect F, not I. We arrive at the 
sole orthorhombic group Fdd not only from group 
14, md= F4,dm (mosaic VII from IV, mosaic VII 
from V), but also from 14, cd= F4,de (mosaic IX 
from VI). 

The lower systems contain (in classical crystal- 
lography) no symmetry elements (apart from d 
planes) that possess figures related by them at more 
than two levels, so it would appear that only group 
Fdd allows of multicolor groups with cells neither 
square nor hexagonal. But such groups (not space 
groups) do exist and are often found in mosaics; in 
addition, we find in them sets of 5, 7, and more 
colors forbidden by classical crystallography. Sev- 
eral such are shown, the first five (XVI to XX) be- 
ing in colors, and the others being shown in the text 
as arrays of plates to be treated as colored in ac- 
cordance with the numbers (the numbers of the levels 
in the other interpretation). 

Mosaics XVI and XVII are similar in pattern, 
being centered with a symmetry plane (Cm) and hav- 
ing odd numbers (3 and 5) of colors; the second is 
forbidden by classical crystallography, but it is ob- 
vious that similar mosaics with 7, 9, and more 
colors are possible. Figure 5 shows schematically 
(in black) a seven-color one. Mosaic XVIII is ana- 
logous as regards symmetry and number of colors 
to XVI, but it consists of rectangular plates instead 
of rhombs. There are 6 colors in mosaics XIX and 
XX. 

The symmetry of the first is Cm, while that of 
the second is Im; Fig. 6 shows a mosaic analogous 
to XIX but having eight colors, while Fig. 7 shows 
eight colors and no symmetry element apart from 
the triclinic F lattice. Figure 8 shows a cell of a 
mosaic consisting of six colors which has a primi- 
tive triclinic cell but an obvious law of succession 
in the colors. 

These mosaics (not of space-group type) arise 
from special causes; lattices I and F\canbe referred 
in the monoclinic system to ones centered only on 
one rectangular face, on account of the absence of 


*See correction at the end of the paper. 
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the requirement that the cell should have a right 
angle between the third axis and the centered face. 
Color symmetry replaces the continuous third di- 
mension by a sequence of colors, which thereby 
eliminates the question of the various possible 
angles between the third axis and the rectangle (b 

x a) formed by the first two. The impossibility of 
making these angles different causes the third di- 
mension to be perpendicular (or in Some other way 
specified) relative to the ba plane. This third axis 
does not have to be the direction of the shortest 
translation out of the plane of the other two if its di- 
rection has been specified; the ca (or cb) plane may 
be centered in two or more ways, not merely one. 

Consider first an even number of centered 
points, especially 2. In this case the diagonal of 
the side rectangle is divided into 3 parts (8 levels), 
which enables us to construct mosaic XVI (three 
colors), which is the one most commonly encoun- 
tered in inlaid floors. The side "rectangular" face 
may be centered by its four corners, in which case 
we have five levels and arrive at mosaic XVII con- 
taining five colors. In both cases the cell of the 
mosaic is centered, but the two side faces are un- 
centered. 

In general, a mosaic with an even number of 
colors may have an odd number of centered points 
on face A; then, if face C is centered, the basic 
theorem of Bravais lattices shows that face B will 
be centered also (mosaic XVIII, Fig. 6). The last 
two examples have a symmetry plane each, where- 
as the sole symmetry element in Fig. 7 is lattice F. 

The 14 Bravais lattices arise when we select 
six elementary parallelepipeds on the basis of the 
principal directions in crystals, these figures cor- 
responding to the setting rules of classical crystal- 
lography (i.e., to the six systems); we then find that 
the vertices of such cells cannot take up all the lat- 
tice nodes. There is no reason why the third direc- 
tion in a monoclinic cell (one special direction) 


Fig. 6. The eight-color mosaic AYIgcm. 


should be perpendicular to the plane of the first two 
when the second direction has been taken perpen- 
dicular to the special one; on this basis we get only 
two monoclinic cells, to which all others can be re- 
duced. Replacement of the third direction on color 
principles essentially makes that direction other 
than perpendicular, one result of which is the exist- 
ence of cells centered in several ways. Colored mo- 
saics thus illustrate Bravais lattices 15, 16, etc., 
which are not reducible to the basic 14. 

This means that the lattice symbols cannot in- 
dicate only the faces centered in the sense of the 14 
Bravais lattices (once and at the center of the face). 
The multiple centering must be allowed for and 
must appear in the formula. For instance, the 
formula for XVI will be Allpc, for XVII and XVIII 
it will be AIVpC, and for XIX itis AVBC. The 
formula for Fig. 7 is (F) AIJBC. Indication of F is 
unnecessary in the last two cases, but the formula 
for XX, (1) AIlpp, requires the (1), for centering in 
the volume does not exclude centering of the faces 
(even number of nodes). The triclinic scheme of 
Fig. 8 has the formula AVplicI, The standard 
rule for the 14 Bravais lattices is replaced by a 
more general one for the color case: an odd num- 
ber of centering nodes in two faces of the cell im- 
plies an odd number of nodes in the third face also. 

Similar ideas appear in the early part of this 
paper in the third method of deducing the 46 color 
groups; Table 1 contains the side-centered triclinic 
cell Al11, monoclinic cells centered on the face 
normal to the special direction (Am11 and Ab11), 
and even the monoclinic Fm11 face centered all 
around. All of these cells are impossible (are re- 
duced to simpler ones) if the third edge is not per- 


Fig. 7. The eight-color mosaic 
(F) Alllgc, 


Fig. 8. The six-color mosaic AY pil. 
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pendicular to the mutually perpendicular edges of 
the monoclinic cell, so even the two-color sym- 
metry groups contain more than the 14 independent 
Bravais lattices (ones not mutually reducible). 
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This paper contains an error in the uniform de- 
formation of the tetragonal mosaics into orthorhom- 
bic ones. Mosaic IV should be compressed along 
the diagonal of the square, white mosaic V should 
be compressed along the sides, not the other way 
around. The same error occurs in the curvilinear 
outlines in mosaic XI, which thus illustrates group 
R3m instead of R3. 

These faults have been corrected in the table 
in our second paper on the color groups, which ap- 
pears on p. 487. In addition, new mosaics having 
only rectilinear boundaries in their components are 
shown there for all mosaics having curvilinear 
boundaries in the original paper. 
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A theory of the author is used to set up the simplest relations and tables for 
the certain and probable determination of the sign of structure products. The 
graphical and tabular material here presented may be successfully used to de- 
termine the signs of the structure amplitudes for crystals with cells containing 


up to 50 or 70 effective atoms. 


1. The theory of structure-amplitude relations 
was developed by the author in papers published in 
Doklady Akademii Nauk SSSR in 1954-1955. In this 
article, we shall formulate very simple practical 
rules for determining the sign of structure ampli- 
tudes. 

2. The theory developed enables us to consider 
relationships involving m amplitudes. In practice, 
however, calculations involving more than four am- 
plitudes must be entrusted to a computer. The 
equations and rules governing relationships involv- 
ing four amplitudes are set forth in the papers men- 
tioned above. 

Here we shall consider only the simplest case: 
relationships involving three amplitudes. Theory 
shows that such a relationship exists for the ampli- 
tudes of three reflections with indices h,h,hg (H), 
kjk kg (K), and kj —h,, k,—hy, kg—h3 (K—H). 

3. This relationship consists in the following. 
The value of the structure product FuoikiK_H: 
where F is a unitary structure amplitude, is posi- 
tive definite for large amplitude values. The 
structure product is certainly positive if it islarger 
than '4 in absolute magnitude. However, structure 
products smaller than de may also be certainly po- 
sitive. Exhaustive data for the certain determina- 
tion of a positive sign for eens ts are given by 
the graph presented in Fig. 1. Along the axes lie 
the values of Fry and ee and each curve corre- 
sponds to a specific value of PK. The following 
rules govern the use of the graph: 
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(a) Choose the curve that corresponds to the 
given value of lFi_nl; plot the point with 
coordinates |Fy], |Fx| on the graph; if 
this point lies above the curve, the struc- 
ture product will certainly be positive. 

(b) Even if the structure product is not certain- 
ly positive, it still has a probability W,>¥4 
of being positive. 


Graphs of W, are given in Figs. 2 and 3 for 
several cases from the author's formula 


Here o2= D(Z2/ Z*), where Z is the number of 
electrons in the cell and Z; is the number of elec- 
trons in the atom. If all the atoms are identical, o 
= 1/VN, where N is the number of atoms in the cell. 
It is convenient to call 1/o? the effective number of 
atoms in the cell. 

The figures show curves of W..as a function of 
(fae yl- In order to keep the graphs two-dimen- 
sional, we had to set lFyl = = [file It is noteworthy 
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that a nearly asymptotic approach to W,= 1 takes 
place only for large values of lf}. 

4. The value of o is a very important constant 
of the crystal. If it is known, the structure can be 
solved by direct methods. 

Since the distribution of structure amplitudes, 
as Wilson demonstrated, is Gaussian, 31.74% of the 
experimentally determined amplitudes are greater 
thano, 4.56% are greater than 20, 0.27% are 
greater than 30, and 0.01% are greater than 40. 
Analysis has shown that structures with values of o 
smaller than 0.1 cannot be solved by direct methods 
while the solutions of those with 0 > 0.17 are trivial. 

5. To calculate F from experimental intensity 
data, we must know the curve of the mean atomic- 


05 — ——e oe — — a. 
CE OS VEE TE TE ET ia 
Nef 100 x-F,-F, 

Fig. 3. 


temperature factor. If the number of reflections 
with which the investigator is operating is small, 
the curve of the atomic-temperature factor must be 
taken from published data. For a number of reflec- 
tions of the order of several hundred, one should de- 
velop the reflection sphere in the region and find the 
F’ curve as a function of (sind YA. According to 
theory, this will also be the mean f -curve given on 
a relative scale. 

With the help of this curve, the unitary struc- 
ture amplitudes may be found apart from a constant 
factor. The constant factor is chosen in such away 
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that the number of amplitudes larger thano, 20, fully applicable to crystals. Its effectiveness is de- 

etc., coincides with the theoretical values indicated termined by the number of atoms in the cell, and not 

in paragraph 4. by the number of crystallographically nonequivalent 
6. What we have discussed exhausts the me- atoms. Second, for reliable results a three-dimen- 


thods of determining structures with o > 0.17-0.18 sional experiment is necessary. However,in simple 
(up to 30 identical atoms in the cell). We make just cases (0 > 0.2), investigations may be successfully 
two further comments. First, the direct method is carried out on the basis of projections. 
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Kinematic electron scattering is discussed for single crystals, which give spot 
patterns; the integral intensity for the hkl spot from an ideal crystal is given by 
formula (8b), which becomes (9a) and (9b) for the angle of regular reflection. 
Ideal crystals are very seldom encountered; the specimen is usually a mosaic, 
whose principal characteristic is the angular distribution function over a certain 


range. 


Formulas (18a) and (18b) give the integral intensity for the reflections 


from a mosaic film; the factor arising from the mosaic structure is of the form 
dhkj/@, in which dy is the interplanar distance for that reflection, and @ is 


the effective angular spread of the crystals. 


tween the kinematic-scattering formulas (which apply for thin blocks) and the 
dynamic-scattering ones (which lack the proportionality to thickness or volume 


characteristic of the kinematic ones). 
cases. 


Introduction 


A necessary step in a diffraction study of a 
crystal structure is conversion from the observed 
intensities to the structure amplitudes. The kine- 
matic scattering theory is entirely satisfactory for 
x-ray work, but dynamic effects are very much 
more common in electron diffraction; nevertheless, 
the kinematic theory is very important in the latter 
area, because, very often, the specimen is such as 
to show that scattering. 

Here we examine kinematic electron scattering 
by various types of specimen and deduce the inten- 
sity formulas for the various types of pattern. 

The principal condition for kinematic scatter- 
ing is that the scattered beams should be weak rela- 
tive to the incident beam; secondary scattering can 
then be neglected, as can energy loss in the inci- 
dent beam. 

The solution to Schrédinger's equation in the 
kinematic approximation [1,2] gives the intensity 
scattered in a direction s, relative to the incident 
intensity, as 
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WiOi irae NiOre gt JAN.-FEB., 1956 
The relationship is discussed be- 
The mosaic factor is the same in both 
J (8) = Jo] B (8) P/r?, (1) 


in which ® (s)is the scattering amplitude, which is 
defined by the Fourier integral of the potential g(r) 
of the object [3]: 


a 


P(s)=K \ 6 (ryets” do, Ae = se, (2) 


i? 


As applied to a single atom, the above condition, 
which gives rise to (2), is simply the Born approxi- 
mation; the amplitude of (2) for a crystal is deter- 
mined by that of a single elementary cell, p);7 (the 
structure factor), and by the Laue interference 
function |D]?: 


J (h) = Jo| Pasi? | D (h) |? /7?, (3) 
in which h = s/27, and 
Drie fam, we (4) 
a 
Here the atomic temperature factor isf gp Fide 


exp{ —B(sin 6/2 )*} , f@ being the atomic factor in 
absolute units [4]. The interference function is 
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sin? tA,h, 
Dil ttee lly Gee a (5) 


; sin® nah, ’ 
i=1, 2,3 ( 


in which hj are the coordinates of reciprocal space, 
Aj are the linear dimensions of the crystal, and aj 
are the cell repeat distances. The sine in the de- 
nominator may be replaced by its argument if we 
consider ipl’ near a given node hkl and reckon ° the 
hj from its center; then each of the factors in (5) at 
the peak (hj = 0) takes the value 


si wA,h, / A, y2 
(ma,h,)? ee a | 2 (6) 


and the integral is 


SES 1G. 

* sin? TtA;h, 1) A, (7) 
— dh, =. 

\ (ta;h,)° 4 a? 


Intensity Scattered by an Ideal Crystal 
[5 | 


Consider the integral intensity Iy,z of spot hkl; 
we have to integrate (3) over two directions in the 
plane of a screen on which reciprocal space is rep- 
resented in a scale LA, X being wavelength and 
L@r) being the specimen—screen distance. From 


(7), 


It = \J (ly hgh s) Lidh,Lidhy 


A 


Ay/ 
2,2 
aya 


1, 


sin® tAgshg 


= Jo)? | Oxi? 


(8a) 


(Tagh3)? 


Let S = A,A, be the area of the crystal (volume V 
= SA3) and = a@1a,a3 be the volume of the cell. Then 


9 . 
* sin? wAghs 


(8b) 
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Taking the In,7 corresponding to the angle of true 
reflection, i.e., to intersection of the reflection 
sphere for node hkl at the maximum (h; = 0), we 
have from (6) that (8a) and (8b) may be replaced by 


(9a) 


or 


hkl f 42 (9b) 


J)S is the incident intensity received by area S, so 
the ratio of Ihkz to J)S is proportional to the square 
of the structure amplitude, to the crystal thickness, 
and to the wavelength. Formula (9a) illustrates the 
proportionality between I,,j and scattering volume 
characteristic of kinematic scattering. 


The curvature of the reflection sphere general- 
ly allows not more than three reflections to appear 
in the scattering by an ideal single crystal. Real 
specimens give spot patterns containing all the re- 
flections in a given zone, which can be explained 
only in terms of a mosaic of single crystals if the 
scattering is of three-dimensional type [6]. There- 
fore, (9a) and (9b) hardly ever correspond to actual 
experimental conditions, but they reflect the essen- 
tial processes in an individual ideal block of the mo- 
saic. The condition that (9b) equals one gives us 
[5] the critical thicknesses for ideal crystals as A4 
© (A |@p),/2|)"1, which correspond to transition 
from kinematic scattering to dynamic. These thick- 
nesses fall as the atomic number of the material in- 
creases, because f, increases and hence, also, the 
@hkj of (4); but they increase as the structure be- 
comes more complex (with the number n of the atoms 
in the cell) because this reduces the mean value of 
the unit amplitudes ®p,j, and hence of ®p 7 (see [5)). 
Simple structures (e.g., of cubic metals) give low 
critical thicknesses, e.g., about 50 A for Au, and 
about 200 A for Al. Moderately complex structures 
(n of 20-25) consisting of light atoms give Aj © 1200 
A, while for the same structures with medium- 
weight atoms Aj is 500-600 A, and for heavy ones, 
250-400 A. 

We assume that the beam is ideally monoener- 
getic and collimated in deriving (9a) and (9b); a 
real beam has an energy spread of about 0.1%, 
while the divergence at a single crystal does not 
exceed 0.0001 radian, so these have entirely neg- 
ligible effects on (9a) and (9b). 


Intensities of Reflections from a 
Mosaic Single-Crystal Film 


A real specimen giving a spot pattern is a mo- 
saic single crystal. The part of the crystal ex- 
posed to the beam contains many! ideal crystallites 
that are not strictly parallel (have a certain angu- 
lar spread). 

Oscillation or rotation at a certain angular 
velocity are used in x-ray work with microcrystals; 
the conception (which does not always correspond 
to the actual working conditions) is meaningless in 
electron diffraction, for the crystallites are at rest 
but the assembly of them has a set of orientations 
in space, which replaces oscillations (i.e., a setof 
orientations in time). 


’ 


‘For instance, the beam might have an area of 1 mm? = 10! ?, 
the film being 500 A thick and the mean volume of a crystallite 
200 x 200 x 200 A? #10" A; then there are about 5° 10° crys- 
tallites. 
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This set of spatial orientations has the follow- hkl, but it would be very difficult to consider this. 


ing consequences. First, not all the crystals are A particular and rather different case of nonuni- 
in the position corresponding to passage of the re- form distribution is a texture, when any orienta- 
flection sphere through the maximum of the inter- tion is possible around one of the crystallographic 
ference function in the third direction, so it is ne- axes. This case is considered separately in the se- 
cessary to integrate over that direction. Second, cond part of this paper. 
not all of the crystallites are in reflecting positions; General Formula. The mosaic film has 
some of them do not reflect at all, and the propor- /f (@) as its distribution function; importance at- 
tion participating may vary from one reflection to taches only to the part of f (a) corresponding todis- 
another. placement of the reciprocal-lattice nodes perpen- 
Distribution Function. Let M (Fig.1) dicular to the reflection sphere (plane of the screen), 
be a reflecting plane; all of its orientations (and so integration in the plane of the screen gives for 
hence those of the crystallites in the mosaic) may the integral intensity, by analogy with (8a), 


be described in terms of the function for the dis- 
tribution of the normal N to this plane over the 
angles a, and a3. There is also a distribution ina, 
(Fig. 1), but this is of no significance in relation to But here G is not |D(h;)|?, as it was in (8a), for itis 


5 Av As 
Fyxt = FQ)? | Onn |? = G. (11) 


2,2 
a4 


the intensity, because these rotations leave the also governed by f(a). We must find G by taking 
plane in the same position relative to the beam. Al- the integral of the interference function in the third 
so, @, does not alter the angle? between the plane direction over all possible positions given by f(a). 
and the incident beam, though the spot is thereby Figure 2a indicates the mutual disposition of f(a), 
displaced through an angle a, on the screen. Only the reflection sphere in a position # , and the inter- 
a, determines whether the given plane M lies in the ference function | D3(a@ — B)|? for crystallites in po- 
reflecting position, so only the component distribu- sition a, the last (as also f) being a function of 
tion function f (a) (discarding the subscript in @,) is angle. The intensity given by the crystallites in the 
significant in the intensity calculation. The film range @ to a +da is proportional to f (a)da (name- 
contains n identical crystallites, so the normaliza- ly to the number in that position) and to the value 
tion condition for f (a) is for the interference function at its intersection with 


the reflection sphere at point a—f. The integral 
\ /(@) da --n =V'/V, (10) over all positions of the crystallites gives G: 
in which V is the volume of a crystallite and V' is G(3) = \ D (« — 8) |? f (a) do. (12) 
the volume of the film. There are f(a)d _ crystal- 
lites in the range a to a+da. We can usually as- 
sume that f(a) is the same for all reflecting planes, 
especially if the crystal is of isometric habit; in 
general, f(@) should be taken to be dependent on 


eg | 


Reflection sphere 


Reflection sphere 


Screen 


Node 000 


Fig. 2. (a) Mutual disposition of f(), |D3|?, and 
reflection sphere; (b) the particular case when the 


Fig. 1. Scheme for the angular distribution of the crystallites 
origin is at the center of the reflection sphere. 


in a mosaic film. 
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The dependence on 8 shows that G(6), and hence the 
intensity of (11), is dependent on the position of the 
mosaic crystal relative to the reflection sphere, 
i.e., on the angle of rotation of the specimen. It is 
convenient to take the origin at the center of the re- 
flection sphere (i.e., to put 8 = 0; Fig. 2b) inorder 
to evaluate G; then (12) becomes 

G v D (a) |? f (a) da. (13) 
But here we should bear in mind that this fixes the 
reflection sphere at some definite point in the dis- 
tribution function; the dependence on the angle of ro- 
tation is not seen. The element da in (12) and (13) 
is related to the linear dimension dhg in reciprocal 
space by 


dhs 


la — m 


sche dhs dirt, (14) 
as Fig. 2b shows; Hhk/ in (14) is the distance to the 
given reciprocal-lattice node and dhk] is the inter- 
planar distance. 

Uniform Distribution. Here we calcu- 
late G for the very simple case of a uniform dis- 
tribution in angle (Fig. 3). The angular range is 3- 
4° (& 1s radian) for real mosaic films that give spot 
patterns. A uniform distribution makes f(a) con- 
stant within this range and equal to n/ (a, — a,)=n/a, 
in which a = a, — a, is the angular width of the range. 
Then the G of (18) becomes constant for all positions 
of the reflection sphere far from the edges of the re- 
gion (Q,, G9): 


Cees 


a 


\ sin? tAshs ae (15) 


(magh3) 
Then (7) and (14) give 


(16) 


This (reduced to one crystallite, n = 1) is less than 
the A3/a 3 of (6), because dy; j/a@ is less than A3. 
The intensity of a reflection from a mosaic film is 
therefore less than that from an ideal single crys- 
tal of the same size. Then from (11) and (16) 


A AgAg en je) 
ap etn . 
Zee a 
asa; 


Fn SF gh? | Daxi |? 


(17) 


But A;A,A3 = V and nV = V' from (10), in which V' 
is the volume of the film, so finally, 


2 Pi 
pony (18a) 


lo 


Dp jet 


= 2 
Tiwi a Joh Q 
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or 


2 td nyt 
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Dart 
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(L8b) 


which resemble (9a) and (9b); here, t is the mean 
thickness of the mosaic film,* which is not neces- 
sarily the thickness A; of the blocks (it can be 
greater or less), and § is the area of the film. 

The dpi of (18a) and (18b) may be elucidated 
by estimating the fraction q of crystallites partici- 
pating in reflection. The interference region has 
a half-width of 1/ A; (width 2/ A3); the angular width 
for the hkl reflection is then (2/A3)qHhkI = 
= 2dpit/ Ag, 80, 


= 2dnx | Ag 


a 


(19) 


For instance, A; = 200A and a = ‘4, give q close to 
lford®5A (near reflections) andq © 0.1 ford 
0,7 A (far reflections). The proportional rela- 
tion of q and of the Ink] of (18) to dhk 7 is explained 
by the fact that constant 1/ @ produces a constant 
linear width in the interference region, not a con- 
stant angular width, the latter falling in proportion 
to dhkj (Fig. 4). 

Nonuniform Angular Distribution. 
A uniform distribution is an ideal case, but it can 
be extended to a nonuniform distribution of density 
jf (a) subject to certain assumptions. 

It is usual for f(@) to have a much greater an- 
gular width than |D;|?, so the G of (12) (for the real 
shape of the reciprocal-lattice nodes) has almost 
the form of f(a); hence the f (@) for a real mosaic 
film may be derived from the observed intensity 
variation at the screen of the diffraction apparatus 
as the specimen is rotated very slowly (the same 
result is attained from a series of photographs at 
small angular intervals, e.g., 14°). The most 
characteristic case is the following: the spot 
gradually appears, the intensity rises to a maxi- 
mum, and then there is a fall to zero. The dis- 
tribution function has the bell-shaped form shown 
in} Figs sel, 2... and.o. 

We consider two extreme cases of this type of 
J (a); the first (Fig. 5a) is that where the range be- 
tween @, and @, [limits of range where f (a) > 0] is 
much greater than the width of the interference re- 
gion for all reflections. The region where the in- 
terference function is finite is much narrower than 
the angular range of finite f(a), so we can take 


*By S in (18b) we should understand the area of the beam, if this is 
less than that of the film; the V of (18a) is then the volume of the 
film exposed to the beam, not the total volume. 
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f(a) as constant within that region and extract it 
from the integral in (13): 


G = f(a) \ SER de (20) 


(Taghs)° 

Strictly speaking, the f (qa) extracted in this 
way will vary from one order of reflection to an- 
other (Fig. 6), because the reflection sphere has a 
certain curvature and the @);7 are different. This 
can be neglected if f(a) varies sufficiently slowly. 

The natural approach is to operate at maxi- 
mum intensity. The bell-shaped function of Fig. 
5a may be taken as having a peak height equal to 
twice the mean; then, from (10), 


H(X)max = 22 f(@) da = 2, (21) 
where a@'= a, — a, is the region in which f (a) is 
defined, i.e., the angular range between emergence 
and extinction of the interference spot for the ro- 
tating sample. From (20), (21), and. (16) it follows 
that, in this case, 


2D dyy 
As (22) 
oS 


and the effective a in (18a) and (18b) is a'/2. 


| D(a) I? 


flee 


Fig. 3. Uniform angular dis- 
tribution of crystallites. 


Fig. 4. Reduction in proportion 
of crystals participating in re- 
flection with increase in order 
of reflection. 


More rapid variation in f(a) (narrower angu- 
lar range) may cause (Fig. 6) the f( 1) to differ 
substantially. Samples of this type do occur, and 
then it is found that the near and far reflections at- 
tain their maximum intensities at slightly different 
angles of rotation. The best approach then is to 
take several photographs at similar angles and take 
the maximum intensity for each reflection; this cor- 
responds to the maximum (constant) value f (@)max, 
so the G of (11) will also be constant and may be 
evaluated from (22). 

The other limiting case is that where the limits 
of finite f (a) are much narrower than the interfer- 
ence region for all reflections (Fig. 5b), i.e., 
there is very little angular spread in the crystal- 
lites. Then (G6) shows that the interference function 
is close to A3/a% throughout the range (a,, @,), So 
we extract it from the integral and get from (10) 
that 


Xs 


G=4 \ f(#) da = nAs/a. (23) 
3 a 
Substituting (23) into (11), we get 
if @ 2 
Taste haleaa| ts (24) 


in which t is the thickness of the film. The blocks 
may form a single layer (t = A3), in which case we 
return to the ideal crystal of (9a) and (9b). This is 
entirely to be expected, because the mosaic ap- 
proximates to an ideal crystal as a decreases. 

The intermediate case (interference region of 
about the same width as the distribution function) is 
the most complicated; the intensities are then very 
inconvenient to use. It is clear that here G will ap- 
proximate to (23) when d is large, and to (22) when 
dis small. 


flOmar e/a." 
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Fig. 5. Nonuniform angular dis- 
tribution: (a) f(a) much 
larger than |D,|?; (b) the con- 
verse case, 
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Fig. 6. Effects of curvature of 
reflection sphere. 


Change in Type of Scattering with 
Specimen Thickness 


The kinematic formulas are applicable only 
when the blocks are thin; the dynamic ones apply to 
thick blocks. Here, we consider the relation be- 
tween the two sets of formulas as the thickness in- 
creases. 

Ideal Single Crystal. The dynamic 
theory [7, 8] (see also [2,6,9]) gives the integrated 
intensity at the reflection angle as 


Fess Preint eae (25) 
Here, I, denotes the incident intensity, namely J)S, 
where S is the area of the crystal. The vhk], H, 
and kof [2,6-8] correspond to our vyy = 47/1® 7/2 , 
H = A; (thickness of crystal), and k= 27/A. Then 
(25) becomes 


Tat =IoS sin? {2 | SE (26) 


J 


The argument of the sine is simply the quantity 
whose square determines the reflection intensity in 
the kinematic theory. Blackman [8] has shown that 
we may replace the square of the sine in (26) by the 
square of the argument for A; small, which gives 
(9b), which is a particular case of the general ex- 
pression (8b) for the kinematic scattering for h,=0. 

Formulas (8a), (8b), (9a), and (9b) are applic- 
able if As is less than the critical thickness Aj (see 
above). 

Mosaic of Large Blocks. This case 
has not been considered in the dynamic theory, but 
the integral intensity has been calculated for a 
crystal rotating with an angular velocity w [6,7]. 
Then (after averaging with respect to thickness), 


I TPL 
Py Oe 


Tne = (27) 
Rotation is not used in electron diffraction, but the 
mosaic structure produces an angular distribution, 
so w may be replaced by the effective angular 
range a (see above). Now, & ~ A/2dp 7 [note also 
the change in notation from (25) to (26)], so (27) 
becomes 


Lh pt 


ere (28) 


Tiny ) | Prag 
Ee Sa 


The linear relation to volume (thickness) charac- 
teristic of the kinematic theory does not appear in 
(26) and (28). The mosaic structure (which is in- 
dependent of the type of scattering in each block) 
in both cases gives the same factor dpxj/a [cf. 
(18a) and (18b)], so the other factor dependent on 
the structure of the specimen (the analog of the 
Lorentz factor for x rays) remains unchanged when 
we pass from the kinematic theory to the dynamic 
one; the square of A]6/Q| is merely replaced by the 
first power together with loss of the volume (thick- 
ness) dependence. The general formula will be 
intermediate between (18b) and (28) if the film con- 
sists of blocks of various sizes, some of which 
scatter kinematically and others dynamically. The 
spot reflections are replaced by Kikuchi lines if 
the blocks are very large. 


Thick Pilms-of SmalliBlocks. | Tet 
the size of the blocks be subcritical; the scattering 
in each is then kinematic, so (18a) and (18b) imply 
that Ink is proportional to t. The analogous form- 
ula for an isolated block is the rise in In,jz with Ag 
given by (9) up to Ink] ~ J)S, which means that we 
have the transition to dynamic scattering; (9a) and 
(9b) are inapplicable for As= Aj. Also, (18a) and 
(L8b) involve the assumption that A; < Ad, but here 
sufficiently large t can give the physically meaning- 
less result Ink] > J)S. The point here is that (18a) 
and (18b) require correction for secondary extinc- 
tion for t large, i.e., for screening by the layers 
of blocks preceding the one in question. At this 
point the treatment is no longer purely kinematic, 
because allowance must be made for the loss in in- 
tensity in the primary beam as it penetrates into 
the material; but no account is taken of the phase 
relations considered in the dynamic theory (prim- 
ary extinction). In addition, absorption cannot be 
neglected for thick layers. 

Secondary scattering can also occur [9-11] at 
the high scattered intensities encountered for thick 
films, no matter what the type of scattering within 
the blocks. 
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Conclusions 


Kinematic scattering has the scattered intensity 
proportional to the volume (thickness) of the speci- 
men and to the squares of the wavelength and the 
structure amplitude. A major characteristic of the 
Specimen is the angular distribution of the crystal- 
lites (mosaic structure). 

An ideal single crystal and a mosaic film with 
only small angular spread will have scattering de- 
scribed as appropriate by (8a), (8b), (9a), (9b), and 
(24); here the relative moduli of the structure amp- 
litudes are related to the relative integral intensity 
by [ne] =e 
spond to occurrence of only a few reflections from 
a zone) are hardly ever encountered in practice. It 
is most usual to find all the reflections in a given 
zone appearing simultaneously, on account of a 
range in @ commonly about 3-4° (mosaic). The 
factor for the mosaic structure is dy, 1/a, so from 


(18a), 
| = 
ial Sul ere nis 


A film of large blocks each having dynamic scatter- 
ing will have a factor for the structure of the same 
form as for a film of small blocks that scatter 
kinematically; in particular, for a mosaic film of 
large blocks, 


But these cases (which corre- 


(29) 


rel 


1 
| Dri | — Ro dice (30) 


It is very desirable to employ the intensities of 
spot patterns in structure studies, but these are 
sensitive to the distribution function and to the 
thickness, and there is always the possibility ofse- 
condary scattering, so it is always necessary to 


make a careful study of the type of scattering toes- 
tablish whether a given formula is applicable, the 
most common case being that described by (18a) 
and (18b), i.e., (29). 


LITERATURE CITED 


1. M. Born,Quantenmechanik der Stossvorgange, 
Z. Phys. 38, 803-827 (1926). 

2. M. von Laue, Materiewellen und ihre Inter- 
ferenzen (Leipzig, 1948). 

3. B. K. Vainshtein, Structural Electron Dif- 
fraction [in Russian] (Author's abstract, 1955). 

4. B. K. Vainshtein, "Atomic scattering factors 
for electrons," ZhETF, 25, 157-168 (1953). 

5. B. K. Vainshtein, "Kinematic scattering of 
electrons by an ideal crystal, " DAN SSSR, 

104, 537-539 (1955). 

6. Z. G. Pinsker, Electron Diffraction [in 
Russian] (Izd. AN SSSR, 1949). 

7. H. Bethe, Theorie der Beugung von Elektronen 
an Kristallen. Ann. Physik, 87, 55-129 (1928). 

8. M. Blackman, "On the intensities of electron- 
diffraction rings," Proc. Roy. Soc., A173, 
68-82 (1939). 

9. Z. G. Pinsker and B. K. Vainshtein, "Dynam- 
ic scattering of electrons, "Izv. AN SSSR, ser. 
fiz., 14, 212-222 (1950). 

10. J. M. Cowley, A. L.-G. Rees, and J. A. 
Spink, "Secondary elastic scattering in elec- 
tron diffraction," Proc. Phys. Soc., A64, 609- 
619 (1951). “a 

11. B. K. Vainshtein, "Some features of spot 
electron-diffraction patterns associated with 
secondary scattering, " Trudy Inst. Kristal- 
lografii, 12 (1956). 


SOVIET PHYSICS ~-CRYSTATLOGRAPHY 


VOLE tf NOE JAN.-FEB., 1956 


DERIVATION OF WORKING FORMULAS FOR THE ELECTRON 
DENSITY AND STRUCTURE AMPLITUDES FROM THE SYMMETRY 
AND ANTISYMMETRY OF TRIGONOMETRIC FUNCTIONS 


M. A. Porai-Koshits 


Kurnakov Institute of General and Inorganic Chemistry 
Translated from Kristallografiya, Vol. 1, No. 1, 


pp. 27-48, January-February, 1956 


Original article submitted October 1, 1955 


The derivations are considered for the various space groups in all settings on 
the basis of the symmetry and antisymmetry of trigonometric functions. Com- 
binations of these functions are derived having even and odd h, k, and /; these 
correspond to particular symmetry elements. The formulas for a set of sym- 
metry elements are derived from the general formula simply by eliminating 
terms contrary to the required symmetry. The relations between the Fourier 
coefficients are used only in the final stage [conversion of the limits of summa- 
tion from (—%», +) to 0 +~]. Of the 230 space groups, only CH: D3 “he Ol, 


8 
and On 


are excluded in this method of finding the formulas. The method can also 


be applied to the formulas for the structure amplitude for crystals in the lower 


systems. 


Introduction 


The International Tables [1] give formulas for 
the structure amplitudes of all 230 space groups, 
but these are not the working formulas, for they re- 
quire rather laborious manipulation to give the 
practical forms. Lonsdale's [2] simplified formu- 
las are much more convenient, for these take ac- 
count of the changes and simplifications arising 
from multiple values and combinations of the hk . 
Lonsdale's book also gives formulas for the elec- 
tron density, but these are only in a general form 
not convenient for practical use. 

Lonsdale's formulas are reproduced in the new 
English edition of the International Tables [3]; two 
settings of the unit cell differing in choice of origin 
are given for 23 of the 230 space groups, whereas 
it is very common to encounter different choices of 
coordinate system in practical analysis, not merely 
for these groups, either. The most pronounced 
changes in the formulas occur when the origin is 
transferred to a point characterized by another set 
of symmetry elements. Belov [4] has shown (for 
groups Dj,, Dh}, Ds,, Djn, etc.), that transfer of 


22 


the origin from an inversion center to certain other 
positions leads to a marked simplification in the 
formulas. 

There are very few copies of Lonsdale's tables 
and of the new edition of the International Tables in 
the USSR, and there is no corresponding work in 
Russian, if we ignore the obviously out-of-date 
Handbook on X-Ray Analysis [5], which was pub- 
lished in 1940. 


The formulas are also much altered if the di- 
rections of the coordinate axes are changed. The 
choice of coordinate system is often related to as- 
pects of the habit for various reasons of chemical 
crystallography, etc., so the setting of the crystal 
very often differs from the generally accepted one. 
This is especially so for monoclinic crystals; the 
space groups have at times to be considered in their 
body-centered, face-centered, and even stranger 
aspects. 


This means that the practical worker has usual- 
ly to deduce the electron-density formulas separate- 
ly for each particular case on the basis of the par- 
ticular setting of the crystal. 
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The strips method is presently the most com- 
mon in the summation of Fourier series; most of 
the mechanized techniques of electron-density com- 
putation employ this method, as does one of the 
techniques for the structure amplitude [6,7]. How- 
ever, it implies that the two- or three-dimensional 
distribution is considered as a set of parallel one- 
dimensional ones, so it is inevitable that the com- 
putation takes account only of those symmetry ele- 
ments that correspond to such a distribution, name- 
ly symmetry planes perpendicular to the coordinate 
axes, diad symmetry axes parallel to the coordi- 
nate axes, and inversion centers. 

The working formulas for the electron density 
thus involve referring the actual space group of the 
crystal to a subgroup belonging to one of the lower 
systems; only the latter require detailed considera- 
tion. 

The usual derivation of the electron-density 
formulas is based on the relations between the Fy 
(Fourier coefficients) for the various signs of the 
indices with subsequent transformation of the trigo- 
nometric functions of the formula.A much shorter 
derivation may be based on the symmetry and anti- 
symmetry of the trigonometric functions; this read- 
ily gives the formulas corresponding to particular 
symmetry élements in different combinations 
simply by elimination of trigonometric terms con- 
trary to the required symmetry. The relations be- 
tween the Fourier coefficients are used in explicit 
form only in the final stage, during conversionfrom 
summation between limits —~ —~ + to the same be- 
tween limits 0 ~~. 

The same method can be used to derive formu- 
las for the structure amplitude for all space groups 
in the lower systems, on account of the analogy be- 
tween the electron density and the structure ampli- 
tude. 


Symmetry and Antisymmetry of 
Trigonometric Functions 


We may consider the terms in the Fourier ex- 
pansion for the electron density as density waves in 
three-dimensional space, so we may consider the 
symmetry of the terms separately. The density 
wave represented by cos 27hx has m symmetry 
planes perpendicular to the X axis and passing 
through the origin and the middle of the cell edge a 
(Fig. la, b). Conversely, sin27hx has a plane of 
antisymmetry w, i.e., it transforms the density at 
point xyz into one equal in magnitude but opposite in 
sign at point xyz (Fig. lc, d). Further, cos 27hx 
has an m plane at x = ‘4 when h is even and an w 


Sin 2Mhzr 


COS 2MhI&L 


Fig. 1. Symmetry and antisymmetry elements of cos 2mhx and 
sin 2mhx for h = 2n and h = 2n—1., 


plane when h is odd. Conversely, sin 27hx has an 
ut plane at x = ‘4 when h is even and an m plane 
when his odd. These properties have been used [8, 
9] to shorten the strips. 

A further aspect of the symmetry is that the se- 
cond half of the cycle repeats the first when h is 
even; the repeat distance is half the cell period. We 
denote a translation of this magnitude by t and call 
it a half translation. The second half of the cell re- 
produces the first with the sign changed if h is odd 
for cos 27hx as well as sin27hx; this may be called 
an anti-half translation and be denoted by 1 . Table1 
summarizes these features of cos 27hx and sin 27hx. 

These properties are sufficient for the discus- 
sion of all space groups in the triclinic, monoclinic, 
and orthorhombic systems apart from Ci? = Fdd 
and D}, = Fddd, which contain d glide planes; inad- 
dition, the setting of the coordinate system causes 
various symmetry elements ‘to pass through points 
displaced from the origin by 4 of a period. These 
two features force us to take account of additional 
symmetry features of the trigonometric functions 
when h is even; but it is then laborious to derive 
the formulas via the properties of the functions, 
and the device appears extremely artificial. These 
space groups are therefore excluded from the gen- 
eral discussion. 

Other excluded space groups are Th = Fd3, Of 
= Fd3m, and O08 = Fd3c, which reduce to Fddd! 
after removal of the higher-order symmetry ele- 
ments and their equivalents. The other space 
groups of the tetragonal and cubic systems have d 
glide planes, as orthorhombic subgroups contain 
combinations of symmetry elements that do not in- 
clude such planes. For instance, groups pia 
= 142d, Té = 143d, and O}} = Ia3d reduce to 1242424; 
group Cie = 14md reduces to Imm; group Di, 
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TABLE 1. Symmetry Properties of cos 2mhx and sin 27hx 


h even h odd 
Function Rae sh Metis VAM Nek x=Y¥, | Ax='o 
cos 2rhe my (x) may, (x) | t (x) |mo (x) Ur), (@) | gle) 
sin 2rhx uly (x) my, (#2) | b(m) Jug (e) } my, (~) | a(x) 


= 14/amd reduces to Imma; group C}%, = I4cd re- 
duces to Iba, and group Di}, = ]4/acd reduces to 
Ibca. These groups, and all the others, may thus 
be discussed solely on the basis of the symmetry 
and antisymmetry features dealt with above. 
However, it may be that for those of these 
groups that have inversion centers it will prove con- 
venient or necessary to preserve the inversion cen- 
ter at any point having one or more components 
equal to 44 (or %, %, or %). The corresponding 
formulas then clearly cannot be deduced on the 
basis of the above properties; this applies to groups 
Ds2}, Dif» Dh, Th, and O. The present method 
of deriving the electron-density formulas thus ap- 
plies fully to 220 of the 230 groups and, with cer- 
tain restrictions, to 5 more groups. Only five 
groups are ruled out altogether: ete [ets ce Of, 
and O}- 
Superposition of the Operations m, w, 
iin SyGl a 


These elements of generalized symmetry are 
sufficient to represent all the other symmetry oper- 
ations encountered in the space groups of the lower 


COs 


systems. We may combine the functions Zit, 


2chz, 


oe Qnky, a 2xIz, (and hence these elements) to ob- 
tain a function having any of the possible sets of 
symmetry elements. 

It is obviously simplest to construct a function 
having two or three mutually perpendicular sym- 
metry planes [functions of symmetry m(x)m(y) or 
m(x)m(y)m(z)]. The other elements (rotation and 
screw axes, glide-reflection planes, inversion cen- 
ters) may be formed by superposing these simple 
combinations. For example, a diad rotation axis 
parallel to the Z axis may be formed by superpos- 
ing a function m(x)m(y) (which thus has a symmetry 
axis along the line of intersection of the two planes) 
and a function w(x)ut(y) (and thus a diad axis along 
the line of intersection): 


m (x) + m(y) + w(x) - w (y). (1) 


This leaves only the common element, namely, the 
diad axis parallel to Z. An inversion center may 


be formed similarly by superposing four functions 
having inversion: 


m (a) m (y) m (2) + m (a) we (y) wm (2) 
-+ we (x) m (y) we (2) + wm (a) wm (y) m (2). (2) 


All other symmetry and antisymmetry elements 
present in each function individually are eliminated 
by the superposition. 

The situation is rather more complicated with 
glide-reflection planes and screw axes. 

We assume that the glide-reflection plane is 
perpendicular to the X axis and has its glide along 
the Z axis, which means that the function should 
not change in sign or magnitude after two succes- 
sive operations (reflection at a plane perpendicular 
to the X axis and displacement through half a trans- 
lation along the Z axis). 

This is satisfied by combinations of trigono- 
metric functions of symmetries m(x)t(z) and 
wi (x) 1(z); but each alone has excess symmetry ele- 
ments: the first has a symmetry plane and a half 
translation t, while the second has an antisymmetry 
plane and an anti-half translation }. Superposition 


m (x) t (z) + w (x)2 (2) (3) 


removes these unwanted elements and leaves only 
the glide-reflection. Transfer by half a translation 
along the YZ diagonal may be considered as the re- 
sult of two successive operations t(y) and t(z) or 
xy) and 4(z); similarly, an anti-half translation 
along the YZ diagonal may be considered as the 
combinations t(y) - 3(z) or 4(y) - t(z), so a diagonal 
glide-reflection plane with a half glide is repre- 
sented by 


m (x) [t(y) t(z) +2 (y) + 2(z)) + w (a) [t(y)2 (2) +2 (y) t (2). (4) 


A 2; screw axis parallel to the Z axis may be 
represented as a superposition of a function of sym- 
metry 2(z)t(z) and one of symmetry Z(z)1(z), where 
zis a diad antiaxis: 


[mm (x)-m (y) + w (x) we (y)] t (2) 
+ [mm (x) w (y) + w (x) m (y)} 2 (2). (5) 
All of these relations are used below to derive 


the working formulas for the electron density in 
terms of Fourier series. 


General Working Formula for the 
Three-Dimensional Electron- Density 
Distribution and General Formula for 
the Structure Amplitude 


In the general case, 


p (wy2) = — >) Fa (hkl) cos 2m (ha + ky + Iz) + 
hkl 
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+23) Fp(hkl)sin2n (hat ky +h), 6) 
hkl 


where Fa (hkl) and Fp(hkl) are, respectively, the 
real and imaginary parts of the structure ampli- 
tude.! We transform the cos and sin sums via the 
usual trigonometric formulas to get: 


fos) 


o (xyz) = Sere G, (hkl) 
ae 


X [cos 2xhz cos 2nky cos 2xlz—cos 2xhzx sin 2ky sin 2nlz 


— sin 2rhz cos 2xky sin 2nlz — sin 2rhz sin 2xky cos 2nlz] 


es Ee >; Fx (hkl) [— sin 2xhz sin 2xky sin 2nlz 


hkl 


—oco 


+ sin 2xhx cos 2xky cos 2alz 


+ cos 2xhz sin 2xky cos 2n/z + cos 2nhz cos 2rky sin 2n/z]. 


For brevity we denote cos by C and sin by S, with 
the triple products representing in the first posi- 
tion h with x, the second k with y, and the third l 
with z. Then 


o (xyz) = vz F gibkl) WEC — CSS — SCS — SSC] 


BE Pik) 885 1 SGC GSC +. CGS)... \Ga) 


hkl 


—O 


It remains to complete the conversion to summation 
from 0 to ~ for each of the indices; Fa (hkl) 
= Fa (hkl) and Fp (hkl) = —Fp (hkl) , from which we 
get the relations of Table 2 [formula (6b)]. These 
relations are basic to the derivation of detailed 
working formulas for the electron density that ap- 
ply in the presence of the various symmetry ele- 
ments. 

The structure amplitude in the general case is 
given by 


F (hkl) = y fj cos 2x (hx + ky + 12) 


j=1 


N 
+ iD, f; Sin 2x (hx + ky + Iz), (7) 


Ta 


and so may be put in a form analogous to (6a): 


N 
=>) fICCE CS S- SCS SSC, 


j=l 


F (hkl)= 


N 


Li hl Sos ESCO CSC+ CUS 


j= 


(7a) 


Electron-Density Formulas for Single 


Symmetry Elements 


The presence of a symmetry element gives rise 
to relations between the structure amplitudes for 
different reflections, which enables us to simplify 
and manipulate the Fourier expansions for the elec- 
tron density. Each of the eight terms in (6b) bears 
a certain combination of m, w, t, and 1; summa- 
tion is equivalent to superposition of these com- 
binations, so only the symmetry operations thatare 
in common are retained. To each real symmetry 
element in the structure there corresponds a cer- 
tain set of terms in the sum of (6b); in deriving the 
working formulas, it is sufficient to note that of the 
eight terms listed in Table 2, we must retain only 
those that are in accordance with the symmetry of 
the structure. 

Consider the effects of single symmetry ele- 
ments on (6b). 

1. Mirror plane perpendicular to the X axis. 

a) The m plane passes through the origin; sym- 
metry condition m)(x). The formula for p (xyz) 
must retain only the terms that contain cos 27hx as 
a factor, so 


oO 


= Di(Meccece AnssCSS)] 


hki 
0 


0 (“yz) = 


+ [BoscCSC + Bacs€ES)}}. (8a) 

b) The m plane intersects the X axis at x = w/e 
symmetry condition m,/,(x). This means that the 
trigonometric part of the formula must retain only 
those functions that contain cos 27hx with h even 
and sin 27hx with h odd: 


4 (wyz) =3 >) [Aeon CC — Acgs€SS] 


hkl 


0 
for n=2n 


+ [BescCSC + BecsCCS}} + 


la jlowance has been made in (6) for the requirement F *(hkD) = 
= F(hkl), i.e., that a function expanded as a Fourier series must be 
real, 
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TABLE 2 
ee  ____ ae 


p (yz) = 


iS Pe 
hkl 
if) 


ly 

eee ae 

Acog (hkl) = FF 4 (uk) + Fy (hkl) + 
+h (hKT) + Fy (hkl) 


(MkO) = 2 [FP , (akOY = Fy (ARO)| 
matey y 


x4 
Cann iii g 
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Awe 
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Meigs (ARO) = 7 By (UK) 
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A ggc (R01) = Ags: (Ohl) = 0 


Lose (OhO) = - Ags, ( 


= Aga, (000) = 0 


Ags, (00) = (01) = 


oo 


ao 
SM Agge CCC = fngg CSS hae SOT me Aged SC| +7 > (-Bsgs 555 + Bsce SCC + Bege€SC + Bees CCS! 


Kk 
vu 


Ip 
Bygg (hkl) = + [Fp (Akl) + Fp (hkl) + 
+ Fp (hkl) = Fp (hkl)] 


2, Eee A Bagg (Okl) =0 


B ges (il 10) = Bes ‘Ss (URO) = 
= Begs (000) = U 


Begs (WU) = — 


te 
i 
q 
— Fg (hkl) — Fp (hkl)| 


| 
Bary (RED) = > [Fg (akO) — F p (RKO) 


Beco (hkl) = = [F g (Akl) + Fp (hkl) — 


+ 

i te) 
Bee Ml) = ze [Fp (hOl) + Fp (h0l)| 
Bgpp (ORL) = 0 
Bye (HOO) = = Fp (hOW) 


j 
By 
(Oh0) = Beye (001) = B gee (000) =0 


Bee 


Hl 


Boge (kl) = [Fp (hkl) — Fp (hkl) + 


wa| —» 


Fp (hkl) — F p (hkl) 


j - = 
Boise (AKO) = 7 [Fp (RKO) -- Fp (hkO)| 
1 ‘ 4 hoo 
Boge (OK) = 7% [Fp (ORL) — F p (Ok/)| 
ie: (AVL) = 0 
Boge (Ok) = z Fp (OKU) 
Bose (hOV) = SC (002) = Bese (000) =0 
IVp 
| 


Bogs (hkl) = = [Fg (hkl) — Fp (hkl) — 


4 


— Fy (Akl) + Fy (hkl)] 


1 Bs 
Boos (h0l) = = (Fy (hOl) — Fp (h01)| 
1 = 
Becs (Okl) = 7 [Fg (Ohl) + Fp (Ok1)| 
Bees (hkO) = 0 
1 


Boog (00l) = =F p (002) 
Becg (200) = Boog (0K) = Begg (000) =0 
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ee ory ee 
hkl 
for aot 


+ [—BgssSSS + BsocSCC}}. (8b) 


2. Glide-reflection plane perpendicular to the 
X axis with the glide along the Z axis. 

a) The plane passes through the origin; sym- 
metry condition m,(x)t(z) + W(x) 1(z). Either of the 
functions cos 27lz and sin27lz has symmetry t(z) 
for 1 = 2n and 1(z) for 1 = 2n—1, so the formula 
must contain only terms in cos 27hx for 1 = 2n and 
in sin 27hx for 1 = 2n —1; for l = 2n we retain 
terms CCC. Css. CSC. and CCS. while for l= 2n 
—1we retain SCS, SSC, SSS, and SCC. 

b) The plane intersects the X axis at x = Je 
symmetry condition m4 /4 (x)t(z) +uy/4(x)4(z). For 
hes /-—)onawe retain CCC, CSs. CSC; and CES; 
while for h + 1 = 2n —1 we retain SCS, SSC, SSS, 
and SCC. 

3. Glide-reflection plane perpendicular to the 
X axis with the glide along the YZ diagonal. 

a) Then plane passes through the origin; sym- 
metry condition 


Tidy (%) (L(y) £(2) + 2 (y) 7 (2)] FH Wo (2) EE(Y) 2 (2) + 7 (y) €()]- 


The formula must contain cos 27hx for k + 1 = 2n 
and sin 27hx for k + 1 = 2n —1, so the terms in CCC, 
CSS, CSC, and CCS are retained only for k + 1 =2n, 
and those in SCS, SSC, SSS,and SCC for k+ 1 = 2n 
cad be 

b) The plane intersects the X axis at x = wh 
symmetry condition 


ma, (2) (£(y) £(2) + 2 (y) 2 (2)] mes, (2) [4 (y) 2 (2) + 2 (y) t(2))- 


The formula contains CCC, CSS, CSC, and CCSfor 
h +k +1= 2n and SCS, SSC, SSS, and SCC for h+k 
ba= 2tie—eles 

The various cases of a symmetry plane of a 
given orientation show that the electron-density 
formulas are extremely similar; the only differ- 
ences lie in the conditions on the indices of the re- 
flections (Table 3). 

The trigonometric terms in the formulas split 
up into two groups. A plane perpendicular to the X 
axis causes the functions containing cos 27hx to be 
active for one combination of the indices and those 
containing sin 27hx for another. 

4. Diad rotation axis parallel to Z axis. 

a) The 2 axis passes through the origin; sym- 
metry condition my(x)mo(y) +U9(x)Wy(y), which is 


satisfied by functions containing as factors cos 2mhx 
- cos 2mky and sin27hx - sin27ky. Consequently, 


p(ty2) = 8 D {tAcce CCC — Asse SSC] 
hkl 
0 


+ [— Bsss SSS + BecsCCS}}. (9a) 
b) The 2 axis passes through the point x = pe 

y = 0; symmetry condition m, /4 (x)mg(y) + Uy /4 (X) 

‘ul (y), which is satisfied by functions containing 

as factors cos 27hx - cos 27ky and sin 27hx - si 

‘ sin 27ky for h even, and by ones in cos 27hx - 

* sin 27ky and sin27hx - cos27ky. Then, 


8 < 4 
G\(@22) — => {[Aece CCE — Agse SSC] 
hkl 


0 
for h=2n 
+ [— Bsss SSS + Becs CCS}} 


8 co 
if Te Sy {[— Acss CSS — A scs SCS] 
Akl 
for h=2n—1 
35 [Bscc SCC + Bese CSC}}. (9b) 
c) The 2 axis passes through the point x = ¥, 
y = ‘4; symmetry condition m4 /4(X)my /4(y) + 
+ UL; /4(x) Wy /4(y), Which gives the same functions 
but with different relations between the indices: 


for h+k = 2n, functions CCC, SSC, SSS, and 
CCS; 

for h+k= 2n —1, functions CSS, SCS, SCC, 
and CSC. 

5. Screw axis parallel to the Z axis. 

a) The 2, axis passes through the origin; sym- 
metry condition 


[719 (%) Mg (Y) + Wy (X) Wy (y)] t (2) 
Ht [rg (©) Mey (Y) A Mg (X) Mo (Y)] 7 (2). 


This means that for 1 = 2n, i.e., t(z) symmetry, 
we have CCC, SSC, SSS, and CCS, while for 1=2n 
—1, i.e., for 3(z) symmetry, we have CSS, SCS, 
SCC, and CSC. 

b) The 2, axis passes through the pointx =, y = 0; 
the symmetry condition has m)(x) and Uh)(x) re- 
placed by m4/4 (x) and wy/4(x), so the formula con- 
tains the same combinations of the trigonometric 
functions, but for h+l = 2n andh+l1 = 2n —1 in- 
stead of h = 2n and h= 2n—-1. 
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TABLE 3, Formulas for the Three-Dimensional Electron-Density Distribution for Various Dispositions of Single Symmetry Elements 


See 


I. Symmetry Plane Perpendicular to X Axis 


Conditions on the indices 
i —* ! 
ian 
Functions present ited A(hkZ) and B(hkZ) m plane c plane nplane 
etween 
in formula in Fourier series a aT ative) OP 
F(hk ) 
yl) Neb 4 x -0 NasT ly azell Ket fy 
JOC, lege 2 a VP (ied) bP (DAD 
= OSS |, | F (HB) = Nigga Thy [P, HED - Fy (heb 
| ae =~ I (hil) ne ee 2] hki — afl h = 2n ] = 2n Af P= hb = 2e| ep ktle= 
JOASG, Bec = Va fh, (Uh) -| I), (hkl)| = nN =2n 
COS |i, 
Bog 2 Ya [By (kl) — Fy, (HAD) | | 
[= SCS, Agog = Va [M4 (aKI) -|- I 4 (hk 1)| 
| SSC] 4 Ase = as ar (hkl) — diay (hkl)| 
Ne ee) oa ee ne eg pera ‘1, (hkl) 
i |— SSS, +e Fogg = la [Pp (LA) + OF), (hk) | None h=2n—1 | J =2n—1[h+l= k+l= h+k+l= 
SCCl, = — (hkl) =2n— 1) =2i—1 = 2n—1 
: Byo go = Ms (Fp (hkl) — F,, (AAD) 
Il, Symmetry Plane Parallel to Z Axis 
Conditions on the indices 
Relation ‘ 
Functions present Laeien A(hk?) and B(hk?) 2 rotation axis 2, screw axis 
in formula F(hk1) in Fourier series 


GGG Accg = Va Fg (AK) + Fg (kK DI 
— SSC , pelts 
I ESCRRE) aN ia LE (hkl) =F gid] 
| Es = Mkt) _. | hkt — all h=2n htk= l=2n he 2a hp kh l= 
|— ss Boss = Va Ey (kl) | By (lke | ssh =2n 
COST, 
Boog = Vo [Fy (hil) — Fy, (hed) 
[= CSS, Aogs = Va U4 (hkl) + F4 (hk d| 
— SCS], ae 
} ‘gag = Ye [Pg (akl) — Bg (hk1)| 
“yey I Al = , 7 oe 
I SCC, Ok) = | Boog = Ma Fp (hkl) + Fp (tkD| | None | h=2n—4] athe |l=mmot) htie= |htkti= 
CSC], |= — ¥ (hkl) =2n—1 =2n—1 | =2n—1 
Boso = Ve [Fp (hkl) — Fp (he) 


*In addition, F(hk/) = F *(hk7), i.e., Fa(hkl) = Fa(hk?) and Fp(hkl) = -F,(hk?). 
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Ill. Inversion Center 


Relation between 
F(hk2) 


Functions present 
in formula 


F (hkl) = F(hkl) | &C© 


{ 


A(hkl) and B(hk/) in Fourier series et Beaty 


=", [F4 (hkl) | Wy (hil) + Fy (hkl) + FA (hiD)| | 
SCS, —SSC| , Angs = Va [P (akT) | By (hkl) -- 2, (hid) — Fy (hkl) 
clgeg = Ma Lg (ah) = Fa (hkl) + By (kt) — By (hk d)| 


chase = Ya [Fy (AKL) — Fy (ahh) =F (bRL) EB y (hkD)| 


Conditions on the indices 


ne=al/y, 
w—0, y—0, 


waa" 4, 
W"/i, 
e—) ai) a—() z= "I4 


hkl — all 


| | 


[— SSS, SCC, |F (hkl)= —F (hk1)| 


I] 


c) The 2, axis passes through the point x = ey 
= vi the symmetry condition causes further change 
in the condition on the indices in (9b), h = 2n being 
replaced by h+k+J1 = 2n andh=2n-—-1byh+k+l 
=2n-1. 

A symmetry axis parallel to the Z axis thus 
causes functions containing cos 27hx - cos 27kx or 
sin 2Thx - sin 27ky to be active for one combination 
of the indices, and those containing cos 27hx 
- sin 27ky or sin27hx - cos 27ky for another. 

6. Inversion center. 

a) Inversion center at the origin; symmetry 
condition : 


reg (2) 17 (y) 12g (2) EM (tt) tly (y) tty (2) 
sf Uy (7) Mg (Y) Wey (Z) Ar Wl (FH) Wey (Y) Mo (2)- 


Working formula: 


8 oO 
o(zyz) =P ~ [Acce CCC 
hk | 
0) 


=A css (GSS — Ascs Cio) = Agsc SSC} 5 (i 0a) 


b) Inversion center at x = 4, y= 0, z= 0; the 
symmetry condition has m (x) and W(x) replaced 


by m4 /4 (x) and uy/,(x), and, therefore, cos 27hx be- 


comes sin 27hx when h = 2n — 1 and. conversely, so 


Bagg = Ma Py, (Uhl) +P) (HRD) 2), (hkl) fe By, (akD)| 
CSC, CCS), Ryo = Vall (ikl) | By (hkl) — Fy, (kl) — Fy, (aki) 
Boge = Val Py (Akl) — Fy, (hkl) +P), (hkl) -~ Fy (kd) 


Boog = Ma [ly (hkl) — By (lel) — Fp (ake) + Py, (aD) 


None h = 2n— | 


=a 


8 ao 

¢ | 

pea 
hkl 


[eee CCE 


) 
for h=2n 


=< A, SS CSS. a As Ss SCS = Asse SSC] 


foe} 
8N RECS 
+ 7 DI Bass SSS 
hikei 
0 
for h=2n—L 


+ Bee SCC + Bese CSC + Bees CES (LOb) 


c) Inversion center atx = 4, y= %, z= 0; the 
conditions on the indices become h +k = 2n for the 
first sum, and h+k= 2n —1 for the second. 

d) Inversion center atx = 4%, y= 4, z= y: 
the conditions on the indices become h+k+J/ = 2n 
for the first sum, and h+k+Jl = 2n —1 for the se- 
cond. 

An inversion center thus causes one group of 
functions to consist of combinations containing an 
even number of sines (or CCC alone) and the other 
group to contain an odd number of sines. 

Table 3 collects these results, which are pre- 
sented in separate parts for symmetry planes, sym- 
metry axes, and inversion centers. 

Each part has two horizontal sections: I, for 
the terms of the Fourier series that appear for the 
first (even) combination of the indices, and II, for 
those for the second (odd) combination. 


Nh f-h = 2n)h | kh4-1=2n 


firctee teT asx 
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Each section gives the functions that are re- 
tained and the conditions (combinations) for the in- 
dices for all types of symmetry element in all po- 
sitions. It also gives the relations between the 
structure amplitudes arising from the symmetry 
elements and the values of A (hkl) and B(hkl) for h 
=0, k~ 0, andl =~ 0 implied by these. 

The general formula of Table 2 may be used 
with the data of Table 3 to write in expanded form 
the electron-density formula for any symmetry ele- 
ment and any setting of the crystal. No particular 
difficulty arises in permutation of the letter sym- 
bols in the functions and indices that is necessitated 
by a setting differing from that of Table 3. 


Working Fourier Series for the 
Electron Density for the Various 
Space Groups 


Tables 2 and 3 form the basis of the derivation 
for the space groups considered here; we retain in 
(6a) only terms that do not conflict with any of the 
symmetry elements in the group, so we have to take 
from Table 3 the functions characterizing the indi- 
vidual symmetry elements (permuting C and S in 
accordance with the orientation of the element). 
These are compared to leave only the ones present 
in all. The A(hk/) and B(hkl) are simultaneously de- 
rived from the relations between the structure amp- 
litudes. 

The following are some examples. 

1. Space group D§, = Pnnn (Fig. 2). 

a) Origin at the inversion center, i.e., at the 
intersection of the three planes. 

The inversion center requires us to retain only 
[CCC, -CSS, -SCS, —SSC]a for all hkl; the n plane 
perpendicular to the X axis, only [CCC, -CSS],, 
[CSC, CCS]p for k + / = 2n, and [-SCS, -SSC]a, 
[-SSS, SCC]p for k +7 = 2n —1; and the diad axis 
parallel to the Z axis and passing through x = 4, y 
= /, requires us to retain (CCC, -SSC]a, [-SSS, . 
CCS]p for h +k = 2n and [-CSS, -SCS], ,[SCC, CSC] p 
forh+k=2n-—-—1. All three conditions are satis- 


fied by 
CCC for h+-h=2n  andk +1 = 2n (a), 
—CSS.~ for 2+ i= 2n—! andhk + 1 = 2n (bd), 
—SCS for h+k& = 2n—landk +1] = 2n—! (0), 
—SSC for ht+A=2n andk+l=2n—l1 (a), 
so 


oO 

8N ; F 

6\ (7/2), = Tv b> elec cos 2nhZ cos 2aky cos 2nlz — 
A‘L 
O(a 


) 


oo 
sy . . 
sarah BA Acss¢0s 2nhe sin 2rky sin 2nlz 
hkl 
(4) 


ios) 
8 as ‘ . 
Say dA scs sin 2xhe cos 2nky sin 2alz 
hkl 
Ofc) 


oo 
8N , ; ‘ 
a >, Agsgc sin 2xhe sin 2xky cos 2nlz, 
Akl 
O(d) 


in which (a), (b), (c), and (d) denote the above con- 
ditions imposed on the indices. 

b) Origin at the intersection of the diad axes. 

The inversion center at 4, 4, 44 leaves [CCC, 
—CSS, -SCS, -SSC], forh+k+/ =2n and [-SSS, 
SCC, CSC, CCS]p for h+k+ 1 = 2n —1; the plane 
perpendicular to the X axis and passing through x 
= ¥, removes half of each set of functions, leaving 
[CCC, -CSS], for h+k+J = 2n and [-SSS, SCC]p 
forh+k+l =2n-—1, while the diad axis coinci- 
dent with the Z axis necessitates [CCC, —SSC], and 
[-SSS, CCS]p for all hki. Then all elements are 
satisfied by 


picyy=5 D> 


ARK A=0 
(h--k+1=2n) 


Acer cos 2nhz cos 2xky cos 2alz 


Bgsg sin 2ahz sin 2exky sin 2elz. 


2. Space group C$, = C2/c (Fig. 3). 

a) Origin at the inversion center. The glide- 
reflection plane perpendicular to the Y axis passes 
through the origin and has its glide in the direction 
of the Z axis, so a permutation is necessary in the 
formulas of Table 3 (Y is replaced by X). Thec 
plane leaves [CCC, -SCS],, [SCC, CCS]p for 1 = 
= 2n and [-CSS, -SSC]a, [-SSS, CSC]p for 1=2n-1. 

The inversion center lies at the origin, so the 
formula can retain only functions in the brackets 


} \ { 


Fig. 2, Space group Den = Poon, 
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Fig. 3. Space group C$, = C2/c. 


with subscript A. The centering of the cell imposes 
a condition on all reflections present: h +k must be 
even. These three elements are sufficient to gen- 
erate the others, and hence the conditions are suf- 
ficient. The working formula is 


co 
8 
6 (zyz) = - > [4ece cos 2nhz cos Achy cos 2alz 
t 
h=2n, (=2n 
— Ases sin 2nhe cos 2xky sin 2alz] 
co 
8 I ? Pale . 
anon >; [— Acss cos 2ahe sin 2aky sin 2alz 
hkl 
0 
h=2n—1, |=2n—1 


-— Asso sin 2nie sin 2xky cos 2alz). 


b) Origin at the intersection of the rotation 
axis with the c glide plane. The latter necessitates 
the same functions subject to the same conditions. 
The inversion center is atx =0, y= 0, z= 4, soit 
leaves group A functions for J = 2n and group B for 
1 =2n —1. Then, 


oO 
8 A 
0 (xyz) = + > [Acco cos Anhz cos 2rky cos 2xlz 
hkl 
i) 
h=2n, l=2n 


— Ages sin 2nhe cos 2xky sin 212] 


co 

+7 > [— Hssg sin 2xhx sin 2rky sin 2nlz 
hk 1 
0 

h=2n—1, l=2n—1 


+ Besc cos 2xhs sin 2ky cos 2x/z). 


It is equally simple to find the expansion when 
the origin is placed at the intersection of the screw 
axis with the c plane or the n plane, or at another 
inversion center. 

3. Space group ne = Ja3. The orthorhombic 
subgroup of group Ia3 is D*j, = Ieab (Fig. 4). It is 
most convenient to take as basis the c plane per- 
pendicular to the X axis, the 2, axis parallel to the 
Z axis, and the inversion center. 

We assume that the origin is placed at the in- 
version center, the c plane passing through the ori- 
gin and the 2, axis being displaced by 1), repeat dis- 
tance along the X axis. Table 3 gives directly that 
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I 


oO 

S 

Oo . ; . 5 

a (ry2) = ca » Be Vey, co COS Qe COS Lnky cos 2zlz 

hkl 
0 

h=2n, l=2n 

ao 
8 oN 5 : 
— Acss cos 2nhyr sin 2ahky sin 2alz 


( 

h=2N—1, l-=2n 

lao} 

8 o wae 

a Ages sin 2n/ur cos 2rhky sin 2alz 
Ik | 


0 
i=2n, (=2n —1 
ioe} 
8 ae eee 
aera alsse Sin 2ahe sin Zrky cos 2xdz. 
hkl 


0 
h=2n—-1, l=2n—-1 


There is no need to examine the transforma- 
tion of A(hkl) and B(hkl) for each group separately. 

No relation between the F(hkl), F(hkl), and 
F(hkl), etc., can reduce any of the A (hkl) or B(hkl) 
to zero, because this would-imply loss of the cor- 
responding trigonometric function and hence change 
in the symmetry of the distribution. Then equality 
in the absolute values of two amplitudes, e.g., 
|F (hkl) and |F(nkl)|, can lead only to simplification 
(not loss) of the corresponding coefficients: 
Accchkl) = ¥.[F hkl) + F(hkl)] will be replaced by 
Accc(hkl) = F(nkl) or Acgg(hkl) = 4[F(hkl) — F(okl)) 
by Acgs = F(hkl) [the F(hkl) reflection is obliged to 
have the sign corresponding to this change in the 
formula]. 

On the other hand, we know that the amplitudes 
for the eight reflections that differ in the signs of 
the indices always split up into two independent sets 
of four that differ in absolute value: 


| F (hkl) | = | F (hkl) | = | F (hkl) | = | F (had) | 
and 
| F (hkl) | = | F (hkl) | = | F (hkl) | = | F (hkl)| 


(we assume that the special axis is the Y axis ofthe 
monoclinic crystal). Similarly, the absolute values 
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TABLE 4. Coefficients for Fourier Expansion of the Electron Density 


I. Monoclinic System, Y Axis the Special Direction 


nnn ee EUUE ENED SEnEEEREESSSRRRERREIET SEERERERREEREEEEnened 


Indices ACEC Acss Asus Assi Risss Esco Rese Bocs 
I eas | sae I I t 
hkl x lPy (hkL) | lk (hkl) “| ly (Uh) ly (hbly— | x lF p(hkl)—— | IF p (hkl) — | sl (hkl) + | a [Fp (hhl) + 
| By (ikl)| | Pa (b)| | By (al) ~— Py (hhi)| — Fy, (hil)| Py, (ohl)| +P 4, (hbI)| + Fy, (at) 
es } Te AiO) BS (hk0) 0 
hhO x Pa (hb) 0 0 sy Fg (KO) 9 ny Pia (AA ZR 
f 
} te tp — [Fp (AOI 
hOl Tl, (HOI) 4 0 ey [Ay (hO1) = 0 0 z [FP (Ol) — 0 j [Fp (AO?) + 
++ F 4 (h0l)| — Fg (h0l)] — F p (hOL)} + F ,, (h0l)| 
: ee I 
ai a F 4 (Od) ik (Okt) 0 0 0 0 x Fp (Ohl) = F p (Okl) 
| 
; I 
hO0O z Fy (00) i) 0 0 0 . F ;, (h00) 0 0 
We ; ; 1 oe OE 
OKO % F , (OKO) 0 0 0 0 0 ra p (OKO) 0 
| RS 1 
O0L Te F , (001) 0 0 0 0 0 0 4 F (002) 
duet 
000 oF 4/008) 0 0 6 0 0 0 0 
Il. Monoclinic System, Z Axis the Special Direction 
Indices SCCC Noss Ases Assc Bsss | Boe | Bose Lees 
bhi ! ere ; fies [ae Lice big ads sear (ngbess 
h ay Mg (ad) sy Hy (URE) FE) Ug (RL) ap Tg ad) == | Ly ha) =P Pg | oy Fg kd) | ap LF (kd) + 
4 Fy (hki)| 4, (hkI)| — Fy (hkl)| —- Fy (hkl)| — Fy, (hkl)| — I, (hkl) + Fy, (ike)| + Fy (kel) 
(ee (ie ie ! i she 
hkO GE g (ek) + vu 0 GE ag (hu) ~ 0 ZOUF y (hav) -- - [Fy (hk) 4- 0 
+ iy (hk4))] —Fy (hkO)] BEN (hkO)] ae (hkO)] 
| 
ho meee CHAN, i 5 I 
ee ps ON) 0 x Pa (RO) ) 0 ap F yy (L02) 0 Ls F , (AOL) 
tien 
Oke — KF, (OkL) es (URL 1) 0 1) 0) one i 
x a y Fa (kl) + Fy (Uke) Z F 1 (Okt) 
hOU : Lf, (hOO () 0 ee 
WW 7 A (LOO) 0 0) G fy (00) 0 0) 
ile 
OKO si Py (OKO) Uv U) 0 0 i} zits (OkV) 7) 
A ; I 
O01 2 F 4 (000) 0 0 0 ) 0 1) F F p (002) 
fiy%: 
000 ak ‘000) 0 0 0 ) 0 0 0 
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Ill. 
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Orthorhombic System 


Indices Accc Acss Asgcs Asse Rgss Bsce Bese Becs 
hkl F , (hkl) Fy (hkl) Fy (hkl) Py (hkl) F (hkl) BF, (hkl) F p (hkl) F p (hkl) 
| | 
ial ee) 0 0 <P, (hkO) 0 +. Py (ik0) < Fy (k0) 0 
hOI a (01) ) I i, ae 
a ro + A { ole (01) 0 () J Fy, (hOL) 0 ap (102) 
{ 
; -— fF Ak] ; ; ' n | { 
Ohl eA x Fa (OL) 0 0 la 0 = Fp (Uhl) ay Fp (Ok1) 
2 
| Lp gnc } , | 
se aa , 0 0 re Z F (200) 0 0 
| 
f | 7 anya b 
OkD = I 4 (ORO) 0 0 0) t i) O Fp (Oh) 0 
| 
| I 
{ [i A 
OU! See (O0/) 0 0 0) | 0 () v ee p (O0L) 
Hi 
{ 
(00 3 F , (000) 0 () () | ( 0 6 0 


of all eight amplitudes are identical for any group 
in the orthorhombic system. 

The transformed expressions for Accc, ACSS, 
etc., will thus be the same for all monoclinic 
groups; the same is true for the orthorhombic 
groups, and so these can be written down in the 
necessary form in advance (Table 4). 

The space group of a crystal of higher sym- 
metry is referred for the purpose of electron-den- 
sity calculation to the highest subgroup in one of the 
lower systems, so the relations of Table 4 apply for 
any symmetry group. The higher groups that are 
referred to monoclinic ones are usually oriented 
with the special axis along the Z axis (not Y), so 
both possible settings are considered. 


Working Formulas for the Structure 
Amplitude for the Lower Systems 


The general series formulas of (6a) and (7a) 
are very similar, for they contain the same com- 
binations of functions with the same signs. The 
functions present in the real part of F(hk/) coincide 
with those multiplied by the A coefficients in p (xyz), 
while those in the imaginary part of F(hkl) coincide 
with those multiplied by the B in p (xyz). The same 
combinations are lost in the two cases when a sym- 
metry element is introduced. 


The loss of a function from p(xyz) occurs be- 
cause the symmetry of the function must satisfy the 
requirements on the symmetry of the crystal. The 
F(hkl) lose functions on account of summation over 
atoms related by symmetry elements; but here the 
only functions that can persist are those whose sym- 
metry coincides with that of the array of atoms, 
otherwise the repeated functions would have (for the 
given fj); opposite signs and so would cancel out. 

Thus, the formula for p(xyz) as a Fourier 
series, must resemble that for the structure amp- 
litude for all space groups; Table 3 may be usedfor 
both. 

This gives a somewhat curious air to Nowacki's 
book [10]on the derivation of electron-density for- 
mulas on the basis of those for the structure ampli- 
tude, For instance, for group Dj = P2,2,2,, that 
derivation takes about four pages, whereas the final 
result can be written down directly from the analogy 
with the initial formula. 


Plane Symmetry Groups 


The present approach can also be applied to 
two-dimensional cases (projections and sections); 
difficulty in discarding the unwanted dimension can 
be avoided by drawing up separate tables for the 
general formula (12) (Table 5) and for the combina- 
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TABLE 5. General Formula for a Two-Dimensional Electron-Density Distribution 


airy) Ne. Cos 455581+ 4 1668 +B, SC] 
~ 8 He aK 
0 
ly I, 
A,-c (ik) =, [PF 4 (hh) -- Fy chk) | Bg (hk) =o (F p(hk) + Fp (hk)] 
Avg: (RO) = oF 4 (RY) By (hO) = 0 
Via (Oh) = 1 oF , (Oh) Bug (Ok) = 1 oF p(Ok) 
Fy) eae cr LID) Byg (00) =0 
a MW 
Aloe (eh) == 2 fy (Ah) — Fh 4 thk)| Bo (hk) = 35 [F p (hk) — Fy (hk) | 
A sg (h0) = 0 Bg (RO) = oF p (h0) 
Ag (OK) =O. Bay (Ok) = 0 
Ags (OO) = 0 | Bs (00) =0 


TABLE 6, Formulas for the Two-Dimensional Distribution for Various Dispositions 


of Symmetry Elements 


I. Line of Symmetry Perpendicular to the X Axis 


Conditions on the indices 


Functions | Relation rr ; ii 
present between A(hk) and B(hk) ees gone 
in formula | F(hk) Le apiec'O es a0 Realy 
a J 
[} (CCly | FUR) = [Ape (hk) = Fy (hk) oak | k= Qn | ha 2n |hek= 
= F (hk) =n 
iGsiy * Bog (hk) = Fy (hk)| AML 
ll h + 
55a. BUR) | Ane Gms Fy UR) sialon 
Pegi 5) None |h=2n—//kA=2n—1] =2n—1 
iSO * Bop (hh) = F p (hk) 


Il. Inversion Center 


——— ee ee 


I] 


Functions | Relation Conditions on the indices 
present between A(hk) and B(hk) 
in formula F(hk) as eon peels 
1 
| 
eg (hk) = 12 [Fg (hk) + 
(CC, = SS)eq F(hk)= Fk), - F , (hk)] h, k 
Ags (hk) = 3.5 [F 4 (hk) — h=2n jhik= 
| — F ek) Al =n 
Beg (hk) =o [F p (hk) + 
(CS, SC], | F (hk) = ~~ Fp (hk)) 
=—F(hk)| Bee (hk) =o [Fp (hk) — | None | k= |atk= 
(hk) n—1 n— 1 


*In addition, F(hk) = F*(hk), ive., F (hk) = Fy (hk) and F, (hk) = -F (hk). 
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The structure of the addition compounds formed by two amines with the halides 


of divalent cobalt is considered. 


which proves molecular, is determined. 


balt is 4, the coordination polyhedron a tetrahedron. 


The crystal structure of CoCl, ° 2C,HoN, 


The coordination number for the co- 
The authors confirm that 


the a- and 6-modifications of compounds of divalent cobalt (blue and violet) 


are not cis and trans isomers. 


The first has molecular tetrahedral structure, 


while a chain structure with coordination number 6 for the cobalt must be as- 


cribed to the second. 


The structure of the addition compounds of 
general formula CoX, - 2A, formed by two am- 
monia molecules, or two molecules of an organic 
amine, with halides of divalent cobalt has attracted 
the attention of many authors. Biltz and Fetken- 
heuer [1] found that the diamine of cobalt chloride, 
CoCl, - 2NH3, existed in twoforms. A similar ob- 
servation was made for the dipyridinate CoCl, - 2 Py 
{2]. Various views have been published regarding 
the structure of these compounds. Biltz [1] and 
Hantzch [2] thought they dealt with cases of cis - 
trans isomers, as in the compounds of divalent 
platinum. Cox [8] used chemical and crystallo- 
graphic methods to study the two forms of the di- 
pyridinate of cobalt chloride: the stable, violet, a 
form and the unstable, blue, 8 form. Considering 
the dimensions of the unit cell in the violet form, 
he came to the conclusion that the lattice contained 
plane CoCl, - 2Py molecules with the pyridine 
groups in the trans position. From the magnetic 
properties of the two forms of CoCl, - 2Py, Mellor 
concluded [4] that the violet form had a chain octa- 
hedral structure and the blue a tetrahedral struc- 
ture. X-ray studies of Cs3;CoCl, [5,6] and x-ray 
[7] and electron-diffraction [6] studies of Cs,CoCl, 
showed that these compounds contained isolated 
CoCl, tetrahedra. Electron-diffraction examina- 
tion of the dihydrate CoCl, - 2H,O indicated that 
this compound had a chain octahedral structure [8]. 
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Since the structure of the diamine dihalides of 
cobalt thus seems to be in doubt, we undertook an 
x-ray investigation of the addition compound of two 
molecules of p-toluidene with cobalt chloride, 
CoCl, - 2H,NCgH,CHs3. 

Quite large single crystals of this compound 
are obtained on evaporating a mixture of solutions 
of anhydrous cobalt chloride and p-toluidene in 
ethyl alcohol. Rose crystals of the dialcoholate 
CoCl - 2C,H,N - 2C,H;OH were found at the bottom 
of the beaker, while at the sides there was a gradu- 
al growth of blue plates of CoCl, - 2C;H.N. 

Found%: Co 17.21; Cl 20.67. 

Calculated %: Co 17.12; Cl 20.60. 

The crystals form concretions parallel to the 
extension. From optical data, the crystals are bi- 
axial, the angle between the optic axes being 2v 
* 90°. The refractive indices are Ng = 1.701, Nm 
= 1.652, and Np =1.610. 

The dimensions of the unit cell were deter- 
mined from oscillation x-ray diffraction photo- 
graphs and then improved by photographs obtained 
on a KFOR x-ray goniometer. The dimensions 
were: a=12.3+ 0.05 A, b=4.594 0.01 A, c 
= 26.1 + 0.1 A, 6 = 93°45'+ 5", The density, de- 
termined pycnometrically, was 1.483, which gives 
for the number of formula units in the unit cell n 
= 3.83 * 4, The density determined by x rayswas 
1,55. 
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The Fedorov symmetry group was determined 
by indexing the x-ray goniometer photographs. The 
indices of the reflections satisfied the following 
conditions: 

hkl only withh +k+1 = 2n; 

hO0lZ only with h = 2n and 1 = 2n; 

0k0 only with k = 2n, 
which characterizes the Fedorov groups C$, and ene 

The structural investigation was carried outfrom 
goniometric photographs taken on the KFORwith Mo 
Ka radiation. The intensities were estimated by com- 
parison withablackeningscale. Only the polariza- 
tion and kinematic factors were taken into account. 
Some 556 independent reflections were obtained, 
taken from five layer lines (h0l, hil, h2/, h3l, h4d). 

From these data we constructed.a series of in- 
teratomic vectors and the projeetion of the electron 
density on the XOZ plane. The Fedorov group and 
atomic position were established by analyzing the 
three-dimensional F? series and the electron-den- 
sity projection. 

Maxima lying in the XOZ and X‘44Z planes indi- 
cate the presence of two-fold rotation and screw axes 
and thus unequivocally give the Fedorov group Coh 
= [2/a. 

The cobalt atom occupies a fourfold particular 
position with one degree of freedom, lying on the 
two-fold rotationaxes. The chlorine, nitrogen, and 


carbon atoms occupy an eightfold general position. 
The coordinates of the atoms are 


Coz oe i, OF oa y, 0; 


3 


I a ae! 
Tar eae 


Dates 


Et). : 7y = 03386. 


The general eightfold position for the Fedorov 
group Con in the I2/a arrangement is 


poe reg { 
Ryes Eyed, ey ee 


a] 1 1 ra 
tb Gy = Oly 1 2; yt, Y, Zz, 


aoe | { 
Die eens aa 


1 1 4 1 = 
zt, yy, y FA ay ee Y; z. 


Gla == 0738447 — 02052 ONS4s 
NE ae== (MISS 7 = OHS, 2 = 0.052, 
Ge t= 04067 == 0 ere — L096). 


The structure of CoCl, - 2p:-— HyNCgH,CHs3 is 
molecular. The cobalt atom lies in the center of 
an almost-regular tetrahedron, two vertices of 
which are occupied by Cl atoms and two by Natoms 
(Fig. 1). As seen from Table 1, the Co—Cl and 
Co—N bonds are predominantly covalent. It is ex- 


Fig. 1. Projection of the electron density of the CoCl, * 2H,NCgH,CHg, unit 
cell on the XOZ plane. 
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TABLE 1, Interatomic Distances in A 
F Between 
Asie neighboring | Val 1 
ei rin alen ngle 
the molecules . sotahins aii aa 
molecules 
Co — Cl 2.24 Cl —N3.40 | C]—Co—N 109° 
Co —- N 1.92 Cc) — CL 3.50 | Cl—-Co—Cl 114° 
N= €, 1,38 N—Co—N 111° 
Ci — C1 3.70 Co—N—C, 102° 
Cl—- N 3.40 
VN —N 3.15 


tremely probable that tetrahedral structure must 
also be ascribed to other addition compounds of co- 
balt halides of composition CoX,* A, and bluecolor, 
including the pyridinate. 

As regards the violet modifications, we should 
probably assign these a chain structure with coor- 
dination number 6, analogous to the structure of 
CoCl, - 2H,O [8]. 

It follows that the a- and 6-forms of the com- 
pounds CoX,° A, are not cis and trans isomers. 

This view regarding the structure of the two 
forms of cobalt dipyridinate was discussed earlier 
on the basis of magnetic measurements by Mellor 
[4]. In support of the fact that the blue color of 
compounds CoX,° A, is connected with the tetra- 
hedral configuration, we may cite the fact that, 
with 1,10-phenanthroline, which occupies two neigh- 
boring coordination sites, cobalt chloride gives a 
blue crystalline substance of composition CoCl,— 
C1.HgNp. 

The tetrahedral disposition of the Cl and N 
atoms around the Co in the ditoluidinate of cobalt 
chloride enables us to solve the problem regarding 
the nature of the bond in the so-called "normal" 
complexes [9, 10, 11], such as the ammoniates 
and aminates of the salts of divalent metals. It is 
usually considered that in these compounds the 
metal-addend bond is ionic—dipolar. In this case, 
the dipole should occupy a position of minimum po- 
tential energy relative to the metal ion, i.e., the 
metal ion should lie on the line of the vector of the 
dipole molecule. For the symmetrical molecules 
NH3 and C;H;N the direction of the vector of the di- 
pole molecule coincides with that of the metal —ni- 
trogen bond. In unsymmetrical molecules, such as 
aniline or p-toluidene, this is not so.. Our value 
for the Co—N—C, angle differs from 109° only with- 
in the limits of experimental error; hence, the four 
atoms coordinated around the nitrogen are disposed 
at the vertices of a tetrahedron, as in the salts of 
ammonia. The directional aspects of the nitrogen 


valences in the additional product CoCl, * 2C7HgN 
further confirm the covalent nature of the metal— 
nitrogen bond. 

Thus, from the crystallochemical point of view, 
there is no difference between normal complexes 
and interstitial complexes, which latter include the 
compounds of tervalent cobalt, platinum, chromi- 
um, etc., that are stable in aqueous solution. 

The authors consider it their pleasant duty to 
thank M. A. Porai-Koshits for consultation in the 
work. 
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An x-ray study has been carried out on crystals of indigo and thioindigo by 
two-dimensional electron density projections and by linear and planar sec- 


tions of the three-dimensional Fourier series. 


It has been shown that simple 


ao bonds and classical double bonds are not present in the indigo molecule. In 
the indigo molecule there is a single system of t-electron interactions. It 
has been established that there are intermolecular hydrogen bonds between 
the carbonyl and imino groups in indigo molecules. 


There is considerable interest in the current 
literature in the structure of one of the longest- 
known organic dyes, indigo, and of related sub- 
stances. The structural formula of indigo which 
correctly expresses the positions of the atoms in 
the molecule (1) was first established by Baeyer 
[1]. However, the indigo molecule has a number of 
properties which indicate that this formula does not 
express the actual character of the bonds and the 
distribution of the valences within the molecule. 
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These properties are: the absence of geometrical 
isomers, the intensity of color, the great stability 
of the molecule, etc. Since these properties are 
not explained by formula (1), many other formulas 
have been proposed for indigo in which the basic 
skeleton is maintained but the nature of the bonds is 
changed. Some authors connect the special proper- 
ties of indigo with an o-quinonoid zwitterionic struc- 
ture [2,3], while others connect them with the exist- 
ence of intramolecular hydrogen bonding between 
the carbonyl and imino groups [4-8]. A number of 
authors draw particular attention to intermolecular 
interactions [9-11]. 
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It seemed interesting to use the x-ray struc- 
tural method to elucidate some of the peculiarities 
of the indigo structure, since this method permits 
the determination of the spatial distribution of atoms 
and molecules in the crystal, and the interatomic 
distances and bond angles in the molecule. 

Recently, Eller [8, 12] published a number of 
short reports which cited x-ray structural data on 
indigo, thioindigo, and selenoindigo. We deciphered 
the crystal structures of indigo and thioindigo at 
the same time as Eller [13]. 

Crystals of indigo and thioindigo suitable for 
structural analysis were obtained by sublimation at 
atmospheric pressure. Crystallization from solu- 
tions did not give good results because of the ex- 
tremely limited solubility of indigo and thioindigo 
in normal solvents. 

The lattice constants,space groups, and num- 
bers of molecules in the unit cell are given in 
Table 1. They are in agreement with those pub- 
lished by Eller. It follows from these data that the 
indigo and thioindigo molecules in the sublimed 
crystals exist in the trans forms. But the sublima- 
tion was carried out under vigorous conditions, so 
that it would be expected that one isomer could be 
converted to the other. 

Our x-ray phase analysis showed that the dyes 
obtained by oxidation of the leucocompounds had the 
same structures as those obtained by sublimation. 
Consequently, it is concluded that the molecules of 
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TABLE 1. Lattice Parameters of Indigo and Thioindigo Crystals 


| Space | 
Substance | a | b | c | ay SOUP 
i = 
| 
Indigo Oe 2a eee eos) 4 P4Jc PD 
Thioindigo 7.82 | 3.92 | 20.34 2°40)" PAlc Z 
| 


indigo and thioindigo exist in the trans form in nor- 
mal crystals. 

The orientation of the thioindigo molecules in 
the unit cell was determined by geometric analysis 
and analysis of intensities. The experimental values 
of the structure amplitudes were obtained from in- 
tensities of reflections of the zero-layer lines on a 
photograph taken with rotation of the crystal around 
the b axis. The photographs were taken with cop- 
per radiation in an RGNS camera. In all, 141 re- 
flections of type h0l were obtained with maximum 
values of h= 10 andl = 26. From the coordinates 
for the chosen structural model, the signs of the 
structure amplitudes were chosen, and further 
structure decipherment was carried out from an 
electron-density projection on the (010) plane of the 
unit cell of the crystal. 

The method of successive approximations was 
used to refine the x and z coordinates. Four ap- 
proximations were carried out. The inclination of 
the molecules to the (010) plane was determined 
from the deformation of the benzene ring in the pro- 
jection, assuming that all the atoms were coplanar. 
In Table 2 the values of the x and z coordinates of 
the atoms of the thioindigo molecule are given in 
fractions of the lattice parameters. 

An electron-density projection along the b axis 
is shown in Fig. 1, and the values of the interatom- 
ic distances and bond angles are shown in Fig. 3. 
In Fig. 3 the experimental and calculated structure 
amplitudes are compared for certain regions of 
sind. 

It was established by geometrical analysis that 
the plane of the indigo molecule in the crystal is at 
a considerable inclination to all three coordination 
planes of the unit cell. Hence, three-dimensional 
electron-density series were used for solution of 
the crystal structure. 

Photographs were taken in copper radiation 
with an RGNS camera with rotation of the crystal 
around the b axis, and the layer lines hdl, hil, 
h2l, h3l, and h4l were developed, 707 reflections 
in all. The signs of the h0/ structure amplitudes 
were first determined from the structural model 
found from intensity analysis, and were then re- 
fined from the atomic coordinates in the electron- 
density projection on the (010) plane. The method 


TABLE 2. Atomic Coordinates (x, z) in the 
Thioindigo Molecule 


No.of No. of 

atom x/a z/c atom x/a z/c 

Ci 0.454 | 0.023 Ce, 0,068) 0.4177 
0,513 |0,C94 C; 0.071] 0.113 

HOON ODI aC  |\Os2u7l) OxXos 

Gy 0,370 | 0.188 s 0.230} 0,009 

C; 0,226 | 0.214 O 0.656} 0.108 


Fig. 1. Electron density projection of the thioindigo crystal along 
the b axis. 


Fig. 2, The thioindigo molecule: values of the interatomic dis- 
tances and the bond angles, 


of successive approximations was used to deter- 
mine exact values of the x and z coordinates. It 
was established from the size of the benzene ring 
that the plane of the molecule coincided, or almost 
coincided, with the (210) plane. The signs of all 
the structure amplitudes were determined from the 
x and z coordinates [from electron-density projec- 
tion on the (010) plane] and the coordinates y = 1 

— 2x, which correspond to the position of the mole- 
cule in the (210) plane. A section of the three- 
dimensional electron-density series on the plane 
(210) (Fig. 4), and ten linear sections parallel to 
the y axis, passing through centers of atoms, were 
then constructed. The atomic coordinates, deter- 
mined from the projection and sections of the elec- 
tron density, are given in Table 3. 


Figure 5 depicts the molecule of indigo and in- 
cludes bond length and bond angle data. In Fig. 6 
is given a comparison of the experimental and cal- 
culated structure amplitudes of the h0l type for 
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Fig. 3. Comparison of the experimental and calculated structure amplitudes for some regions of sin 6 in thioindigo. 
Solid line — calculated; dashed line — experimental. 


Fig. 4. Section on the plane (210) of the three-dimensional 
electron density series for the indigo crystal. 


TABLE 3, Atomic Coordinates (x, y, z) in the Indigo Molecule 


atom x/a y/b z/¢c x/a y/b zl¢c 


Cy 0.049 | 0.907 | 0.020 Cy, 0.373 | 0.2641 | 0.087 
C, 0.103 | 0.793 | 0.142 C, 0,293 |0.432 | 0.004 
C; 0.197 | 0.593 | 0.138 Cy 0.208 |0.602 | 0.029 
Cy OL2HT* WOL420 5) 08222 N 0.124 |0.793 | 0,038 
C; 0,364 | 0.256 | 0.196 O 0.059 | 0.847 | 0,219 


certain regions of sin@. It has been established 
that all the atoms of each half of the molecule (ex- 
cept nitrogen and Cy and Cz and the center of sym- 
metry) lie in a single plane (with a precision of a 
few thousandths of an Angstrom unit), and this 
plane is defined by the formulas 
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Fig. 5. Indigo molecule, showing interatomic distances and 
bond angles. 


The divergence of the atoms N, Cy, Cz, and the 
center of symmetry from this plane are, respec- 
tively, 10.07 A, 40.05 A, 40.03 A, and -0.03 A. 

An analysis of the interatomic distances in the 
indigo molecule shows that there are no simple o 
bonds (except the bond to hydrogen) nor classical 
double bonds. The bond lengths are intermediate 
between single and double in both the benzene ring 
and in the heterocyclic ring. The fact that the 
C()—N and N—C,) distances are shorter than those 
characteristic of an aliphatic single bond (1.47 A) 
may be explained by formation of 7 -bonds between 
the unshared pair on nitrogen and its neighboring 
atoms. The very weak basicity of the imino group 
in indigo confirms this view. There is but a single 
system of 7-electronic interactions in indigo. The 
results obtained show that there is no reality in the 
proposed o-quinonoid zwitterionic structure. 
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Fig. 6. Comparison of the experimental and calculated structural amplitudes of the h0J type 
for some regions of sin@ in indigo. Solid line — calculated; dashed line — experimental. 


Less conjugation in the heterocyclic ring was 
observed in the thioindigo molecule. The inter- 
atomic distances between carbon and sulfur (1.81 A) 
show that the C;—S and S—Cg, bonds cannot be con- 
sidered to be intermediate between single and 
double. The sulfur atom either plays no part or 
plays a negligible part in the 7-bond interaction. 

The N...O distance within the indigo molecule 
is 2.94 A, i.e., approximately equal to the sum of 
the intermolecular radii of nitrogen and oxygen. The 
position of the hydrogen atom in the iminio group, 
like those of the other hydrogen atoms, was not de- 
termined. But, to judge from the distance between 
the oxygen and nitrogen atoms, and from the rela- 
tive position of imino and carbonyl groups, there 
are no internal hydrogen bonds in the indigo mole- 
cule, or in any case they are very weak. This con- 
clusion is confirmed by a study of the vibrational 
absorption spectra of indigo, its tetrahalogenoderiva- 
tives, and thioindigo [14]. 

In the indigo crystal, each nitrogen atom is at 
a distance of 2.88 A from an oxygen atom in the 
neighboring molecule, related by a glide plane. This 
indicates the existence of intermolecular hydrogen 
bonds between the imino and carbonyl groups of in- 
digo molecules. Each molecule is linked by hydro- 
gen bonds to four neighboring molecules. Therefore, 
bonds appear between all molecules whose centers 
of symmetry lie on the same plane perpendicular to 
the a axis of the crystal. It is natural that all these 
hydrogen bonds in the indigo crystal should increase 
its thermal stability. It may also explain the higher 
melting point of indigo (8390-392°) relative to thio- 
indigo (288°). 

Some authors have observed in x-ray structur- 
al studies of molecular crystals that the heights of 


the electron-density maximum for the same atom 
depends on its position in the molecule [15, 16]. A 
similar phenomenon was observed in this study of 
indigo and thioindigo. The electron density at 
maxima decreased steadily as the atom was moved 
away from the center to the periphery of the mole- 
cule. The maximum difference in heights of 
maxima was 17% in indigo and 28% in thioindigo. The 
difference in heights of maxima of atoms at equal 
distances from the center of the molecule did not 
exceed 3% of the average for thioindigo and 5% for 
indigo. It may be suggested that the decrease in 
electron density from the center to the periphery of 
the molecules observed in the x-ray study of indigo 
and thioindigo crystals is connected with thermal 
vibrations of the whole molecule about the center of 
gravity. 

The further the atoms from the center of the 
molecule, the greater their limits of vibration and 
the greater the blurring of the electron density. In 
the indigo crystal in which the molecules are linked 
by hydrogen bonds, thermal vibration may not be 
as strong as in thioindigo and, therefore, the de- 
crease in the height of electron-density maxima is 
not so great for indigo. 

On the basis of the x-ray structural study, it 
can be concluded that if cis isomers of indigo and 
thioindigo exist, the distance between oxygenatoms 
in their molecules would be about 2.0-2.1 A. Such 
short distances between nonbonded atoms with the 
same excess charge cannot exist in molecules. For 
there to be larger distances between the oxygen 
atoms in the cis isomers, there would have to be 
strong distortion of the bonding angles or destruc- 
tion of the coplanarity of the molecule, and this 
is energetically unsuitable. This may be the cause 
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of the lack of geometrical isomers in indigo and 


thioindigo. 


Recently, Brode, Pearson, and Wiman 


[7] reported the preparation of the cis isomer of 
thioindigo and some of its derivatives, but the re- 


sults of their study are not completely unambiguous. 


wn 
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The lattice parameters of methylene blue [CygHigNgS)C1- - nH,O have been 
refined to a = 9.68,, b = 31.86), c= 7.07,A, and B = 97°11'. An electron- 
density projection on the plane (001) has been constructed from F-series 
data and by isomorphous replacement of the chlorine atom by bromine. From 


this the ionic structure of methylene blue 


has been confirmed and it has been 


established that the atom of the cation closest to halogen is a sulfur atom. 


The structures of dyes have seldom been 
studied. A systematic x-ray structural study would 
be of considerable interest for the elucidation of the 
relation between structure and properties of dyes. 

In the present work, we have studied by x-ray 
methods the structure of the thiazine dye, methyl- 
ene blue, which is much used in practice. 

In x-ray and physicochemical papers, two 
forms of the formula have been suggested: 


Cl 
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Cl 


Absorption spectroscopic studies favor the ion- 
ic structure [1]. In this paper it was suggested that 
on hydrolysis the OH” group became attached to the 
central nitrogen atom of the molecule and the H* to 
a nitrogen atom of an amino group. The charge lo- 
calization in the complex ion was not established. 

In many chemical papers it has been suggested that 
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the positive charge is averaged over the three nitro- 
gen atoms. 

Taylor [2] has carried out x-ray and optical 
studies of the methylene blue chloride, bromide, 
and iodide, and for the selenium-substituted brom- 
ide. He established that all the crystals except the 
iodide were twins. The twinning plane was parallel 
to the direction of elongationofthecrystal. Taylor's 
data are cited in Table 1. They establish that the 
chloride, bromide, and selenium-substituted brom- 
ide are isomorphous. 

Structure determination was carried out only 
for the iodide I, and this was by trial and error. The 
projection of the structure of iodide I on the (001) 
plane is shown in Fig. 1. The figures represent 
atomic coordinates along the Z axis. 

All the atoms of the complex ion [C,gH;sSN3]* 
lie in a single plane and the ions are associated in 
chains. The chains are distributed in layers paral- 
lel to the plane (012) and at a distance 345 A apart. 
Taylor suggested that all the compounds in Table 1 
were isostructural, i.e., they had a layer structure 
and consisted of the ions [C;gH,,N3S]* with the methyl 
groups as points of contact. The nearest neighbors 
of the halogen atoms are the methyl groups and the 
water molecules. 

To elucidate the structure of the methylene 
blue dye [CjgHjgN3S]‘Cl° - nH,O, we determined the 
structure by the electron-density series method. 
The dye crystals were obtained from the colorim- 
etry laboratory of the Scientific Research Institute 
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TABLE 1. 


Substance 


Chloride 
Bromide 
Iodide I 

Iodide II 


Seleno-bromide 


Figs 


for Organic Semifinished Materials and Dyes. Sev- 
eral recrystallizations were carried out by V. P. 
Glushkova to obtain the best x-ray quality. X-ray 
photographs were taken of the industrial and re- 
crystallized products. 

The crystals were deep blue needles with a 
bronze metallic luster, 1-2 mm long, and 0.2-0.3 
mm in thickness. The crystals were very brittle, 
and most appeared as fragments of needles. 

Refinement of the lattice constants was carried 
out from an x-ray oscillation photograph taken on 
an RKU camera constructed by M. M. Umanskii. 
The following values were obtained: a= 9.68,, b 
= 31.865, c= 70 GN B = 97°11". The photograph 
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Monoclinic Ch P2,/c 
Monoclinic | C?,—P2,/c 1.41 
Triclinic | Ci—pi 1.68 
Triclinic/ | Not deter. 1.64 
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Monoclinic | C3, —P2/c 157 


was indexed in terms of the space group Cey= P2,/c. 
Ninety-eight reflections were revealed by develop- 
ment of the zero-layer lines relative to the c axis 

in Cu radiation. Three x-ray photographs with ex- 
posures of 18, 36, and 70 h were taken for more 
precise estimation of the intensities. The intensi- 
ties were estimated quantitatively by comparison 
with a standard blackening scale. 

Projections of the interatomic vector functions 
on the planes (001) (Fig. 2) and (100) (Fig. 3) were 
constructed for methylene blue chloride. 

Since the number of atoms in the unit cell is 84 
(excluding the hydrogen atoms), there is consider- 
able superposition of maxima in the projection on 
(001), as was to be expected. The peaks belonging 
to the heavy Cl and S atoms were revealed, but it 
was not possible to distinguish between them be- 
cause of their very similar diffracting power. 

The isomorphous replacement of chlorine by 
bromine was used to find the coordinates of the Cl 
atoms from the F* series. The approximate co- 
ordinates of the Cl (Br) and S atoms were deter- 
mined by a comparison of the peak heights of the in- 
teratomic vectors for the chloride and bromide: 


Sep=0208.8 Yo, — 0eal, 


xg = 0.21, ys = 0.19. 


The projection of the interatomic vector function on 
the plane (100) showed that the structure was layered. 
These coordinates of the halogen and sulfur 

atoms permit the determination of the signs of 78 
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structure amplitudes from a comparison with the 
moduli of the structure amplitudes for the chloride 
and bromide with uniform (hk0). The electron-den- 
sity projections were found by using these signs. 
The electron-density distribution diagram for the 

HC ratiaNas] Cl. - nH,O crystal in the projection of 
the (001) plane (Fig. 4a) confirmed the Cl and S co- 
ordinates found previously from the F* series: 


Cl (0.94; 0.21), 
S (0.23; 0.19), 


and allowed some conclusions to be drawn about the 
molecular structure of methylene blue: the complex 
ion is extended along the y axis. The distribution of 
the electron-density maxima corresponding to the 


projection of this structure model is shown in Fig.4b. 


The distribution of regions of negative electron 
density correspond to the projections of the central 
"holes" in the benzene and heterocyclic rings. It 
must be taken into account that severe superposi- 
tion of maxima occurs, especially in the part of the 
projection, as a result of projections of molecules 
from different layers. The false maximum at (x 

= 0.21, y = 0.25) probably occurs from a break in 
the series and the probable errors in determining 
the signs of some of the structural amplitudes. 


A geometric analysis was carried out to con- 
firm the structural model. A model of the complex 
ion [CygH;gN3S]* (Fig. 5) was constructed from 
literature values for interatomic and intermolecu- 
lar radii, taking account of errors arising from the 
directionality of covalent bonds. The observed 
electron density distributions (Figs. 4a and 4b) 
agreed with the data from the geometrical analysis 
(Fig. 6). 

Thus, by the successive use of the following 
methods: construction of the interatomic vector 
function, isomorphous replacement, construction 
of electron-density projections, and geometrical 
analysis — the principles of the structure of methyl- 
ene blue were established for the (001) projection. 

Analysis of the S—Cl distance even in the ap- 
proximate model of the structure in the (001) pro- 
jection, makes possible the conclusion that the S— 
—Cl bond is not covalent, since the closest S—Cl 
distance, 2.8 A in projection, is considerably 
longer than that of a covalent bond (2. 03 A). Con- 
sequently, the chemical formula of methylene blue 
should be given in the form [CygH;gN3S]'Cl” - nH,O. 

It is also seen from Fig. 5 that the sulfuratom 
of the cation is closest to the chloride ion, and not 
the methyl groups as Taylor suggested. This is 
possibly explained by placement of the positive 
charge on sulfur. An earlier study of complex 
metal thiocyanates [3] established that the electron 
density was redistributed between the sulfur and ni- 
trogen atoms in the covalent SCN group and in the 
SCN” ion as a result of which the nitrogen atom had 
an excess of electron density in comparison with 
sulfur — (N~1) and (S~2), respectively. A similar 
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Fig. 4b. 
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Fig. 6. 


redistribution of electron density has also been es- 
tablished in the molecules N,S, and As,S, [4]. Thus 
our work has not confirmed Taylor's suggestion 
that methylene blue chloride and iodide are iso- 
structural. 

It is possible that the trial-and-error method 
used by Taylor can lead to an incorrect structural 
model when applied to a complex model with 75 
atomic parameters. 

A further study of the structure of [CigH,sN3S]° 
* Cl” + nH,O by the construction of a three-dimen- 
sional electron-density section will allow for re- 
finement of the interatomic distances and bonding - 
angles, and will reveal more precise details of the 
structure. 
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During electron-diffraction study of a-Fe films subjected to cementation at 
400 to 480°C in a current of CO or at approximately 650 to 700°C in a current 
of acetylene, a new cubic carbide of iron with lattice spacing a = 3.878 

+ 0.002 A was observed. Following structural analysis by means of 2 and® 
series, and also minimization of the confidence factor for the new carbide 

(of probable composition Fe,C), a structure based on the slightly distorted 
lattice of a-Fe was established. There was a single "molecule" of Fe,C in 
the unit cell. The Fedorov group was 1 1C in (a):000; 4Fe in (e): rez 
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Despite the large number of investigations de- 
voted to the structure of the crystal lattices of va- 
rious interstitial phases, the results obtained are 
by no means complete. This is especially true of 
the carbides of iron. It is evidently generally ac- 
cepted that at temperatures above 721°C austenite 
is a stable carbide. Below this temperature we 
know martensite and cementite as metastable 
phases. The only stable phase in this temperature 
range is the cubic @ solid solution with a negligible 
carbon content. . 

However, other carbides of iron undoubtedly 
exist. Investigations concerning these seem to pro- 
ceed in two different directions: the study of pro- 
cesses and structure during the tempering of 
quenched steels, and the structure of iron catalysts 
used in the synthesis of hydrocarbons. We may 
take the existence of iron carbides richer in carbon 
than cementite (Fe3;C) as established. The well- 
known carbide of Hagg, and also the hexagonal carb- 
ide are assigned the composition Fe,C [1,2,3]. 
There is also evidence of the existence of a carbide 
with formula Fey)C9. For the hexagonal carbide we 
have lattice spacings a = 2.74 and c = 4. 32 A. For 
the analogous iron nitride Fe.N, the spacings are 
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The carbon had tetragonal coordination with Fe—C distance 


a=2.77 and c = 4.42 (a which disagrees with the 
usually given atomic radii: carbon 0.77 and nitro- 
gen 0.71 A. The incomplete nature of structural 
data on the carbides of iron is demonstrated by the 
absence of any indication of lattice spacings for 
Hagg's carbide, although the list of interplane spac- 
ings for reflections ascribed to this phase was first 
published more than twenty years ago and repro- 
duced in other work. All data on the carbides of 
iron are usually based on Debye photographs from 
samples constituting multiphase systems. There 
are no reliable determinations of the chemical com- 
position of the separate phases. All this is insharp 
contrast to the remarkable achievements of struc- 
tural x-ray diffraction in other fields. 

In addition to purely structural results, there 
are a considerable number of chemical and physico- 
chemical investigations into the interaction of iron 
and its oxides with CO and H, (synthesis of hydro- 
carbons) which have been reported.! Various 
authors have suggested the possible formation of 
the following compounds: Fe,C, Fe.C, FeC, FeC,, 


1) iterature indicated in review [4]. 
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Fe3C, and Fe,C3. These results, however, require 
structural verification. 

We note that the structures of long-established 
carbide phases are also insufficiently well known. 
In particular, the position of the carbon in the aus- 
tenite, cementite, and martensite lattices has 
never been reliably established by direct experi- 
mental determination. The data for austenite and 
martensite are not well founded, since the propor- 
tion of x-ray scattering introduced by the carbon 
atoms is of the same order as the error in meas- 
uring (and calculating) the intensities of the diffrac- 
tion photographs. For cementite, the position of 
the carbon atoms in the lattice has been established 
only from geometrical considerations. 

The inadequacy of structural data for the carb- 
ides of iron, together with the extreme importance 
of an all-around knowledge of these phases and the 
difficulties found in x-ray examination, led us to de- 
velop a new method of studying them by electron 
diffraction. After we had obtained experimental 
data on the new carbide, some other papers were 
published on the electron-diffraction study of iron— 
carbide structures, using a completely analogous 
method. In these papers, however, only the forma- 
tion of cementite and the hexagonal carbide was 
mentioned, and no new structural data were pre- 
sented [5]. 

In order to obtain iron carbides, we subjected 
iron films with some degree of orientation to treat- 
ment in the gas phase. Treatment in a current of 
CO or CO +H, was used, as well as cementation in 
a current of acetylene. 

In developing the method of study, we tried to 
create conditions favorable toward as detailed as 
possible a structural examination of the carbide 
phases formed. Results obtained for the structure 
of nitrides illustrate the possibilities created by 
electron-diffraction analysis of films prepared by 
treating iron in a gaseous medium [6]. The origi- 
nal layers of single-crystal or polycrystalline iron 
were obtained by condensation from the vapor onto 
the surface of heated rocksalt crystals. After con- 
trol photographs were taken of films stripped from 
the salt crystals, samples of crystals bearing iron 
films were placed in a furnace for cementation. The 
carbon monoxide was produced by decomposing 
formic acid, and was purified and dried. Before 
passing the carbon monoxide, the system was 
flushed with pure, dry hydrogen. Interaction with 
the CO was effected at temperatures from 400 to 
~520° for various periods from 5 min to 5h. Treat- 
ment with dry acetylene was effected at 400 to 500° 
and 650 to 700° for 15 min. Then the carbide films 


were stripped from the salt crystals in the usual 
way and placed in the electron-diffraction camera 
for study. Structural study showed that, apartfrom 
oxides of iron formed by oxygen accidentally drawn 
in from the air, this treatment gave three carbide 
phases, including cementite and the hexagonal carb- 
ide. As regards the third phase, we were able to ob- 
tain electron-diffraction photographs containing only 
the lines of this phase on treatment in a current of 
CO at 400 and 480°, and also on 15-min treatment in 
a current of acetylene at 700° (Figs. 1-5). In some 
samples obtained by treatment with CO, we noted 
partial orientation. Apart from this, with this 
treatment we obtained some films containing com- 
pletely nonoriented crystals of the new phase. We 
must assume that the degree of orientation of the 


Fig. 1. Electron-diffraction photograph of a two- 
phase film containing Fe,C (mosaic single crystal) 
and oriented cementite, 


Fig. 2. Electron-diffraction photograph of a part- 
ly oriented film of the carbide Fe,C. 
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Fig. 3. Electron-diffraction photograph of a part- 
ly oriented film of the carbide Fe,C, 


crystalline film of FeyC containing NaCl lines. 


samples does not depend on the composition of the 
cementing gas, but is determined by the structure 
of the original iron film and various other factors. 

Electron-diffraction photographs of the new 
phase clearly demonstrate its cubic symmetry. The 
network of reflections observable in the diffraction 
picture, where the main pattern appertains to Fe3Q,, 
is immediately indexed by means of indices hk0. 
The electron-diffraction photographs of Figs. 2 and 
3, with partial orientation of the scattering crystal- 
lites, show texture typical of cubic phases: the di- 
rections [111], [110], and — to a smaller extent — 
[100] lie perpendicular to the plane of the film. 

In determining the lattice spacing of the new 
phase, we used electron-diffraction photographs 


Fig. 5. Electron-diffraction photograph of a poly- 
crystalline film of Fe,C; reflections from NaCl 
are noticeable in the form of spotty lines. 


containing NaCl lines as well as those of the new 
phase itself. The lattice spacing proved to be a 

= 3.878 + 0.002 A. An extremely important result 
is the constancy of this quantity on determining the 
lattice spacing of the new phase in a series of ex- 
periments carried out under different conditions 
over a number of years. 

The experimental intensities of the electron- 
diffraction photographs were estimated on the well- 
known ten-point scale. 

Electron-diffraction photographs obtained in 
the series of experiments using carbon monoxide 
(Figs. 2 and 3) correspond to partly oriented 
samples. The diffraction pattern is formed by the 
superposition of more intense arcs on the continu- 
ous Debye ring. Fairly perfect orientation of the 
crystallites produced short arcs, and on all the 
rings (though less so for planes with a large recur- 
rence factor), the sections between the intense arcs 
could be used for estimating intensities. Error in 
intensity measurement is nevertheless increased by 
these features of the diffraction pictures. 

The electron-diffraction photographs of the 
nonoriented polycrystalline specimens (Figs. 4 and 
5) obtained in experiments with acetylene also 
proved unsuitable for exact intensity measurements. 
In some photographs the superposition of NaCl lines 
caused confusion; in others, there was serious fog- 
ging. It should be noted that, in order to determine 
the position of the light carbon atoms in the pres- 
ence of iron (especially for composition Fe,C), 
microphotometer intensity measurements are re- 
quired if precise results and good agreement be- 
tween theory and experiment are to be achieved. 


DETERMINATION OF IRON CARBIDE Fe,C STRUCTURE 51 


Our Fourier synthesis based on experimental 2 and 
@ values has the advantage of being less sensitive 
to the accuracy of intensity measurements. Table1l 
shows the corresponding averaged data obtained by 
using a large number of electron-diffraction photo- 
graphs. 

Considering the data of Table 1, we may con- 
clude that our carbide has a lattice based on y iron. 
In favor of this supposition is the value of the lat- 
tice spacing a = 3.878 A, which is in reasonable 
agreement with those of austenite (3.62 A) and the 
y phase of Fe,N (3. 79A). We also note that reflec- 
tions with nonmixed indices in general have greater 
intensities than the rest. We calculated one-dimen- 
sional sections of the 6* series along the [110] and 
[111] directions. The first section gave a quite de- 
finite result: a strong maximum in the center of the 
face. As regards the second, in the [111] direction 
there was a slight maximum near the center of the 
octant A A AY 

The sensitivity of the ®* function was insuffi- 
cient to fix the coordinates of the end of the Fe—C 
vector reliably. 

It should be noted that the data presented do 
not contradict a lattice of three iron atoms in the 
unit cell, disposed in the centers of the faces,which 
corresponds, for example, to the group O(c). Con- 
sidering the possible composition and structure of 
the carbide on the basis of such a lattice with acar- 
bon atom at the vertex of the cube, we arrive at 
the extremely unlikely Fe—C distance of 2.74 A. 
Hence, without further investigation, considering 
this lattice as extremely loose, we reject this vari- 
ant. 

Just as improbable is the structure on the a- 
Fe basis. 

Turning to the possible distribution of the car- 
bon atoms in the fcc lattice of iron, and to the 
composition of the present carbide, we must also 
reject the simplest solution: a carbon atom in the 
octahedral position (center of the cell). Sucha 
structure should give extremely weak reflections 
with mixed indices, rapidly diminishing with in- 
creasing hkl; here, the scattering power 4ffe+/o 
compares withfcq. The fairly large intensity of re- 
flections with mixed indices compels us to re- 
linquish this model. We note furthermore that the 
considerable intensity of the 110 and 210 reflec- 
tions may be determined by the considerable scat- 
tering power of carbon at low angles. 

It is clear that, on placing one carbon atom in 
a tetrahedral vacancy, we should expect a similar 
scattering picture. 


TABLE 1. Experimental Electron-Diffraction Data for the Carbide 


Fe,C 
kkl | dpxr ind aes akl | dpyyinA Peep Si 
uk nk in A =Viexp [Pt hk in =Viexp [PF 
100 | 3.875 4.8 422 0.79 17.1 
110 | 2.74 ils si 430 0.775 2,4 
JA 2A 2od 34.5 500 0.775 3.7 
200 | 1.935 24.6 431 0.76 2.0 
Z1O A732, 4,7 333 0,746 12.8 
211 | 1.58 Use 544 0.746 sa 
22 ON nef 25.6 520) 0.72 2.6 
SOOM R29 5.8 021 0.708 1.8 
2PM Ny B28) 2.3 440 0.685 9.2 
310 | 1,224 Dl 522 0.675 1.9 
341 1.169 20.3 441 0.675 1.9 
222) es 10.4 530 0.665 1.9 
230 | 1.074 4.4 433 0.665 io 
231 1.035 Sarl 534 0.655 6.8 
400 | 0.968 9.4 G00 0.646 8.9 
410 | 0,94 0 442 0.646 Oa 
322 | 0.94 4.7 611 0.627 4.6 
441 | 0.943 3.9 620 0.612 7.6 
330 | 0,913 Aer? 533 0.85 6.0 
331 | 0.888 17.3 622 0.583 5.0 
420 | 0.866 13.9 
424 | 0.845 3.0 
332 | 0.8255 4.8 | 


Is it possible to have a number of carbonatoms 
greater than unity in the unit cell? The symmetry 
of the cubic groups ane Ti o!," and OF admits 
placing three carbon atoms in octahedral vacancies 
(centers of edges). In this case, however, the re- 
flections with unmixed indices are strengthened, 
while the others will again be determined by the 
scattering of only one carbon atom. Clearly, by 
placing four carbon atoms in octahedral vacancies, 
we obtain an NaCl-type lattice with corresponding 
extinctions. 

The only solution which we can suggest is based 
on the fact that the iron atoms also take part in the 
formation of reflections with mixed indices, which 
is possible with a slight rearrangement of the ideal 
fec lattice. This rearrangement may be called the 
"implantation" of a carbon atom into a tetrahedral 
vacancy, Since, for a period a ~ 3.88 A the dis- 
tance from the center of an iron atom to the center 
of the nearest tetrahedron is 1.675 A. Thus we 
find a carbide of composition Fe,C with carbon ina 
tetrahedron and four iron atoms displaced from the 
ideal positions along the space diagonals. This 
structure is described by group ce The carbon is 
in (a) :000 and the iron in €) :\zzrq). In order to 
check this model and determine the value of x, we 
examined the ® series from the experimental struc- 
ture amplitudes. The model proposed has no cen- 
ter of symmetry; hence, calculation of the three- 
dimensional series requires knowing the phases of 
the structure amplitudes. So as not to introduce a 
substantial arbitrary element into the calculation in 
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view of the sensitivity of the phase values to the pa- For the parameter value in question, the mini- 
rameter taken for iron, we decided to limit the in- mum Fe-—C distance in the lattice is 1. 78 A, for 
vestigation to the centrosymmetrical projection of the "average" distance in austenite 1.81 A, and for 
the structure on the cube face. cementite, 1.86 A. The small Fe—C distance and 
In calculating the theoretical structure ampli- the tetrahedral coordination of the carbon atom in 


tudes, we met some indeterminacy of the data for 
the scattering power of the carbon atoms. In fact, 
the literature contains no less than four tables for 
i ne calculated by different authors. Table 2 gives 
the corresponding values of the atomic amplitudes 
for electrons. The most serious discrepancies oc- 
cur for small angles of reflection. We note that 
this fact introduces difficulty into accounting for 
the cationor anion state of the carbon atom in cal- 
culating the theoretical @. In calculating ®, we 
used the values II in Table 2. As calculation of the 
structure amplitudes showed, the signs of ®y,7, 
except for two weak reflections, did not vary over 
a wide range of values of the parameter x. Figure 6, 
representing the projection of the lattice potential 
of our carbide on the cube face, gives an unequivo- 
cal solution for its structure; the four maxima near Fig. 6. Projection of the potential of the FegC 
the centers of the quadrants correspond to iron lattice on the cube face, calculated from ex- 
atoms, and the maxima in the vertices of the cell perimental S.v ales, 

to carbon (Fig. 7). The parameter x for the iron 

atoms, as determined from this projection, is 0.26. 

By using the method of minimizing the confidence 


factor R= >)|| ME) —| >| / >| PE | we established 


that the most probable value was 0.265. For this, 
the R is quite large (0.28), which of course is 
chiefly due to insufficient accuracy in estimating 
the experimental reflection intensities. The over- 


all picture of the relation between theoretical and Or Fe 
experimental structure amplitudes is shown in 
Pigs .8. Fig. 7. Structure of the carbide Fe,C. 


TABLE 2, Atomic Scattering Amplitudes of Electrons by Carbon f =(Z—f pi? from 
Various Data 


| 
sin %/2-10® | 0.0 | 0.05 | 0.7 | 0.15 | 0.20 | 0.25] 0.30 | 0.40 | 0.50 | 0.60 | 0.70 | 0.80 | a.90 | 41.00 
10 Dd | 3,33) 2,5 | 2) | 1.67) 1.25) 4.00 |0.833/0.714/0.62510.555/0.500 
13 52 35 18.75 10.55] 6.40] 4.30]3.06 |2.34 |1.84 [4.49 |1.25 
Wigcaee 23.8 | 21.9]19.9/16.8 |14.9 |12.2} 9.7 |6.54)4.5 3.10 |2.35 |4.84 [1.49 |1.25 
Te s* 80 |60 {41 25.9 [17.2|12,2 | 7.22] 4,83 
Lyi *%2* 23.6 |21.5}18.05]14.9 }12.0] 9.65) 6.26] 4.29 13.10 


*X-ray amplitudes by Hartree, Internationale Tabellen z, Bestimmung v. Kristall- 
strukturen, No, II, 571 (1935), 

**Taken from tables of V. K. Vainshtein, “Atomic scattering factors of electrons," 
ZhETF, 25, 2, 157-168 (1953), 

***X-ray amplitudes by MacWeeny. “X-ray scattering of aggregates of bound 
atoms,” Acta Cryst., 7, 2, 180-186 (1954), 

****X-ray amplitudes by I, Berghuis, I. Iberta, et al,, "New calculations of atomic 
scattering factors," Acta Cryst., 8, 8, 478-483 (1955), 


DETERMINATION OF IRON CARBIDE Fe,C STRUCTURE 53 


~ NBER IBDOSU SS EGS SSSLRN 


Fig. 8, Relation between theoretical and experimental structure amplitudes, 


the lattice possibly indicate that this carbide is a 


chemical compound and distinguish it from solid so- 


lutions such as austenite. This hypothesis is also 
supported by the constancy of the lattice spacing of 
the carbide Fe,C, which indicates a constant com- 
position. 

In conclusion, we express thanks to former 
graduate students of Gor'kii State University, E. 
Krasotskaya, M. Grinberg, and G. Pinsker, and 
also to ourco-workerat the Gor'kii Physicotech- 
nical-Research Institute at the University, N.V. 
Gudkova, who took part in obtaining the experi- 
mental material. We also express thanks to E. L. 
Lapidus, who carried out a large part of the com- 
puting work. 
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The crystal structure of the boron trifluoride pyridine complex has been deter- 


mined. 


The parameters of the monoclinic lattice are:. a= 17.71, b= 5.89, c 


= 14.34 A; B =118°42'. The number of molecules per unit cell was z= 8. The 
space group was C2), — P2,/c. The method of statistical equalities was used to 
determine the signs of the structure amplitudes and also the method of three- 
The interatomic distances found for 
the compound F3;B > NC;H; were: B—F, 1.41; B—N, 1.53; F—F, 2.30; N-—C, 
1.36; C—C, 1.39 iG and the bond angles were {FBF = 109.5°, {NBF = 108.5°, 


dimensional electron-density sections. 


X C,NC, = 117. 5s X NC,Cs3 = WAS J C,C3C5 = ILS). J C3C5Cy = son Oe. 


tallochemical theory of donor—acceptor bonds is discussed. 


There are many important theoretical and prac- 
tical problems in the field of physical chemistry 
for crystallochemical study. One of the central 
problems of chemical bonding theory is the elucida- 
tion of the nature of donor—acceptor bonds in com- 
plex compounds. The theoretical side of this prob- 
lem provides for an explanation of the part played 
by the unshared pairs of atoms in molecule forma- 
tion and the discovery of the nature of the chemical 
bond formed between the donor and acceptor mole- 
cules. The practical side is connected with broad 
problems of synthesizing materials with a variety 
of technological properties. 

The object of x-ray structural study in the 
present work was the crystalline complex formed 
from boron trifluoride and pyridine, F3B ~-NC;Hs. 
Boron trifluoride and its derivatives are very im- 
portant catalysts which have been studied vigorous- 
ly of late. It is sufficient to mention that there are 
more than a thousand papers and monographs on 
the study of these compounds [1,2]. Of the x-ray 
structural studies, those of Hoard and his co- 
workers must be noted first. They studied F3B 
~——NH3, F3B ~~NH,CH3, F3B ~-N(CHs3)3, and F3B 
—-NCCH;. They cited the results of their studies 
in a review [3] as interatomic distances and bond 
angles (Table 1). 


04 


The crys- 
TABLE 1 

Interatomic 

Coie: F,B<“NH, | F,3B-NU.CH;| F2B<N (CH3).| F2B-NCCH, 
angles 
N—B 1.60 Lm 1.585 1.635 
B—F 1.38 1.39 4.39 deo 
F—F Dene PAPA | DOK Dees 
N—F 2.40 2.40 2.46 2.34 
N —CHy, — 1.50 1.50 os 

A. FBF 411° 410552 107° eA 
ne FBN 107° 108.5° DS 103° 
\ BNC — 114° 105° —- 
1 CNC — = 114° = 
\ CNB — 144 105° -— 


It is known that the nitrogen atom is found in 
molecules in three basic valence states: (1) three 


/ 
o bonds Ne , &.g., ammonia, NH3; (2) two o and 
7 
one 7 bond Ne , €.g., pyridine, NC;H,; and, (3) 


one 0 and 2 7 bonds, N=, e.g., methyl cyanide, 
CH3;CN. The crystalline compounds studied by 
Hoard contain nitrogen atoms in the first and third 
valence states. The current data on the dissocia- 
tion energy of the compounds F3;B—NC;H,, 50.6 
keal [4], F3B -N(CHg)3, 28 kcal; and F;B ~NHg, 
17 + 10°keal, and the energy of formation of the 
crystalline compounds F3B -NCCHg, 26.5 kcal, 
F3B ~-NH,CH3, 42.7 kcal, and F;B ~-NHs, 42 kcal, 
indicate that the greatest effect of the donor —ac- 
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TABLE 2 
Formula | System ainA |binA | cin A g° Cin Fn | z | Space group 
FB — NC;H;| monoclinic 17,71| 5.89 14.34 118°42’| 1.49 | 1,48)| § “C3, = Pave 
ceptor bond is observed for the second valence TABLE 3 
state of nitrogen. It was therefore of interest for | 
crystallochemical correlations to study the struc- aes pres Pa ies ue 
ture of F3;B ~- NCsHs. N 0.125 |0,320]0.379 | 0.254| 0.079 
The crystals were protected against moisture a Mas Wee pee ac Masts 
2 12210; 0.377]0.200] 0.2 
while the x-ray photographs were taken by cover- Cs 0.116] 0.177]0.374 | 0.141 0.873 
ing with a thin layer of paraffin. The results of de- Ce ee eee ode tas neta 
eet f 5 . . 00. Ue . 
termining the crystal class, the lattice parameters, 0.130] 0.427 | 0,385 | 0.361 | 0.093 
: , F, 0.175 | 0.487 | 0.430} 0.421 | 0.902 
the number of particles per unit cell, the x-ray and F, 0175 10.479|0.43010.4131 0.292 
pycnometric densities, and the space group of crys- Ps 0.045 | 0.409 | 0,300 | 0.343 | 0.097 


tals of F3;B ~-NC;H; are cited in Table 2. 

This structure is complex to decipher because 
the general position 4(e) of the space group Cin- 

— P2,/c is doubly occupied by atoms of the molecule. 
‘It is therefore necessary to determine the param- 
eters of 60 atoms (excluding hydrogen atoms). How- 
ever, the determination of this structure is alsoes- 
sential for the extension to the pyridine ring of the 
previously established shortening of the intermole- 
cular radius of the hydrogen in the benzene ring to 
0.74 A in the structure of [CITICgHs]Cl. It can be 
postulated that this effect facilitates the close simi- 
larity of the b axis of the monoclinic unit cells of 
[CITIC,H,]Cl (b = 5.77 A) and F;B —NC;H, (b = 

= 5.89. Al- 

The direct method of statistical equalities, 
which we have used before [6], was used to deter- 
mine the signs of the structure amplitudes. The 
methodology for the decipherment of the F3B <-NC;H; 
structure was described at the Fifth All-Union Con- 
ference on X-Ray Crystallography. Two electron- 
density projections were constructed on the planes 
(010) and (100). Further refinement of the struc- 
ture was carried out from sections of the three-di- 
mensional electron density constructed for k= 0,1, 
2, and 3. To construct the three-dimensional series, 
890 experimental structure amplitudes were used. 
These were taken from films of crystals obtained 
on the RGNS x-ray goniometer constructed by M.M. 
Umanskii. 

The coordinates of the atoms of two molecules 
are given in Table 3. Two projections of the inter- 
atomic bonds in the compound F3B<-NC;H; are shown 
in Figs. 1 and 2. The positions of the molecules in 
the unit cell are shown in Fig. 3 in a projection on 
the (010) plane. The average interatomic distances 
in the structure are: 


The precision of determination of the interatomic 
distances is 0.01-0.02 A , which can be seen from 
the coincidence of the results of the independent de- 
termination of interatomic distances in two mole- 
cules. The bond angles are: { FBF = 109.5°, ¥ NBF 
= 108.5°; { CyNCy = 117.5°; { NC,C3 = 123°; { CyC3C, 
= 19973) CoC. Catal onan: 

In a much earlier electron-diffraction study of 
the pyridine molecule by Schomaker and Pauling, 
the bond angles were assumed to be approximately 
120°. The deformation of the bond angles in the py- 
ridine molecule which has now been established for 
F3B~-NC;H; is in accord with radiospectroscopic 
data on the structure of pyridine [7]. 

The large value of the structure amplitude with 
the index (400) is caused by superposition of the 
pyridine rings at approximately x = 4 and *%. Re- 
finement of the atomic coordinates shows that the 
planes of the pyridine rings of both molecules make 
an angle of 4° with the (100) plane. This circum- 
stance leads to the same values of the coordinates 
for the two molecules, because both intermolecu- 
lar F,—H, distances correspond to the sum of the F 
and H radii. The construction of sections of the 
electron density at different heights shows that the 
lines along which the atoms B,N,C, and H; are dis- 
tributed make an angle of 3° with the (010) plane. 
The closest intermolecular distances are: F4—CJ 
= 2.95 A; Ff-Cj= 2.96 A, and FJ-C¥= 3.04 A 
(primes and double primes are used to differentiate 
between the two systems of molecules). The vari- 
ation in these distances is determined by the B—F— 
—C bond angles, so that the values of the interatom- 
ic radii depend on mutual orientation of the coval- 
ent bonds in the molecules [5]. The plane of the 
pyridine ring is perpendicular to the plane (010). 
They are precluded from being angled because of 
the short distances FJ—C4, F§—CJ, and Fy —H;. 
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Fig. 2, 


An experimental value of the intermolecular 
fluorine radius (rp = 1.34 A) was calculated from 
the intermolecular distance FJ—Ff'. Belov and 
Bokii [8] give an ionic radius of fluorine rp = 1.33 
A, while Pauling gives an intermolecular radius rp 
= 158624: (0), 

A considerable shortening of the intermolecu- 
lar hydrogen atomic radius ry = 0.80 A occurs in 
the structure of F;B ~-NC;H,; from the H—H dis- 
tance of molecules from neighboring unit cells. This 
effect confirms calculations of the intermolecular 
radius ry = 0.74 A in the phenyl ring of [CITI- 
C,H;]Cl. The occurrence of this effect in dif- 
ferent structures indicates that the basic reasonfor 
the shortening of the intermolecular hydrogen radi- 
us is the covalent nature of the C—H bond [5] and 
not thermal vibrations of the atoms in the crystals, 
as was proposed by Cochrane. It is interesting that 
the hydrogen radius in the water molecule, as de- 
termined by nuclear resonance [10], is 0.78 A. 


V. ZVONKOVA 


The accumulation of a large quantity of data on 
crystallochemical intermolecular atomic radii 
shows the necessity for constructing new data tables 
of intermolecular distances containing maximum 
and minimum values of the radii, i.e., values of 
the radii along the direction of the covalent bond 
and at right angles to this. 

It is possible that the presence of H—H and 
F—F contacts in the structure cause the low melt- 
ing point (t = 45°). 

To judge from the interatomic distances and the 
energies of formation, the strength of the donor - 
acceptor bonds increases in the series F;,B~-NCCHsg, 
F3;B~-NH3, F3B*-NC;H;. This indicates that the 
asymmetry and elongation of the unshared pair of 
electrons on the nitrogen atom is much greater in 
pyridine than in methyl cyanide. 

The widely accepted explanation of the nature 
of the donor—acceptor bond is that the nitrogen 
atom donates electrons to the boron atom (i.e., 


N—B) to form a positively charged quadricovalent 
nitrogen atom and a negatively charged quadrico- 
valent boron atom. This is incorrect. It does not 
correspond to the crystallochemical electronega- 
tivity scale, according to which the nitrogen atom 
is more electronegative than the boron atom. Con- 
sequently, the reverse redistribution of charges on 
the boron and nitrogen atoms, B'‘! and N7°2, should 
occur on compound formation. 

A crystallochemical study of the influence of 
the unfilled bonding orbitals of the metal atom on 
the formation of complexes of mercury and thalli- 
um has led to the conclusion that the orbitals of 
these elements are unable to accept more electrons 
than are necessary to form normal covalent bonds 
in molecules [11, 12]. 

The donor-acceptor bond is different in nature 
from the covalent bond. To form a donor—acceptor 


Fig. 3. 
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bond, an unshared pair of electrons from an atom 
of the donor molecule is used. 

The important role of unshared pairs of elec- 
trons in the formation of molecules arises because 
these pairs take part in the hybridization of orbit- 
als and hence take up a definite direction in space. 
For example, the nitrogen atom has the following 


electron distribution: | 2|p|p| p |. The unshared 


pair of s* electrons has spherical symmetry. In 
forming the pyridine molecule sp* hybridization oc- 
curs to give three o bonds with a plane trigonal 
configuration and one 7 bond in the direction of the 
remaining p orbital. Two of the three o bonds are 
attached to carbon atoms, while the third is filled 
by the fifth nitrogen electron, in the pyridine mole- 
cule. Thus asymmetry and elongation of the un- 
shared pair of electrons on nitrogen arises in the 
pyridine molecule. This asymmetry and direction- 
ality of the unshared pair of electrons is a neces- 
sary condition for taking part in the formation of 
donor—acceptor bonds. The directionality of the 
nitrogen valences in the pyridine molecule is shown 
in Fig. 4, where the arrow indicates the direction 
of the unshared pair of electrons, the thin lines 
represent the directions of the 7 bonds, and the 
thick lines indicate o bonds. 

The overlap of the directional unshared pair of 
the nitrogen atom with the vacant orbital of the 
boron atom leads to the formation of a donor—ac- 
ceptor bond between the BF3 molecule and the NC;H, 
molecule. In fact, in the formation of the coordina- 
tion compound, the electronegativity is little 
changed because the unshared pair of electrons re- 
mains to a considerable degree with the nitrogen 
atom. This is seen from the x-ray structural data. 
The greatest changes occur in the BF; molecule: 
the B—F bond length increases from 1.30 A in the 
free molecule [13] to 1.41 A in the complex, and 
the stereochemistry around boron changes from 
planar to tetrahedral. The hybridization changes 
from sp” to sp® on forming the fourth bond to boron. 
The interatomic distances in the pyridine molecule 
in the BF; NC;H,; molecule change within the 
limits of 0.01-0. 02 A, but there is a small change 
in the electron density in the pyridine molecule. 

The results of the x-ray structural study ef- 
fectively supplement data on the energy of forma- 
tion of a compound, since thermodynamic measure- 
ments only give the overall effect without an analy- 
sis of the separate chemical bonds in the molecule. 
At the present time, the role of the unshared pair 
of electrons is being studied hy a variety of physi- 
cochemical methods. A calculation of the dipole 


Fig. 4. 


moment of water by Coulson [14] showed that the 
asymmetry of the unshared pairs of electrons on 
oxygen played the main part in creating the dipole 
moment. The measurements of the dipole moments 
of boron complexes with nitrogen showed that the 


polarization of the bond B_-N was absent — they can 
be explained in terms of the asymmetry of thé un- 
shared pair of electrons in nitrogen [15]. 

Crystallochemical studies have great value in 
elucidating the nature of donor—acceptor bonds. In 
a recent review on the nature of donor—acceptor 
bonds in complex compounds, Nyholm [16] noted 
that the understanding of the nature of the donor— 
acceptor bond in complex compounds, and the elu- 
cidation of factors which affect its stability, are of 
great importance in the stereochemistry, thermo- 
dynamics, and kinetics of complex compounds, but 
the theory has weak points and, in particular, it 
does not explain that many compounds which con- 
tain atoms with unshared pairs of electrons, e.g., 
halides, do not form coordination compounds with 
metal salts. 

However, crystallochemical theory can be used 
to explain these facts. As an example, we shalldis- 
cuss the reason why boron trifluoride does notform 
coordination compounds with organic molecules 
which contain halogen atoms. In the light of the dis- 
cussion, this indicates that the fluorine atom does 
not have the asymmetry and elongation of the un- 
shared pairs necessary to form donor—acceptor 
bonds between the molecules. 

A great deal of experimental crystallochemical 
material on directed covalency in molecules and in 
complex compounds is now available. The direc- 
tions of the covalencies in elements of the Vth, VIth, 
and VIIth groups of the periodic system can be il- 
lustrated by the structures of the molecules PCl,, 
TeCl,, and Cl1F3 (Fig. 5). The halogen atoms can 
form three covalent bonds in a T structure. To 
form the additional bonds, a p* pair of electrons is 
promoted from the ground state of chlorine: 


s°| p?|p?| p | The direction of the T bonds is such 


that the formation of the donor—acceptor bonds 
avoids steric hindrance. The effect of steric hind- 
rance appears to be stronger in donor-acceptor 


58 Ze NA ZVONKOYV & 


PCI, TeCl, CIF, 
Cl Cl F 
Cl Cl 
Cl P Te Cl Fi 
Cl Cl 
Cl Ci F 
Fig. 5. 
bonds than in covalent bonds [17]. The B<-N 


donor —acceptor bond is weaker than the B—N co- 
valent bond. This is seen from the dissociation 
energy of the compound, B(NMe,)3, which is 89.7 
keal [18], a compound in which the B—N covalent 
bond distance is 1.50 iN the sum of the covalent 
radii. 

Finally, it must be remembered that the chemi- 
cal bonds in molecules are usually intermediate in 
type and in their structures molecules display pre- 
dominantly characteristic features of different types 
of chemical bonds. 

Attention is now turned to the fact that three 
equal bonds are formed in the BF; molecule with 
sp” hybridization, while in the ClF; molecule the 
bonds have different strengths, the equatorial bond 
being stronger than the two axial bonds. This is 
evident from radiospectroscopic data on the mole- 
cular structure of CIF; [19], for the Cl—F bond 
lengths are 1.598 and 1.698 A, respectively, and 
the bond angle FCIF = 87°29', This is explained by 
the differences in energy of the 2s, 2p, and 3s 
electronic states. 

The following rule exists: As the filling of the 
outer electron orbitals occurs with a change from 
s* to p*, the directional character of the additional 
valence bonds changes. These are formed by: un- 
shared electron pairs taking part in hybridization 
of the electrons in the molecule. Hence, the sta- 
bility of the donor—acceptor bonds in complex com- 
pounds decreases in order of donor atoms: nitro- 
gen, oxygen, fluorine, and their electronic analogs. 

Since nitrogen and chlorine atoms have the 
same directions for separate valence bonds perpen- 
dicular to the equatorial plane, they should have 
some properties in common when the atoms are in 
combination. This is notobserved in the formation of 
donor-acceptor bonds between molecules, butitis ob- 
served in intramolecular covalent bonds. In deriva- 
tives of benzene in which hydrogen is replaced by a 
chlorine atom or an amino group, shortening of the in- 
teratomic distances is observed, Cl—C to 1.706 A 
[20] and N-C to 1.37 A [21] in comparison with 
Cl-C = 1.78 A and N-C = 1.47 A in aliphatic com- 


pounds. Thus, the C—Cl and C—N bonds are 
strengthened because the unshared pairs take part 
in hybridization. The sp” electrons take part to 
some extent in the three o bonds at the nitrogen 


H 
atom in the group wn — cg because,x-ray data 


show,that the group is planar. The effect of the un- 
shared pairs of electrons on the nitrogen and chlo- 
rine atoms is in the formation of sp hybrid orbitals 
directed at right angles to the equatorial plane and 
parallel to the p orbitals of the carbon atoms of the 
benzene ring, so that they take part in 7 -interac- 
tion with the benzene ring. 

In crystallochemistry, the difference in the 
contributions of the s and p electrons to hybridiza- 
tion is estimated from the bond angles in molecules; 
for example, the larger size of the bond angles at 
nitrogen in comparison with those at bismuth is 
connected with the greater participation of the ni- 
trogen s* electrons in the general hybridization. It 
can be confirmed from an analysis of the internu- 
clear distances that the effect of the unshared pair 
of electrons on nitrogen is much greater in the py- 
ridine molecule than in the nitrile group. 

Crystallochemical studies of interatomic dis- 
tances are an important additional criterion in de- 
ciding the question of the role of the unshared pair 
of electrons in the formation of chemical bonds in 
molecules and chemical compounds. 

The author expresses her deep thanks to G.S. 
Zhdanov for discussing this work. Thanks are also 
expressed to Z.P. Lipina, A.N. Khvatkina, and 
A.N. Abramov for help in taking the x-ray photo- 
graphs, preparing the crystals, and determining 
the pycnometric density. 
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Ionic refraction determinations were made on ions which did not have noble- 
gas structures, using experimental data on crystalline materials containing 
the ions. The determinations were carried out by comparing the molar re- 
fractions of salts containing complex anions, since the additive law is most 
closely followed under these conditions. The ionic refraction values obtained 
were Zn*t = 0,6; Ni?+ = 0.6; Co*t = 0.7; Fe?* = 1.15; Cu?+ = 1,2; Mn** = 1.4; 
Pb** = 9.0; Ga®t= 0.5; Cr’? = 2.1; Fe®* = 3.4 cm; these agreed with values 
calculated by independent methods in all cases. 


Refraction data are used widely in chemistry 
to solve the most diverse problems, included 
among which are structural questions. In the or- 
ganic field, the refractometer method has been ap- 
plied to structural chemistry studies for a consider- 
able time [1], but its application to inorganic chem- 
istry is comparatively recent. Refractometer re- 
sults were successfully used to determine ionic 
radii by Wasastjerna [2] and, in particular, by 
Kordes [38,4], and have been applied by us [5-11] to 
the determination of geometric configurations of 
complex compounds, and to the study of hydrogen 
bonding in inorganic crystals. 

The refractometric method has some decided 
advantages over other structural chemistry me- 
thods; it provides a rapid result, and does not need 
a high-quality specimen. It therefore appeared of 
interest to extend the use of the refractometer me- 
thod to the most widespread and practically import- 
ant class of inorganic substances, the silicates. 
One of the biggest obstacles to this was the lack of 
data on refraction in the crystalline state for many 
of the nonnoble-gas type cations met with in sili- 
cates. The present investigation included a deter- 
mination of the refractions of a series of di- and 
trivalent cations, first among which were iron and 
manganese. 

At the present time there are many different 
systems of ionic refraction determination but, of 
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these, only the ones described by Fajans [12], 
Pauling [13], and Kordes [3] are complete and con- 
sistent. Fajans derived his ionic refraction values 
from experimental data on aqueous solutions of 
electrolytes, Pauling used the relationship between 
refraction and the secondary Stark effect, and 
Kordes worked from the relationship between ionic 
refraction and ionic radii. In principle, none of 
these three methods can give characteristic refrac- 
tion values for ions not of the noble-gas type. In 
fact, since there is no rigorous quantitative theory 
relating the refractions of these ions and their spec- 
tral properties and dimensions, the Pauling and 
Kordes methods cannot, in general, be used to find 
the refractions of nonnoble-gas ions accurately. 
The Fajans method, based on the principle of addi- 
tive superposition of refractions of infinitely diluted 
aqueous electrolyte solutions, also cannot lead to 
the desired result, since the method gives only the 
apparent ionic refractions. In aqueous solutions, the 
ions exert a distorting influence on the water mole- 
cules surrounding them, decreasing their refrac- 
tion. The magnitude of this decrease can be esti- 
mated for noble-gas type ions, since it is known for 
the sodium ion. For ions not of the noble-gas type, 
this starting point does not exist, and so data on so- 
lutions cannot be used to determine refractions of 
ions in the gaseous and crystalline states (these 
states being equivalent where the cations are con- 
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cerned [7]). In Fajan's opinion, data on character- 
istic refractions of ions of this type can be obtained 
only from measurements on crystals. 

1. Not long ago, some American workers [14] 
attempted to determine the polarizability of the iron 
ion using data on crystalline iron halides. However, 
when these authors subtracted the polarizability of 
a halogen ion, derived from results from alkaline 
salts, from the molar polarizability of the ironhal- 
ide, they arrived at a physically meaningless re- 
sult — the ionic polarizability turned out to be nega- 
tive. It appeared that this result could only be due 
to an incorrect method of calculation, the error 
arising in not taking into account the considerable 
reduction of the halogen ion refractions in the field 
of the iron ions. To avoid this error, it is obvious- 
ly necessary, first, to use anionic refractions cal- 
culated from results on salts with cations of iden- 
tical charge and, second, to use complex rather: 
than simple anions, since the former are more dif- 
ficult to distort and, therefore, obey the additive 
law better. These conditions are approached most 
closely in the case of isomorphic salts containing 
chemically identical anions. 

A comparison of the molar refractions of iso- 
morphic salts gives the difference in cationic re- 
fractions, just as a comparison of refractions of 
neighboring members of a homologous series gives 
the value of the homologous difference refraction, 
this being the basic method of determining atomic 
refractions in organic chemistry. Table 1 lists the 
refractions of a series of calcium and lead salts, 
worked out from results given in Winchell's hand- 
book [15]. 

From the table it can be seen that the differ- 
ence between the refraction of the lead and calcium 
salts varies around 7.94 cm®. On adding the 
refraction of the calcium ion (1.40 cm® according to 
Fajans [12]), we obtain an ionic refraction value for 


TABLE 1 
Compound Rp, cm AR, cm? 
PbCO 20.05 
CaCO; 12.61 7.44 
PbSO 22.14 
CaSO, 15.83 6.31 
PbWO 34.54 
CaWOg 22°33 9.18 
PbMo0, 32.05 Pen 
CaMoO 22.54 = 
Pb-Si0, 37.08 
Ca,Si0, 20.92 8.08 
Pb; (PO4)s 56.44 ie 
Ca (PO4)a 35.08 : 


Pb** of 9.34 cm, According to independent meas- 
urements made by Kordes [16] on lead glass, the 
ionic refraction of lead is 9.13 cm*®. Thus, a com- 
parison of the refractions of analogous complex 
salts can give quite accurate results for the refrac- 
tions of ions not possessing noble-gas electron 
shells. 

Table 2 shows a comparison between the molar 
refractions of iron and magnesium salts, calculated 
from the most reliable experimental data. For all 
the substances only one refractive index is given, 
averaged by the formula n = (ig + Ny ° Np) 3 this 
being necessary for calculation of refractions by 
the Lorentz — Lorenz formula. 

From Table 2 it can be seen that the difference 
between the molar refractions of analogous ironand 
magnesium salts averages 2.27 cm®. The additive 
rule is obeyed in this case quite closely, since the 
spread of results lies within the range of experi- 
mental error (£0.10-0.15 em’). Since, for increas- 
ing molecular weights, errors in measurement of 
density have an increasing effect on calculations of 
molar volumes, and hence also on refractions, we 
cut short the Tutton salt series at the rubidium salts 
for the sulfates, and at the calcium salts for the 
selenates. We acted similarly in the comparisons 
shown in the other tables. 

If the complex anions had undergone no change 
in refraction at all on exchanging the cations, then 
to calculate the refraction it would be sufficient to 
add the refraction of the magnesium ion (0.28) onto 
the difference in refractions found (2.27). However, 
the polarizabilities of the complex ions do in fact 
change, although this change may be quite insigni- 


TABLE 2 
Compound = i) Rp,cm AR,cm> 

OC a a ee ee 29,33[27] [4.794 [27] | 12.48 |-o gy 
MEO 6S cate alo ole Ae AOE owe {Pq} |) O98) || 
ANSON WENO! 6 aS oo 146. 5O(?8] | 1.4784[?9] | 41.52 | 5 13 
MSO, (LO | es te 146 .80(29] | 1.449215] | 39.39 | * 
NESSIOA oo 6 6 00 6 6 47 .16[?1] | 4.854 35] | 20.76 2.38 
MeSiO ges 6 hoes 43.75(22] | 4.652 [25}| 16.00 | “ 
K Fe (SO4)2:6H2O0 - | 198.05[7?] | 1,4850[?2] | 56.76 | 5 15 
K2Mg (SOq)2-6H,0 196 .58[??] | 1.4664[?2] | 54.64 ; 
(NHa)dFe (SO4)2-6H,0 . . | 208.86{?2] | 1,4925[%2] | 60.66 | 5 99 
(NH4)2Mg’SOq)2-6H,O . . | 207.78[2] | 1.4744[?2] | 58.46 | “ 
RbpFe (SO4).-6H20 . . . | 207.84 22] | 1.4889[72]| 59.98 | 4 3, 
RbyMgSOq)2-6H20 . . . | 206.18[22] | 1.4713[22]] 57,68 | 
K2Fe (SeO4)2-6H20 . . . | 210.39[°9] | 1.5207[3] | 64.03 | 5 35 
K.Mg (SeO,4)2-6H20 . . . | 208.63[?9] | 1.5033[23] | 641.70 | ~ 
(NH,)sFe (SeO4)-6H,0 «| 220,39(25] | 1.5292[°9] | 68,00 | 5 95 
(NH4)sMg(SeOq)o-6H,O . | 249.4223] |4.5411(29]] 65.75 | * 
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ficant. Because of its small volume, the magnesi- 
um ion exerts a strong polarizing effect on its 
partners, and it therefore is possible that part of 
the difference AR is due to a reduction in anionic 
refraction. To check this possibility, we must 
compare the molar refractions of iron salts with 
the refractions of analogous calcium and strontium 
compounds, since, in the latter, we are dealing 
with underformed anions [24]. Table 3 shows the 
corresponding molar refractions, wholly calculated 
from data given in Winchell [15]. 

From the information shown in Table 8, it fol- 
lows that the difference in refraction between the 
Cat and Fe2* ions is 0.17 cm’, and that between 
Sr*‘ and Fe** is 1.44 cm’, Using the ionic refrac- 
tion values quoted by Fajans (Ca2* = 1.40 cm? and 
Sr = 2.58 em®), we obtain the following refraction 
values for divalent iron: from the Ca salts, 1.23 
em?; and from the Sr salts, 1.14 cm®, giving a 
weighted mean value of Rpg2+= 1.20 cm®. This last 
value differs considerably from that obtained from 
the comparison of iron and magnesium salt refrac- 
tions. The disagreement is due to the lowered re- 
fraction of the anions in the magnesium salts be- 
cause of the strong polarizing influence of the Me2* 
ion. This effect can be allowed for by comparing 
the molar refractions of the magnesium salts with 
those of the analogous Ca and Sr compounds. 

The difference between the ionic refractions of 
Ca2t and Mg?t is equal to 1.40—0.28 = 1.12, while 
the difference between the molar refractions of 
identical anions is equal to 2.54 cm®, as shown by 
Table 4. Since the cations do not alter their refrac- 
tions, the difference of 1.42 cm? (2.54 —1.12) is 
the drop in anionic refraction due to the field of the 
magnesium ion. A similar value is obtained from a 
comparison of magnesium and strontium salts: 
Rgrx —Ryex = 2-30 cm’, which means that the re- 
fraction drop is equal to 3. 88—2.30 = 1.58 cm’®, 
From this, it can be concluded that a reduction in 
anionic refraction takes place in magnesium salts, 
averaging 1.47 cm*®. When this correction is taken 
into account, the results in Table 2 give a value for 
the refraction of the Fe’* ion which agrees well with 


TABLE 3, 

Compound Rp,cmP| aR,cm? 
GaGOsas, bane 12.61 0,18 
HeCOsgy ns 12.43 
GagSiOae. ese 20,92 0.16 
Fe,Si0;. .. .| 20.76 
ehGans wt bette con BT laced 
WeCOy . . . cl td.ad 


that obtained from Table 3 (1.20 cm’). In fact, if 
we subtract from the difference ARfe-mg = 2-27 
cm? the quantity 1.47 cm® (the "magnesium drop"), 
then, on adding the Mgt ion refraction of 0.28 cm®, 
we obtain a refraction value for the Fe** ion of 1.08 
em®. This leads to a final result of Rpett= 1.15cm?, 

Thus, the ionic refraction of iron, worked out 
from a comparison of iron and magnesium salts (al- 
lowing for the refraction-drop correction) agrees, 
within the limits of experimental error, with the 
value obtained from a comparison of iron salts with 
calcium or strontium salts. This circumstance al- 
lows us to make future comparisons using only Mg 
salts, which is very convenient, since they are iso- 
morphic with Zn, Mn, Co, Ni, and Cu salts, and 
have been studied much more extensively than any 
other salts containing divalent cations. _ 

The validity of this method of calculating ionic 
refractions can be checked through a study of zinc 
salts, since the ionic refraction of Zn2+ has been 
determined by a number of authors in various ways. 
Table 5 shows a comparison of the molar refrac- 
tions of analogous zinc and magnesium salts. 

The average difference between these refrac- 
tions of zinc and magnesium salts is 1.82. From 
this, the refraction of the Zn2* ion is equal to 0.63 
cm?®, According to Kruis [27], the refraction of the 
Zn" ion is equal to 0.60 cm®. Kordes, working 
from a relationship between ionic refractions and 


radii [3] gives a value of 0.71 cm®. So our value of 


the zinc ion refraction agrees with values obtained 
by completely different methods, and the method of 


calculation we have described can be applied to the 
determination of the refraction of other ions. 

Table 6 lists experimental data from which the 
molar refractions of copper, nickel, cobalt, and 
manganese salts have been calculated. 

From theseresults (Table 6) we can determine 
the ionic refractions of Cu2", Ni2*, Co2t, and Mn". 
To do this, we need only to subtract, from the 
molar refractions of the metal salts listed, the 
molar refraction of their magnesium analog, and 
then from this difference subtract a further 1.19 
em? (i.e., subtract 1.47 cm?, the "magnesium 


TABLE 4. 

Compound Rp,cn? AR,cm? 
Ca COsP co oa te 12.61 2.62 
MoCOS mis aa 9.99 
GasSiOg a sak o 20.92 2.46 
Mg.SiOg 16,00 
SrCOg tose. 13.87 3.88 
MeGOy rt. . 9.99 
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TABLE 5 
M a 3 
Compound ak. nD Rp,cm | AR,cm 
EARS kOe Ge bey beh eS 28.34 [35] 4.772[}5] 11.79 1.80 
Mo COs oe teenie € 74 cach 28 . 30 624 9.99 
CoS. TH.Oe sn. 143. 78[21] 4738 [25 | 44.37 1.98 
MESO THO: roan o pans oc 146.80 4492 39.39 
K52n (G0,),-6H.0.°". 196.1628] .4859[25] | 56.32 1.74 
KoMg (S0)),,6bs0 a po, 196 58 4664 54.61 
(NH4)oZn (SOq)o-GH,O . . . | 206. 38(26 .4937(25[ | 60.06 1,60 
(NH4)oMg (SO4)2-6H30 207.78 4T44 58.46 
Rb.Zn (SO4)o-6HO . . . «| 205.58 [25] .4897(2°] | 59.44 4-73 
RboMg (SO4)2:6H.O . . . . | 206.18 4743 57.68 
KZn (SeO4)o-6H,O . . . . | 208.80[25] 5212/26} | 63.62 1,92 
K5Mg (SeO4).-6H290 . . . . | 208.63 5033 61.70 
(NH4)oZn (SeO4).-6H2O . . | 247.73[26] 5308[2*] | 67.3 1.61 
(NH,4)2Mg (SeO4)o-GH20 . . . | 219.42 S114 65.75 
TABLE 6 character. Since, however, the compounds met 
oe 7 ‘ with in nature often include, say, trivalent iron, it 
iSite R 
Compound d tie de would be useful to know how to calculate the propor- 
eo aie heey 31.23[5] 1.7391] 12.58 tion of the molar polarizability of the whole com- 
MisclOpe ret sane e 49 _98[22[ 1,7967[22] | 21,29 pound which was due to the iron. As well as the re- 
Se honing Mae K pr eit ‘a fraction of trivalent iron, we also worked out the 
N . Pham wake 4. z ¢ oF -9 . 5 : * ° . 
(NH,),Mn (Se0,)»-6H,0 223.35[2| | 4.5217[°9] | 68.09 ionic refractions of trivalent chromium and galli- 
NGO en ee ae 29 .01]}5] 1.766[35] 12.00 um, to see how far we would be justified in apply- 
cy lesa eat Gene cath er ing the method to the calculation of the refraction of 
K,Co(SOq)2:6H,O. . . 68 [22 41.4892[?8] | 56. ee : : 
(NH,)sCo (S0,)2-6H,0 . . | 206.40[2%| 1.4962(281 | 60.34 positively charged trivalent ions. 
Rb,Co (SOq)2-6H,O . . .| 205.03[2] | 4.4930[29] | 59.57 In Table 8 are compared the molar refractions 
K,Co (SeO4)2:6H,O . . | 208.60/?9 1.5252[23] | 63.95 of aluminum and iron (ferri) compounds, only the 
(NH4)2Co (SeOg)g-6H,O =| 218.40[] | 1.5335]%] | 67,74 most suitable for comparison being chosen, i.e., 
WiSO HIB) Gen: OAL. 14.40[7] 1.4836(°] | 44,18 the al includ ly th tash and i 
K,Ni (SO4)o-6H,O. . . . | 193.99[°2] | 1.4934(28] | 56.44 SERS Dae OS Oe Le ate eee ee Soe 
(NHq)2Ni (SOq)2-6H,O . | 203.91[22| 4.5012(28] | 60.09 members. Rubidium, cesium, and thallium alums 
RbpNi (SO4)2-6H20 . . . | 203.43(?*] 1.4969[°*] | 59.52 have such high molecular weights that, for an error 
i(8 23 2 23 4 ope ° * 

Po ieee heed piace iaeaee a eee in specific gravity of +0.02, the error in calculat- 
i -OH,O . i < : : is A . 
Cen cas 14.83(2) | 1.48[23] 42.14 ing the refractions would be more than a cubic cen- 

KyCu (SO4)2:6H,O. . . . | 196.49[??] 41.4907[22] | 56.88 timeter. 

ea area tel _% oe pa eats ete On the average, the molar refractions of analo- 
Rb,Cu Ae OH sOo,. St f ; : E , ? 3 0 A 

K,Cu (SeO,),-6H,0 . . | 209.66(29] 1.5226(29] | 64.04 gous iron and aluminum compounds differ by 4.98 
(NH,)2Cu (SeO,4)2: 6H,O 220.64[23] 4.5311(23] | 68.27 cm’. 


drop," and add 0.28 em’, the magnesium ion re- 
fraction). Table 7 lists the refraction differences 
Rex— Rvigx and the refractions calculated for the 
cations mentioned above. 

2. This method of determining ionic refrac- 
tions from a comparison of the molar refractions 
of complex compounds can be extended to cover tri- 
valent cations. Of course, the concept of ionic re- 
fraction cannot strictly be applied to highly charged 
atoms, since they form bonds of partially covalent 


To determine the refraction of trivalent iron, 
we need only to calculate what part of this differ- 
ence AR is due to the drop in anionic refraction be- 
cause of the strong distorting influence of the alumi- 
num. This problem can be solved by comparing the 
refractions of identical anions in aluminum and cal- 
cium salts. Table 9 shows such a comparison, with 
values of aluminum and calcium ion refractions 
taken from Fajans (0.17 and 1.40 cm’, respectively) 

It can be seen here that the aluminum ion re- 
duces the refraction of its partners by 1.7 cm’, 
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TABLE 7 
E= element 
S-ompound Cu Co Ni Mn 
BECO ST cere he Ge ee aka = 2.04 — 2,99 
EeSiOGES AR ea ee eee - - 2.64 
ESOre Hs OF roo eee eee 2.75 Uae I ate AS, -— 
SSO) SO 3 8 56 og a Se 2.27 1.88 | 1.80 -— 
(NHa)oE (SOy)2-6H2O . . . . 2.06 1.85 | 1.63 | 2.34 
RbsE (SO,4)2-6H2O. . . 2... 2.40 M089) | W845 162,91 
KE (SeQs)s-6H, On ae ees 2.31 ye Pay \\ ah Aghl _ 
(NH,)2E(SeOa)2-6H,O . . . . . 2.52 eGR) | AI ee) 1) ache 
AN ths ESIC 6 Go 6 SB 2.39 USS |) etek Peas 
Reg 1. ee hale ee Ai 0.7 0.6 1.4 
TABLE 8 
Compound Rp,cm?® AR, cm? 
Ke (SOq)2:12H20 78.35[22] 5.00 
K Al (SOq)o-12H 30 73.35[22| 
NII Fe (SOq)o-412H0 80.35]22] Lei? 
NH4AL (SO4)9-42H 90 75.63|22] 
HOMO cA MIO. Gs ok oF 29 .65[79] 9.22 
AIPO,:-2H,O . . 20,2329] 
TABLE 9 
M . 3 : 3 
Compound i | @p Rp,cm| R.anion | Ar,cm 
CHS 6 a eal) AO MES| yf Hobiskoe| 15,59 14.19 hate 
Alo (SO4)3_. «| 126.3[8) | 1.510% 37.75 12.47 
Ca3 (PO), . .| 107.2)15] | 1.628[25] | 35.50 15.65 1.67 
AIPO,. . . .| 46.19[5] | 1.525)5] | 14.45 13.98 


*From our measurements, 


Thus, the difference in refraction between the Fe®* 
and Al** ions works out at 4.98 —1.70 = 3.28 cm?®. 
Adding the refraction of the aluminum ion, we find 
that Rygst = 3.45 cm?, 

In a Similar manner, the refractions of trival- 
ent chromium and gallium can be calculated. These 
calculations are shown in Table 10. 

According to Kordes [3] and Pauling [13], the 
refraction of the Ga*+ ion is equal to 0.49 or 0.50 
cem®, respectively, which is almost exactly the 
same as our value. So the refractions of trivalent 
cations can be determined quite accurately from a 
comparison of the molar refractions of their cry- 
stalline salts. 

3. In conclusion, we will compare (Table 11) 
the characteristic refraction values obtained by us 
for nonnoble-gas ions, with their values in aqueous 
solutions, determined by Fajans and Joos [81]. 

This comparison of the refractions of "crystal- 
line" and "dissolved" ions shows that the pairs of 


TABLE 10 
Compound | Rp,cm® arom AR | RRis 
KCr ($04)o-12H,0 77.02(%] | 3 67 
KAI (SOq)2:12H,O 73.35 : 
3.68 2.15 
NH,Cr (SO4).:12H,O 79 .32[?] 3.69 
NH, AI (SOq4)2-12H,O 75.63 s 
KGa(SOq}o-12H,0 . . .| 75.441) | 5 og 
KAI (SO4)2:42H,O Seat Wey) : 
2.06 0.52 
NH,Ga (SO4)s-12H,0 77.6771 | 5 oy 
NH,A1 (SO4).-12H2O oOo ; 


values are, in general, similar; the refractions of 
the divalent ions, nickel, iron, manganese, and 
lead are the same, within the limits of experiment- 
al accuracy, the refractions of the Zn", Cu2+, and 
Fe* ions are less in solution, and only with Co2* is 
the refraction of the "dissolved" ion significantly 
greater than that of the "crystalline" ion. 

Since our ionic refractions coincide in two 
known cases (see above) with the "gaseous" ionic 
refractions, we can deduce that some of the differ- 
ences in comparison to "dissolved" ions are due to 
the influence of the water. We know that the mole- 
cules of water around cations have reduced refrac- 
tions. On subtracting the refraction of the water 
from the general refraction of a solution, this fixed 
refraction drop is transferred to the cations. This 
effect operates in the case of divalent zinc and cop- 
per, and for trivalent iron. The discrepancy in the 
case of the cobalt refraction is apparently due to an 
inaccurate determination of its ionic refraction in 
solution, since divalent nickel and cobalt ions are 
very similar in all their properties, and it is im- 
possible to understand why the polarizability of co- 
balt in solution should be twice that of nickel. 

It is an interesting point that the refractions of 
nonnoble-gas ions increase with increase in their 
positive valences. Thus, the refraction of Cu!* 
= 1.08, according to Kodes' data [3], while that of 
Cyt 2 cm?®; the refraction of Fe2t= ier and 
that of Fe**= 3.4 cm®. The increase in refraction 
of the cations is due to a shift in the nature of the 
chemical bonding, from ionic to covalent, because 
of the increase in polarizing action of the nonnoble- 
gas ions as their charge increases. Since the ioni- 
zation of atoms to form cations is accompanied by 
a loss of electrons, and ionization to form anions 
involves a gain in electrons, it is clear that on go- 
ing from an ionic bond to a covalent one the elec- 
tronic polarizabilities (refractions) of cations will 
increase, while those of anions will fall. 
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TABLE 11 
Ion 
Rp 
Znt2 | Nit? | Cot2 Fet? | Cut? | Mnt? Pbt2 | om Crt3 Fets 
| 
Crystal OnGE KORG On Tiel tied 5 2s diez 9.0 On|) Dall Bee 
Solution @.29 || 0.67'| 1.28 | 4.22 | 0-39 14.37, | 9.38[2) | — | OEE 
This change in the type of chemical bonding im- 9. G. B. Bokii and S. S. Batsanov, "The refrac- 


poses a natural limitation on the application of the 
very concept of ionic refraction and, consequently, 
on methods of determining it. The method described 
in the present article can therefore only be applied 
to di- and trivalent ions. 

The method of determining refraction values of 
nonnoble-gas cations derived in the present work 
supplements the available systems of finding ionic 
refractions in the crystalline state, and at the same 
time extends the area of application of refractom- 
etric methods in inorganic chemistry. 
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It is established empirically that the fundamental principles of crystallography 
with regard to constancy of the angles between the faces and the integral nature 
of the ratios of the parameters of the faces are equivalent to the principles of 

chemistry with regard to constancy of composition and multiple proportions by 


weight, these being also a consequence of the atomic structure of matter. 


It 


is shown that constancy and multiple proportions also appear in the actual 
building units of the crystal, i.e., the atoms and ions. 


The fundamental principles of crystallography, 
which were discovered empirically, are the law of 
constancy of the angles between the crystal faces of 
a given substance (Steno, Rome de Lisle), and the 
law of rational proportions, expressing the integral 
nature of the ratios of the parameters of the faces 
in axial units (Hatiy) [1]. These two principles ap- 
pear to be completely equivalent to the two funda- 
mental laws of chemistry, namely, constancy ofthe 
chemical individual (Proust), and the principle of 
multiple proportions by weight (Dalton). All these 
generalizations are essentially the logical conse- 


quence of the real existence of atoms and molecules. 


Outside the atomic theory, all these laws would be 
obscure, incomprehensible, and would be deprived 
not only of inner unity, but also of physical mean- 
ing. In different ways, they express the same 
thing: the discrete nature of matter; chemistry in 
terms of mass, crystallography in terms of spatial 
forms. This unity of crystallography and chemistry 
can be shown in various ways and hardly calls for 
further examination; in view, however, of its par- 
ticularly convincing nature, one of these ways de- 
serves special mention; this is the Shubnikov law 
[2], establishing a connection between the chemical 
formula expressing composition and the crystallo- 
graphic characteristics of a compound. 

Although the foregoing general principles of 
crystallography and chemistry are quite important 
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and reveal the specific chemistry of the atoms most 
distinctly, they can by no means claim to be rules 
which have no exceptions. On the contrary, as 
shown by Kurnakov [8], compounds of constant and 
variable composition are possible, i.e., daltonides 
and berthollides, the latter, while they do not pos- 
sess a Singular point in the limits of the representa- 
tion of the homogeneous individual phase in compo- 
sition—property diagrams, are nonetheless equally 
chemical individuals with the daltonides. Crystal 
chemistry indeed provides an interesting approach 
to the understanding of this problem as a whole. The 
point is that, for an isolated molecule (gas), the 
principles of constant composition and multiple pro- 
portions know of no exceptions. There are also al- 
most no exceptions for typical ionic lattices, in 
which the interchangeability of differently charged 
ions is impossible on the basis of electrostatic con- 
siderations alone. On the other hand, it is a priori 
probable, in the case of covalent and intermetallic 
compounds, in which it is relatively an easy matter 
to effect variation of the components, and, therefore, 
precisely for such crystals, the propagation of indi- 
viduals of variable composition is quite natural. 

As in chemistry, departures from the observ- 
ance of the laws of constancy and multiple propor- 
tions are well known and have been explained, so in 
crystallography it is not difficult to predict such de- 
partures from its fundamental laws, and since these 
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very laws are the consequence of the atomic nature 
of matter, the deviations from these laws, which 
constitute an over-abstraction of the reality, are 
also due to the atomic structure of crystals. 

In fact, crystal faces are far from being plane 
structures, but possess "irregularities," whichare 
by no means negligibly small. Tolanski [4] has 
demonstrated experimentally the existence of a re- 
lief exceeding 1000 A on some crystals, which, 
figuratively speaking, constitutes almost the same 
"irregularity" as does a mountain 1000 m highcom- 
pared with the height of a man. The resulting rela- 
tive constancy of the angles of a crystal, and also 
the repetition of the parameters of the faces, are 
obvious. These generalizations have thus rather a 
limited character; they are statistically true onthe 
average. This,of course, is not to depreciate their 
functions as generalizations which, in themselves, 
reveal the essential features of the solid state. 

The question naturally arises that, if the prin- 
ciples of constancy and multiple proportions occur 
in chemistry as a consequence of the discrete na- 
ture of matter, are they not also applicable in some 
degree to the very building units of crystal chem- 
istry, the atoms and ions. The special examination 
of this question is of considerable interest since, 
as applied to the masses of atoms, it was raised 
more than a century ago by Proust (formation of 
complex atoms by the "polymerization" of hydrogen 
atoms); at that time, it was answered in the nega- 
tive, in view of the clearly nonintegral nature of 
atomic masses, but now, following the discovery of 
isotopes,ithas been answered in the affirmative, 
since the atomic weight of any isotope is an integral 
repetition of the weight of the nucleus. 

By following this simple analogy, it might be 
expected that the principle of multiple proportions, 
as applied to individual atoms and ions, ought to re- 
sult in the multiple relationship of the dimensions 
of all atoms (or ions) to the dimension of the par- 
ticle appearing as the smallest. It is well known, 
however, that such expectation is by no means justi- 
fied. 

There is, however, another possible approach 
to the solution of this problem, which is based on a 
consideration of the comparative "thickness" of the 
various electron levels (K, M,N, L,O,P, etc.), and 
which is particularly simple for particles having 
spherical symmetry, since, in such cases, the 
question is always one of the distance correspond- 
ing to two successive minima of radial electron den- 
sity. In principle, such minima can be calculated, 
a definite solution of the problem being obtained, 
both from the geometrical point of view, and that of 


theoretical physics. In fact, Odiot [5], making use 
of the Hartree function, obtained quite good results 
for a number of chemical elements, close to the re- 
sults obtained by calculation of the radial density 
minima. 

These methods, however, are somewhat compli- 
cated, and the results are comparatively limited. It 
will be shown below that there is another possible 
approach to the same problem, the theoretical basis 
of which will be considered later. Following this 
method, it is possible to show convincingly that the 
principle of multiple proportions actually applies al- 
so to atomic and ionic dimensions. 

We commence with an examination of those par- 
ticles which have been studied most reliably, and 
which are closest to the spherical form in their 
electronic "architecture." This requirement is sa- 
tisfied by ionic radii of the rare-gas type. More 
particularly, we choose to consider the cations. 

The electrons of any ion (and atom) form two 
groups, sharply distinguished in character, the 
electrons of the K-layer, closest to the nucleus, 
(two-electron group), and all the others (eight- 
electron environment in the case of ions of rare- 
gas type). 

The electrons of the K-layer, which are clos- 
er than the others to the nucleus, form a sphere, 
the dimensions of which depend on the number of 
charges of the ion. This may be shown as follows. 

An examination of the more reliable system of 
ionic radii, the well-known Goldschmidt system [6], 
already shows that Me--- ions are quite rare, Me™ 
ions actually do not exist in lattices, and the estim- 
ate of their dimensions is only tentative. Ions of 
period II of the periodic table are examples of rare- 
gas type ions having a single K-shell. The largest 
is the lithium ion and the smallest the nitrogen ion 
(0.1 A). Confining ourselves in all subsequent cal- 
culations to an accuracy of 0.1 A, we may practi- 
cally assume the radius of the last-mentioned ion 
to be zero. 

The dimension of the K-shell decreases linear- 
ly from the lithium ion with increase in size of the 
ion charge 7». Taking the foregoing into account, 
we have 


K =0.2(5—»). @) 


This dimension of the K-shell also represents 
the radius of the ion as a whole, since the dimen- 
sions of the nucleus are negligibly small. With the 
exception of Be’*, equation (1) is confirmed satis- 
factorily (see Table 1 and Fig. 1). 

For all the electron shells occurring after the 
K-level, the dependence on the magnitude of the 
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TABLE 1, Verification of Multiple Proportions 
in Ionic Dimensions [Eq. (4)]. Radii of Rare- 
Gas Type Ions in A 


Element 


calculated 
acc. to (4) 


acc, to 
Goldschmidt 


where the K-level is specially selected, or taking 
(L) into account 


r = 0.2(5 —%) + 0.2n, (3) 


where n is the number of quantum levels in the ion, 
following the K-shell closest to the nucleus. 


Lithium 0.8 0.8 

Saaanen 1.0 1.0 Substituting in this principal quantum number 
Potassium 12 1S m=ntl, we have 

Rubidium 1.4 1.5 

Cesium 1.6 1.6 r == 0.2(5 — y) + 0.2 (m — 1) (4) 
Francium 1.8 ig 

Beryllium 0.6 0.3 and 

see i i on ka Ga) 
Strontium 1.2 1.3 Table 1 gives the results calculated according 
Barium 1.4 1.4 to (4) compared with Goldschmidt's data. The radii 
Pere , i of the ions not given in Goldschmidt's book have 

Ti seats OiG 0.6 been taken from Yatsimirskii [7] for francium, and 
Seundinie 0.8 0.8 from Ormont [8] for radium and thorium. The radi- 
Yteium 1.0 11 us of B’” has been interpolated from B’* and C’’”’ 
Lanthanum 1.2 ee The data obtained by Zackhariasen [9] for actinium 
Actinium 1.4 1.4 have been converted for the Goldschmidt system. 
Carbon 0.2 0.2 The only element to show any pronounced dis- 
aan bie re crepancy is beryllium. All the others, within the 
Titanium 0.6 0.6 eee 

Sale eas ae AS limits of the established accuracy, obey the rule 
Ges 1.0 1.0 which has been found, and this rule is obeyed much 
Thertuin 1.2 1.1 better than, for instance, the rules of Dulong and 


charge is negligibly small, since these shells are 
comparatively remote from the nucleus. They ex- 
hibit a surprising regularity: The "thickness" of 
each of the shells is the same, namely 0.23 AS to 
an accuracy of 0.03 A. Roughly speaking (see dia- 
gram, Fig. 2), an atom or ion "puts on" more elec- 
tron shells of the same thickness as its structure 
becomes more complex. The radius of an ion may 
therefore be expressed by the simple formula 


Petit, Trouton, and Hume-Rothery. 

Similar calculations have been performed, al- 
though with lesser accuracy, for the radii of atoms, 
since the latter, in the case of rare-gas type struc- 
tures, differ from ionic radii approximately by the 
same constant quantity. Here, again, there is the 
repetition of the same quantity 0.20, evidently rep- 
resenting some "universal thickness of the electron 
layer"; as was to be expected, however, the num- 
ber of discrepancies is here greater. Out of 26 
atoms of different elements, there are six discrep- 
ancies exceeding 0.1 A, instead of the single simi- 


Ke NERO =P (2) lar discrepancy in the case of the cations. 
+ ? ? ? ? ? y 
K-layer 
a 29 
Charge 7 4 | 2 | } | “ | 5 ey as Cation 
| t cans \ 
Radius equal | Less than 0 i eee ? é charge 
to 0.2 (5—n) | 0.8 0.6 0.4 OFZ Ox & 7 / 


Fig. 1. Diagram showing the occurrence of constancy and multiple proportions in the 
K-shell closest to the nucleus, The size of the K-shell varies as a function of the ca- 
tion charge (dotted lines), All the dimensions are expressed in an integral repetition 

of 0.2 A. 
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Fig. 2. Diagram of the occurrence 
of constancy and multiple proportions 
in the dimensions of ions. The thick- 
ness of each “electron layer" (quant- 
um level) is 0,2 A. The integral re- 
petition of this level gives the radius 
of any cation of rare-gas type. The 
K-layer closest to the nucleus is spe- 
cially marked. The principal P, O, 
N, M, L layers are marked along a 
vertical. 


Even for rare-gas type cations, however, the 
structure of which is ideally simple and symmetri- 
cal, and for which the ionic radii are most reliable, 
it would be illusory to expect high accuracy, meas- 
ured, for example, in hundredths of an angstrom. 
First, in hundredths, there is nonadditivity, i.e., 
inconstancy of the radii themselves and, second, 
it would be quite unlikely for the volume of the atom 
or ion to increase by an exactly observed increase 
in the number of electron shells of equal thickness. 

The rule which has been found points toa steady 
decrease in electron density with increase in dis- 
tance of the levels from the nucleus, in qualitative 
agreement with wave-mechanics calculations. This, 
however, is unessential for us in the plan in which 
the consideration of the question is developed in 
this communication. It is important to note that for 
all the ideally spherical rare-gas type of ions (with 
one exception), there occurs a quite definite factor 
of multiple dimensions, both for the K-shell closest 
to the nucleus, constructed by repetition of the 
standard unit of 0.20 - and for all the other shells 
(L,M,N,O, P,Q, etc.), built up by an (m — 1)-fold 
repetition of the same unit 0.20 a 


The establishment of this simple rule permits 
us to return once more to the general statement of 
the problem, with which the article began. Summing 
up the foregoing, we may draw the following conclu- 
sions. 

1. Constancy and multiple proportions in crys- 
tallography and chemistry represent a natural con- 
sequence of the corpuscular (atomic) structure of 
matter. 

2. Constancy and multiple proportions also ap- 
pear in atomic (ionic) dimensions, here again,evi- 
dently precisely as a consequence of the corpuscu- 
lar (electronic) structure of the particles (atoms 
and ions) forming the crystal lattices. With the ex- 
ception of the K-layer closest to the nucleus, all the 
other electron levels of cations of rare-gas type 
have a "thickness" of the order of 0.20 i so that 
the dimensions of an ion are determined by the 
(m— 1)-fold repetition of this standard multiple plus 
the radius of the K-layer. 
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The article shows that the thermal expansion of crystals cannot be reduced 

merely to a simple elongation of the crystal in all directions passing through 
any point in its interior. A correct pattern of the thermal expansion of crys- 
tals is given by the transformation of the unit radii of a sphere to the radius- 


vectors of an ellipsoid OB by means of translation vectors AB (Fig. 2). The 
coefficients of thermal expansion of crystals are assumed to be the inter- 
cepts AC and not the intercepts AD. The dependence of the coefficients of 
expansion of a crystal on direction is given by formula (8), by means of 
which it is possible to construct all imaginary surfaces of coefficients of ex- 


pansion of crystals (Figs. 9-15). 


A more complete picture of the thermal 


deformation of crystals calls for an investigation of the surfaces of the coef- 
ficients of thermal displacement BC (Fig. 2) of crystals. 


It is well known that crystals undergo homo - 
geneous deformation on thermal expansion. 
One of its peculiarities is that a sphere, imagined 
to be separated somewhere inside a crystal, is 
transformed by the deformation into an ellipsoid, 
and acube into aparallelepiped. In the homogene- 
ous deformation of plane figures, a circle is trans- 
formed into an ellipse and a square into a parallelo- 
gram. 

Such a definition of the thermal deformation of 
crystals allows the thermal deformation to be under- 
stood in different ways. We shall consider two 
treatments of this phenomenon: one definitely in- 
correct and the other correct; or, better, more 
correct, since it is amenable to further improve- 
ment. We shall commence with the first. 


Incorrect Concept of the Thermal 
Deformation of Crystals 


We assume that the thermal expansion of crys- 
tals occurs in such a manner that every radius of a 
sphere, imagined to be separated inside the crystal, 
is transformed into an identically directed radius- 
vector of an ellipsoid by simple elongation or com- 
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pression. This means that in a plane diagram the 
radius of a circle OA is transformed into a radius- 
vector of the ellipse OB, OC into OD, OE into OF, 
ete. (Fig. 1). 

We assume as principal coefficients of 


expansion, the magnitudes 
ae Mey Chie - DC | 
Oke a) ar 


where a, is the length of the principal horizontal 
semiaxis of the ellipse; @, is the length of the prin- 
cipal vertical semiaxis of the ellipse; r is the length 
of the radius of the circle; and t is the increase in 
temperature. 

If the thermal deformation of crystals actually 
occurred by simple elongation of all the radii of the 
circle, we should have to assume as coefficient of 
expansion of the crystal along any oblique direction 
the magnitude 


PE) eg (ea 


Ct 
OL -t rie & 


where R denotes the length of the variable radius- 
vector of the ellipse. It is not difficult to see, how- 
ever, that the described mechanism of thermal de- 
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Fig. 1. If the deformation of a circle 
to form an ellipse occurs in such a 
manner that any point E on the cir- 
cumference passes to the point F, the 
point G on the side of the square ought 
to pass to the point H. Accordingly, 
the square ought to be transformed in- 
to the curvilinear shown in the diagram, 


formation of crystals, with all its simplicity and 
apparent clarity, is greatly in error, if not from 
the practical standpoint, at least from the theoreti- 
cal standpoint. This becomes perfectly obvious if 
we apply this mechanism not to the deformation of 
a circle, but to that of a square KLMN. In this 
case, it appears that the square is not transformed 
into a parallelogram, but into some curvilinear 
figure, shown in the diagram. To construct any 
point H of this figure, the length of the intercept 
GH must be determined. If we assume that the co- 
efficient ~@' has the value given above, the length of 
the intercept GH may easily be found from the re- 
lationship 


GH  R=s 


On a 


since the length of the radius-vector OG of the 
square is known for each of its points G. If the 
crystal were deformed by thermal expansion in the 
manner indicated, its lattice would be distorted. 
This means that it could have a strictly lattice 
structure at only one definite temperature with 
which, of course, it is impossible to agree. 

We have examined this incorrect treatment of 
homogeneous deformation in such detail because, 
following the footsteps of earlier authors [1], cer- 
tain contemporary crystal physicists [2], including 
in the past the author of the present article [3], have 
unfortunately resorted to it. 

By way of example, we shall cite one of the 
usual discussions based on the incorrect mechanism 


of homogeneous deformation which we have con- 
sidered. We shall consider essentially the deriva- 
tion of one of the most important formulas of crys- 
tal physics, although specifically the same thermal 
expansion of crystals will be borne in mind, 

If we assume that the radius of a deformed 
circle is equal to unity, the equation of the ellipse 
obtained from it may be written in Cartesian coor- 
dinates, referred to the principal axes of the el- 
lipse, as follows: 


2 
ul 


Chal eae = 


oo z. 


we 


Making use of the substitutions 


2, = Reos(ARX;) = + «)c,, 
tom Acostitr,) == (1-4 ec; 


and ignoring squares and derivatives of 
the small quantities @',@,,@,, we re- 
duce the equation of the ellipse to the form 


, 2 2 
Xe) == Cy  MQCp. 


We have obtained this formula bearing in mind 
the deformation of plane figures. If we repeated our 
argument as applied to three-dimensional figures, 
we Should arrive at the formula 


2 2 2 
a! == O4Cy + &yCg + H%5C3, (1) 


enabling us to calculate the coefficient of expansion 
a'from the three known principle coefficients @,, 
@,, @3 for any direction in the crystal defined by the 
direction cosines C,,C9,C3. 


More Correct Concept of the Thermal 
Expansion of Crystals 


This concept is based on the assumption that 
the thermal deformation of a crystal may be de- 
scribed by the system of equations 


, 
Ly = A,X, 


Ly = AyL ) » (2) 
i= ae. 


where xX;,X9,X3 are the Cartesian coordinates of any 
point of the crystal before deformation, xj, xj, x4 
are the coordinates of the same point after deforma- 
tion, and a@,, a), @3 are some constants (positive 

or negative real numbers). It is assumed that the 
coordinates are referred to the system of the 
principal axes of the crystal which theoreti- 
cally may always be found (for monoclinic and tri- 
clinic crystals, the position of the principal axes is 
fixed only for a given temperature). If it is a ques- 
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tion of the deformation of a sphere which we ima- 
gine to be separated inside the crystal, or actually 
ground out of it, we take the center of the sphere as 
the origin of the coordinates. 

-We assume that all the points of interest to us 
are arranged after deformation on the surface of 
this sphere and that the radius of the latter is equal 
to unity. This means that the coordinates x,, X_,X3 
are connected together by the equation 


ai + a +23 = 1. (3) 


Having determined from equation (2) the coordinates 
X4, X9, X3 and substituting the values found for them 
in the equation (8) for the sphere, we obtain the ex- 
pression 


"9 


1 : zs 
geet So, 
ay ay ae (4) 


which represents the equation of an ellipsoid. We 
have shown that the sphere is transformed by equa- 
tions (2) into an ellipsoid. Similarly, it may be 
shown that a cube is converted by them into a paral- 
lelepiped. For this purpose, it is necessary to take 
as starting equation not that of a sphere, but the 
equation for two parallel sides of a cube, and to 
show that by means of equations (2) itis transformed 
into the equations of two planes equidistant from the 
origin of the coordinates and parallel to each other. 

As is known, the equation of a plane may be 
represented in the form 


Ly Lo X32 weed 
Ay, us Ree ua 1, 


where A;,A,,A3 are the lengths of the intercepts of 
the plane with the coordinate axes. If this equation 
is taken as the equation for one face of a cube, the 
cube face parallel to it will have the equation 


zy Ly Ly y 


Ais te AeA 


Substituting in both these equations the values of x,, 
X», X3 found from equations (2), we get 


a 


/ 
%: = | 


eh} 


gies a3A3 


4 

gamer eeat 
Mg Ay an 
' , , 
x Zz, v» 
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These are the equations for two parallel planes equi- 
distant from the origin of the coordinates, which 
was to be demonstrated. 

We shall examine the process of thermal de- 
formation, described by equations (2) in greater 


detail with reference to the plane scheme (Fig. 2). 
We shall take a circle having a radius equal to 
unity. We divide it into 4n equal parts (24 parts in 
the diagram). We connect the points of division to 
the center of the circle; we have a system of radii 
of the circle before its transformation into an el- 
lipse. We now apply to each point of division of the 
circle, and correspondingly to each radius of the 
latter, the transformation (2), assuming x3 = 0.We 
take first the point K with the coordinates x, = 1, 

X, = 0. Multiplying them. respectively by a, and ay, 
we get x{ = a4, x} = 0. This result shows that the 
point K, after deformation while remaining on the 
X, axis, ought to be displaced to the right if a, > 1, 
or to the left, if a, < 1. In constructing the diagram 
we assumed a, = 1.25, a) = 2.00. In accordance 
with this, our point was displaced to the position E 
on the right of K. It follows from the foregoing that 
a,= OE. Similarly, it may be shown that a, = OF. 
As will be seen, the constants a@;, @,, from the 
geometrical point of view, are the semiaxes of the 
ellipse into which the circle is transformed as the 
result of deformation. If we apply our transforma- 
tion to the point A, we find that, after deformation, 
it will have moved to the point B along the dis - 
placement vector AB. In accordance with 


Fig. 2, Transformation of a circle in- 
to an ellipse as described by equations 
(2) and corresponding exactly to the 
observed thermal deformation of crys- 
tals. In this case, the vector OA is 
transformed into the vector OB by 
means of the displacement vector AB, 
which transfers the point A to B. In 
this case, a magnitude proportional to 
the length of the intercept AC is taken 
as coefficient of expansion.in the OA 
direction, 
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this, the radius of the circle OA is transformed in- 
to the radius vector OB of the ellipse. If we com- 
plete this construction for all the points of division 
of the circle, without representing for the sake of 
simplicity the radius vectors of the ellipse, we ob- 
tain the complete deformation pattern shown in 

Fig. 2. It will be seen that this pattern differs sub- 
stantially from that previously considered (Fig.1). 

By using formulas (2), it is easy to calculate 
and construct the diagram for the transformation of 
a square into a parallelogram (Fig. 3). 

Returning to Fig. 2, there is to be seen in it 
the point D, coinciding with the point of intersection 
of the ellipse and the prolonged straight line OA, 
and the point C, which is located at the intersection 
of the same straight line with the normal BC 
dropped from the point B onto this straight line. 
We transfer these and other points of the diagram 
considered to a special diagram (Fig. 4). We de- 
note the angles formed by the straight line OC with 
the axes of the coordinates by ¢;, @, and their 
cosines, which are the direction cosines of this 
straight line, by c,,c,. It is evident that the coor- 
dinates of the point A for OA = 1 are c,,c, and the 
coordinates of point B, in agreement with formulas 
(2) are equal to acy; ayCy. We furthermore denote 
the length of the intercept of the straight line OC by 
GM. 

The diagram shows that the normal equation of 
the straight line GH for this notation may be written 
in the form 


Ca ie ee ee 


Fig. 3, Diagram to show that a square is 
transformed into an oblique parallelo- 
gram by equations (2), 


Since this straight line passes through the point B, 
the equation cited will not be violated if, instead of 
X4,X_ we substitute the coordinates of these points 
a1C1; Cy. After this substitution, the equation 
cited becomes the equation 


a’ = a,Ci + apc}. (5) 


We must here digress slightly from the main 
theme of our article and point out that if, in the 
equation (5) we have obtained,the quantities a,,a, 
are regarded as constants and the quantities a',c,, 
c, as variables, taking into account, of course, the 
fact that the sum of the squares of the direction 
cosines is always equal to unity, 


1{=a+6, (6) 


this equation (5) will define in Cartesian coordinates 
some curve which, after due examination, is found 
to be an oval-like curve, with four points in com- 
mon with the ellipse under consideration (Fig. 2), 
and situated outside the latter with all its other 
points. 

We return to our subject. We subtract equation 
(6) from (5) 


1) ct + (a2 — 1) ¢. 


av—tl=(q 


As we already know, the quantities a, —1 and a,-1 
represent the principal coefficients of expansion Q,; 
and a». If, in agreement with this, we denote the 
quantity a'—1 by a', the last-mentioned equation 
becomes 


OO" me als Lats « (7) 


Fig. 4. Diagram for deduction of Eq. (8), 
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Formula (7) describes the thermal expansion 
of plane figures. If we repeated our considerations, 
having in mind three-dimensional figures, we 
should obtain the equation 


a = %1C + ayc3 + 303" 


(8) 


We now direct our attention to the fact that 
formulas (1) and (8) do not differ from each other 
in outline, despite the fact that the first of them is 
approximate (since, in deriving it we ignored 
squares and derivatives of the quantities @', a,,Q5, 
Qs) and the second is perfectly accurate. This is 
due to the fact that in formula (1), by the quantity 
@' we understood the intercept AD (see Fig. 4), 
which is actually defined approximately by this 
formula, whereas, in formula (8) by the quantity a' 
we understood the intercept AC, which is defined 
perfectly accurately by this formula. 

Taking this into account, and bearing in mind 
the fact that formula (8) has been derived on the 
basis of more exact concepts of thermal deforma- 
tion, we shall henceforth understand by the coeffi- 
cient of expansion @' a quantity proportional to the 
intercept AC and not the intercept AD. 

At first sight, such a definition of the coeffi- 
cient of expansion — in fact not novel — may appear 
to be too artificial, since the point C is situated 
outside the ellipse (i.e., outside the crystal, if be- 
fore deformation it is given the shape of a sphere). 
It is necessary, however, to take into consideration 
the fact that the real transformation of the vector 
OA (Figs. 2 and 4) into the vector OB by means of 
the displacement vector ABoccurs not only 
by extension of the vector OA to the magnitude 
of the intercept AC, but also on account of the 
shift CB, as the result of which the point C coin- 
cides with the point B, where it ought to be. 

Itis very important to note that the formula (8) we 
have derived, being an equation of second degree 
with regard to the variables c,,c,,C3, is an equa- 
tion of fourth degree with regard to the variables 
X1,X_,X3- This may easily be verified by replacing 
the variables c;,C,,¢3,@' by their equivalent quanti- 
ties 


Be a 83 V 2 2 2 
Sar a Sep x x X3. 
a’? a’? ao’? tte 2+ 3 


It follows from this that equation (8), despite 
often-encountered statements, cannot express in 
Cartesian coordinates any surface of a se- 
cond order, in particular the surface of anellipsoid. 

The fact that in individual cases equation (8) 
cannot express an ellipsoid even approximate- 
ly is clear from the following. We know that the 


coefficient of calcite along the principal axis has a 
positive sign, and along directions perpendicular 
thereto a negative sign. This means that in certain 
oblique directions the coefficient of expansion ought 
to be exactly equal to zero, i.e., certain radius 
vectors of the ellipsoid ought also to be equal to 
zero. Such ellipsoids do not exist. 

We should here like to point out that, despite 
the previously often-repeated view [4,5], in the 
thermal expansion of calcite there is not only no 
contradiction of the ellipsoid law, but, on the con- 
trary, a most elegant confirmation of this law, 
since, in this case also, the sphere is transformed 
by equations (2) into an ellipsoid and a circle into 
an ellipse (Fig. 5). 

Form, 


Symmetry, and Antisymmetry 


of Surfaces of the Coefficients of 


Expansion of Crystals 


By means of equation (8), it is possible to find 
from the three known principal coefficients of ex- 
pansion Q;,Q@»,Q3, the coefficient of expansion a! 
for any direction in the crystal, determined by the 
cosines C;,Cy,C3, which affords the possibility, by 
plotting from one point in different directions radi- 
us vectors proportional to @', to construct theo- 
retically all conceivable forms of surfaces of 
coefficients of expansion of crystals for 
1, Q», O3 differing from each other in magnitude 
and sign. 

For crystals of the cubic system, all three 
principal coefficients of expansion are equal to each 
other and all three are positive (Q, = Q, = Q3 > 0). 
In this case, equation (8) assumes the form 


(9) 


a” == 0, (ci + e+ ch) re. 


Fig. 5. Diagram showing that the ellipsoid law for thermal 
deformation remains valid even if one of the principal co- 
efficients has a negative sign. 
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The corresponding surface evidently is a sphere 
with positive radii, a white sphere having the sym- 
metry ~/o +m. 

For Q3 * Qs) = Q4; Gy, G3 > 0, the surface of 
the coefficients of expansion, expressed by the 
equation 


a! 2 oy (CL + C3) + %3C3, (10) 
becomes an oval-like, positive (white) surface of 
revolution, prolate (for a3 > Q,), or oblate (for 3 
< a,) along the X3 axis (Figs. 6 and 7), and, inboth 
cases having the symmetry m+ ~:m. Such sur- 
faces describe the simplest and most frequently en- 
countered case of the thermal expansion of optical- 
ly uniaxial crystals. 

It.was pointed out in the foregoing that, in cal- 
cite, a3 has a positive sign and a, a negative sign. 
In accordance with this, equation (8) for this and 
similar uniaxial crystals will have the form 


a vs — yeh $8) $ ag. 


(11) 


Figure 8 shows the corresponding surface of revo- 
lution having the symmetry m+ ~:m. As will be 
seen, it consists of two ovoid positive (white) re- 
gions and one toroidal negative (black) region. 


Fig. 6. Surface of coefficients of 
expansion of crystals satisfying the 


equation a' = a4(c? + c3) +a4c% 
for a3 > Oy. 


Fig. 7. Surface satisfying the 
equation a" = o4(c? + c3) + a3ck 
for 3 < 4. 
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If all the principal coefficients of expansion are 
positive and all differ from each other (@; ~ Q: * G3; 
Q1,M%,@3 > 0), equation (8) assumes the form 


(12) 


a! = C1 + Oc + M55- 

The corresponding surface having the m- 2:m 
symmetry is shown in Fig. 9. It describes the 
most frequently encountered cases of the thermal 
expansion of the majority of optically biaxial crys- 
tals. 

If one of the principal coefficients, say @3, has 
a negative sign, equation (8) becomes 


(13) 


The corresponding surface (Fig.10) consists of two 
negative (black) ovoid regions, and one elongated 


2 2 2 
A! == H4Ci + Co — HgC3. 


Ef 
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Fig. 8. Surface satisfying the 
equation a" = —a4(c? + cB)ach. 


Fig. 9. Surface satisfying the 
equation a" = a,c? + agc3 + agck. 


Fig. 10. Surface satisfying the 
equation a" = oycP + a,c = O4C%. 
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positive (white) region, having the form of an 
elongated annulus (bread roll). Theoretically, it 
may be predicted that such surfaces, having the 
symmetry m- 2:m, may describe the thermal ex- 
pansion of some biaxial crystals. 

The surfaces described change their form sub- 
stantially when one or two of the principal coeffi- 
cients of expansion become equal to zero. 

For @; = 0, equation (10) becomes 


a” = o303. (14) 
The corresponding surface of revolution, shown in 
Fig. 11, has the symmetry m+ ~:m and repre- 
sents a pair of ovoid (white) regions, osculating at 
one point. 

The same equation (10) for a3 = 0 assumes the 
form 

a! = (C3 + a), (15) 

The corresponding surface of revolution, shown in 
Fig. 12, has the form of a round bread roll (white) 
with two funnel-shaped depressions osculating at 
their apices. The surface has the symmetry 
m:+ ©:m. 

If one of the principal coefficients in equation 
(12) is equated to zero, then (for 3 = 0),.itassumes ; 
the form 


(16) 


2 2 
a! = HC) + C5. 


Fig. 11. Surface satisfying the 
equation a' = ac, 


BS 
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Fig. 12. Surface satisfying the 
equation a" = a4(c? + cB), 


) 
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The corresponding surface, shown in Fig. 13, has 
the form of an elongated bread roll with two funnel- 
shaped depressions, osculating at their apices. The 
surface has the symmetry m- 2:m. 

It remains to consider surfaces in which, for 
equality to zero of one coefficient, the other two 
have different signs. If, in equation (13), we put a 
= 0 and change the sign of c,, it becomes 


a” =— aC] + 903 « (L7) 
The corresponding surface, shown in Fig. 14, rep- 
resents a pair of identical negative (black), and a 
pair of identical, but of another size, positive 
(white) ovoid regions, connected by their apices at 
one point. The surface has the symmetry m: 2:m. 
If, in the last-mentioned equation, the prin- 
cipal coefficients of expansion are antiequal to each 
other, it assumes the form 
a! = 4 (C3 — ci). (18) 
The corresponding surface, shown in Fig. 15, a 
pair of identical white ovoid regions and a pair of 
black ovoid regions, antiequal to the former, con- 
nected together by their apices at one point. The 
surface has the symmetry m- 2:m, and it also 
has the property of antisymmetry (m- 4:m). 


Fig. 13. Surface satisfying the 
equation a" = ayc? + ach. 


ae 


Fig. 14, Surface satisfying the 
equation a = —a,c? + a,c. 


OS 


Fig. 15. Surface satisfying the 


equation a" = a4(c}— cf), 
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We consider the last five surfaces (Figs.11-15) 
constructed theoretically on the assumption thatone 
or two principal coefficients of expansion are equal 
to zero, to be unreal for crystals, but assume that 
it will be possible artificially to produce uniform 
anisotropic materials, the thermal deformation of 
which will correspond to these surfaces. 

In conclusion, we would point out that the pat- 
tern of thermal deformation of crystals we have 
outlined is not final, since we have restricted our- 
selves merely to the investigation of the surfaces 
of the coefficients of thermal expansion and have 
left entirely out of consideration the surfaces of the 
coefficients of thermal shear of crystals. We pro- 


pose to deal with this second problem in a future 
article. 


a 
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A geometrical interpretation is given for the dielectric-constant tensor, which 
gives the Fresnel ellipsoid, and also a geometric interpretation of the recipro- 
cal tensor, which gives the refractive-index ellipsoid. An anisotropic dielec- 
tric may be characterized in terms of these two ellipsoids, in which case a 
simple graphical method gives any two of the vectors D, E, and P from the 


other one. 


The analogy with optics is employed with the refractive indices 


Ng along the ray and ny along the normal to give the ray dielectric constant eg 
(in the direction of the given vector |E) and normal one €y (in the direction of 


the given vector D). 


The following simple relations apply to a homo- 
geneous isotropic dielectric: 


(1) 
(2) 


P= ak; 
cE = E+ 4rP, 


D = (1 + 4nxa) 
where D is electrical induction, P is electrical po- 
larization, E is electric field strength, & is dielec- 
tric constant, and @ is polarizability; € and @ are 
scalars here, so D, P, and E coincide in direc- 
tion. 

Equations (1) and (2) allow of a very simple 
geometrical interpretation. If, for example, we 
mark out vectors D in all directions (keeping the 
absolute magnitude constant) from some point, the 
ends of these vectors will lie on a sphere. The ends 
of the corresponding E vectors will also lie on a 
sphere, so (1) and (2) represent the geometrical 
operation of transformation of a sphere into a 
sphere. In the present case, the radius of sphere 
D will be © times the radius of sphere E. 

An anisotropic dielectric (crystal or homo- 
geneous anisotropic medium) has a tensor relation 
between these quantities: 


P=J7,E, 
D = TE, 


ao 


Double dielectric-constant surfaces are constructed. 


where II, and Il, are symmetric tensors of the 
second rank. The expressions can be put as 


Py =4,E, + OF’, + a 3h 


Cuneo esl Seri ve a eee Rie” cee eees, 


see ena ie) lala ele ae tall ef Walint 9) 


eh 0 16) ey TGF te Beret te de Mian do latte 


in which P,... and D,... aré the components of 
vectors P and D along the coordinate axes. 

fl, and Il, are symmetric (© jk = €4j, ik 
= Q,4); the coefficients of the form €jj; and jj are 
usually termed longitudinal, while €;;, and Qj, (i 
~ k) are called transverse. 

An orthogonal coordinate system may also be 
used to describe phenomena in crystals; the axes 
are chosen to coincide with the symmetry axes of 
the crystal as far as possible. If the crystal has 
less than three mutually perpendicular axes, one 
of the axes is taken as coinciding with one of the 
symmetry axes, the other two passing through 
edges or perpendicular to faces to give an ortho- 
gonal system. Two axes lie in the symmetry plane 
for monoclinic crystals of class m, the third being 
perpendicular to this plane. The coordinate system 
is chosen arbitrarily for a triclinic crystal, in 
general. 
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These tensor relations also allow of a geo- 
metrical interpretation. 


Dielectric Ellipsoids 


Tensor calculus [1] shows that a symmetric 
second-rank tensor may be put in diagonal form by 
transformation of the coordinate system; it also 
shows that such a tensor has a fairly simple geo- 
metrical interpretation. Here we consider only 
the geometrical interpretation of Il, already in di- 
agonal form; there is no difference of principle in 
the consideration of II, and II,. 

Consider the equation of the surface 


Ei-(/7,-) = 1, (3) 


formed by the ends of the variable vectorE in va- 
rious directions. In coordinate form this is 


uly eg ot-lexkss =2r1, (4) 


and is the equation of an ellipsoid (€; > 0); equation 
(4) in canonical form is 


(5) 


The scalar product Il, - E = D of (3) has a simple 
geometrical meaning; if vector E ends at point M 
on the surface of (3), then vector Il; - E = D has 
the direction of the normal to the plane tangential 
to the surface at M (Fig. 1). The absolute value of 
this scalar product is given by 


z | 
| 7.-E| =|D| =: >> 


where ON is the distance along the normal n from 
the center of the ellipsoid to the plane tangential to 
the surface at M. We leave aside the proofs that 
the surface of (8) has these properties, because 
these are given in general form in textbooks on ten- 
sor calculus. 

Consider also the equation of the surface 


D. (Tz". D)=1, (6) 


formed by the ends of vectors D in various direc- 
tions. 

Here, II af is a tensor reciprocal to Il,; if Il, 
takes the form 
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Fig. 1. Determination of the direction 
and magnitude of vector D from the 
Fresnel ellipsoid (section of an ellipsoid 
of rotation by a plane passing through 
the axis of rotation). 


then Iz! takes the form 
eet pat 
ey 
i Page [: 0 - 0 
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Surface (6) has properties analogous to those 
of (3); the scalar product Il;! - D = E has the di- 
rection of the normal to the plane tangential to the 
surface of (6) at M', while the absolute value of the 
product (of the vector E) is given by 


\o"-D\=|El =a, 

where ON is the distance along the normal n' from 
the center of the ellipsoid to the plane tangential to 
the ellipsoid at M' (Fig.2). 

Turning now to (5) and (7), we see that here 
the variables are the components of E and D, the 
constants being €,, &, and &€3. The semiaxes of 
the ellipsoids are 1/Ve&, I/V €,, 1/V€3, and Vey, 
V&, and Vés, respectively. Comparing (5) with 
the equation of the Fresnel ellipsoid 


z 


ae) 
8 


-- 
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where the v are the principal velocities of light in 
the crystal, and comparing (7) with that for the op- 
tical indicatrix 
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Fig. 2. Derivation of the direction and 
absolute magnitude of vector Efor the 
ellipsoid of refractive indices (section 
of an ellipsoid of rotation by a plane 
passing through the axis of rotation), 


== 4. 
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where the nj are the principal refractive indices of 
the crystal, we have, because optics gives 


Cc c 
n=Vé, orien ire 


(where c is the speed of light in vacuo), that (5) de- 
fines the Fresnel ellipsoid and (7) the optical indi- 
catrix. The ellipsoids defined by (5) and (7), or by 
(3) and (6), are called reciprocal, because they are 
defined by tensors II, and nett 

A very important property of these ellipsoids 
is that the ends of vectors D describe the ellipsoid 
of (6) if those of E describe that of (3), the direc- 
tion of D being defined by that of the normal to a 
plane tangential to the ellipsoid of (3) at the end of 
vector E (point M), and conversely: the ends of E 
describe the ellipsoid of (3) if those of D describe 
that of (6), the direction of E being defined by that 
of the normal to a plane tangential to the ellipsoid 
of (6) at the end of vector D (point M'). Figure 3 
shows this construction for a section of an ellipsoid 
of rotation by a plane passing through the axis of 
rotation. 

Let 


Be 3 Si De (8) 
where the * denotes that E * corresponds to D* in 
accordance with ellipsoid (3). Now let a D, equal to 
D** in (6),be equal to D*; if we can show that the 


E** derived from (6),is equal to E*, we have the 
proof of the above assertion. Infact, from 


E* _ Tee Dt — MN ipee . D*, 
and from (8) we have 


E* = I7,° A Whe Z E* = E’, 


Fig. 3. Vectors D, E, and 4nP in an aniso- 
tropic dielectric for the particular case of 
ellipsoids of rotation in section by a plane 
passing through the axis of rotation: 1) 
Fresnel ellipsoid; 2) optical indicatrix. 


which was to be proved. We can also use (6) and 
specify D* to show that this vector and the D** de- 
rived from (3) will be equal, provided that the cor- 
responding vectors E * and E** are also equal. This 
demonstrates the one-to-one relation between the 
vectors D and E defined by (8) and (6). 

These features make it very simple to find one 
vector given the other. 


In addition, 
D=E-}+ 4cP, 


which gives a simple means of finding the vector 
47P (Fig.3). Any homogeneous isotropic dielectric 
may thus be characterized by (5) and (7) as regards 
its electrical parameters. Sections of the ellipsoids 
are used by laying off one of the quantities (D, P, 
or E) in a suitable scale in a given direction and 
then by construction determining the two others as 
to direction and absolute magnitude in the same 
scale. 

The following are some points on the shape of 
the Fresnel and refractive-index ellipsoids. The 
€; are not distinct for all crystal systems, on ac- 
count of the symmetry; the number of distinct &j 
gives us either a general ellipsoid (€; ~ & ~ €3),an 
ellipsoid of rotation (€, = €) ~ €3), or a sphere (€ 
= €) = €3). A prolate Fresnel ellipsoid corresponds 
to an oblate optical indicatrix, etc. 

This is not the only method of considering the 
polarization of anisotropic media or of interpreting 
the dielectric-constant tensor. For instance, 
Shubnikov [2] has treated the tensor relation be- 
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tween two polar vectors in a plane as the trans- 
formation of a circle (whose points are the ends of 
one of the vectors) into an ellipse (whose points are 
the ends of the other vector). This treatment is es- 
pecially valuable in the thermal deformation of a 
crystal, hydrostatic compression, etc. Polariza- 
tion of a dielectric can be considered in this way, 
and the final results must be the same. The two 
methods are thus complementary, and either may 
prove the more convenient for a particular problem. 

The above method of interpretation is suitable 
for a linear dielectric (one with a linear relation 
between field and polarization). A special study is 
needed to determine whether the method can be ap- 
plied to dielectrics of other types. 


Dielectric-Constant Surfaces 


By these surfaces are meant ones formed by 
the ends of the radius vector r from a point O in 
the crystal, the length being equal to the dielectric 
constant. An electrically isotropic body has spheres 
as its surfaces, and these are simple to construct. 

The dielectric constant is a tensor for any 
anisotropic medium; D and E coincide in direction 
only along the principal axes, so the concept of a 
dielectric constant as the coefficient of proportion- 
ality between D and E becomes meaningless. 

We now examine the meaning to be given to the 
dielectric constant of an anisotropic body. 

The refractive index ny is defined in optics as 


where c is the velocity of light in vacuo and vy is 
the velocity along the normal to the wave front. 
This is the refractive index measured in optics. The 
refractive index ng along a ray is defined by 


¢ 6] 
ns = os , 
where vg is the ray speed or the speed of energy trans - 
port. 
It is shown in optics [3] that the velocity along 
the normal is equal to the projection of vg on the 
wave normal: 


vn = Ug-Ccosa, 


where a@ is the angle between the directions of the 
normal and the ray. Then the ray refractive index 
ng is related to ny by 


ng = Ny-Cose. 
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The two indices are clearly the same along the 
principal axes of the Fresnel ellipsoid and optical 
indicatrix: 


ns =—=nn =n. 


Also, the refractive index n of an isotropic medium 
is shown in electromagnetic theory to be related to 
the dielectric constant €: 


n=Vu 


Then we have two € for an anisotropic medium, 
namely €y in the direction of D, and €g in the di- 
rection of E, with 


But 


ng = Ny-COsa, 
so €y Should be related to &g by 


eg = Ey COS? o, (9) 

Condition (9) essentially defines €g and €y. It 
is simple to show that (9) is obeyed if, by &y, we 
understand the quantity that indicates by what fac- 
tor D (as defined by its direction cosines c;,C,,¢3) 
is larger than the projection (in the direction of D) 
of the corresponding vector E ; by Eg we understand 
the quantity that indicates by what factorE is less 
than the projection (on the direction of E) of D 
(Fig. 3). We then set 


See eee 
its ON aw Beco (10) 
pyre D-cos « 
es (11) 
From (10) and (11), 
en Pas: 
Nn = “cost & 


or 


es = En Cos? x, 


which coincides with (9). 
The equations for the surfaces of En and Egare 


as follows. From (10) and (11) we have 
Fale Be D dee ii 
“N =" c0s & E D-E #OD-E’? 
D-B 
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From the construction of the Fresnel ellipsoid and 
optical indicatrix 


Deh 


and, so, 


tn = Dz (12) 
1 
EB? 


(13) 


Expressing in (5) and (7) the components of E 
and D via the direction cosines of these vectors, 
we put these equations in the form 


{ ee e e 
ot (14) 
4 2 2 2 
Be Se Saka Ses: (15) Fig. 4. Double surfaces for eg and éy 
for the particular case of a section of 
Comparison of (12) with (14), and of (13) with (15), surfaces of rotation by a plane through 
gives the axis of rotation: 1) Fresnel ellips- 
oid; 2) optical indicatrix, 
Leech hates sha (16) 
sy ome By em We are indebted to Academician A. V. |Shub- 
d nikov for a valuable discussion of the present 
eg = 8C) + 29C3 + es C9) (17) results. 


where c;,C,, and cs are the direction cosines of D ee 


and E, respectively. Figure 4 shows the method of 1. N.E. Kochin, Vector Calculus and the Prin- 


constructing the surfaces for ey and ég. The di- ciples of Tensor Calculus [in Russian](Moscow, 

electric-constant surface is thus double for an L952): 

anisotropic body; &y and &g are equal only along 2. A.V. Shubnikov, "Symmetry and geometrical 

the principal axes of the Fresnel ellipsoid and op- interpretation of two-dimensional polar 

tical indicatrix. tensors," Izv. AN SSSR, ser. fiz., 13,3,376- 
This method of interpreting the electrical po- 391 (1949). a 

larization can be extended to other phenomena 3. A.V. Shubnikov, Optical Crystallography [in 


described by symmetric second-rank tensors. Russian] (Moscow, 1950). 
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Growth from the vapor state is reported for p-toluidine, with emphasis on the 
relation of tangential growth rate on (001) faces to thickness for elementary 


layers at 1.2% supersaturation. 


fringence, by Senarmont's method, and in reflection. 


The thickness is measured from the bire- 


The results are con- 


sidered to show that there is an adsorbed layer on a growing surface; the 
differences in growth rates are due to surfacé diffusion, which is especially 
important for elementary layers, which grow more rapidly than thicker ones. 
Dislocations also occur at the junctions between parts of a skeletal structure. 


Surface morphology and direct examination 
show that crystals grow via very thin layers. 
Theory shows that more energy is released when 
an ion is attached to an incomplete layer on a NaCl 
crystal than when it starts a new layer, so the 
growth probability for an existing layer is higher 
than the nucleation probability for a new one [1]. 
The current theory of crystal growth from the vapor 
state deals only with the motion of a layer of ele- 
mentary height [2], although observation shows that 
actual layers are tens or hundreds of units high. 

Quantitative empirical relations for layer 
growth are needed as a basis for further advances 
in crystal growth theory. We especially need evi- 
dence on the effects of layer weight on propagation 
rate (tangential growth rate). 

Marcelin (1914) [3] was the first to observe the 
formation of very thin layers on faces for p-toluid- 
ine growing from alcoholic solution. The crystals 
were thin enough to give low-order interference 
colors in reflection, so the thickness of any layer 
was easily measured. The minimum thickness ap- 
peared to be that of two molecules (about 18 A). 

In 1922, Volmer [4] reported layered growth 
on PbI, crystals formed at an interface between 
PbNO; and KI solutions; he found that thick layers 
grow more slowly than thin ones, and he supposed 
that the lower limit of thickness would be of mole- 
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cular order, but he did not measure the thickness. 
Some general observations on growth from solution 
for layers 1600-4000 A thick have been reported by 
Bunn and Emmett [5]. 

We have examined layer growth in p-toluidine 
in the vapor at low supersaturation, which was pro- 
duced in a chamber with a horizontal temperature 
gradient. This consisted of a glass cylinder having 
a heating wire in the base and also a small internal 
heater surrounding a crucible containing the ma- 
terial. This chamber was set up on the stage of a 
microscope, the crystals growing on the flat glass 
top (Fig.1). The bottom heater controlled the gen- 


Fig. 1. Chamber for observing crystals growing 
from the vapor. 1) Microscope; 2) glass cover; 
3) glass base with heater wire; 4) crucible; 5) 
heater; 6) drops of liquid; 7) crystals; 8) thermo- 
couple, 
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eral temperature, while the crucible heater con- 
trolled the supply of vapor. The appropriate hori- 
zontal temperature gradient was thereby set up on 
the cover. The temperature near the surface was 
measured with a copper—constantan couple. Growth 
could be observed around the melting point (44°). 
The layers on (001) faces had rounded outlines at 1- 
5° away from that point (tangential growth rate not 
dependent on the azimuth). Layered growth involv- 
ing liquid droplets was observed at the same tem- 
perature [6], but these could appear on the surface 
only at a certain temperature and vapor pressure. 
A constant vapor pressure was provided by keeping 
the evaporator at a fixed temperature; the tempera- 
ture at the crystal was controlled with the bottom 
heater. The drops and protuberances indicate cer- 
tain definite conditions. The supersaturation was 
calculated from the temperature; in the above case 
it was 1.2%, or about 107. 

Paratoluidine is orthorhombic and has a layered 
structure; the unit cell contains 8 molecules (two 
layers of four) and has the large value of c = 23.5 A, 
The birefringence is very high (0.22). 


The layered growth on (001) faces is very 
prominent. The tangential growth rates were meas- 
ured without reference to the crystallographic di- 
rection, as the layers were rounded. Layers of 
elementary thickness could be observed at x400 on 
account of the large c parameter, as in Marcelin's 
[3] and Kowarski's [7] work. The above conditions 
allowed us to examine layers from elementary 
thickness up to 150c, on account of the high bire- 
fringence. 

The growth rate was measured with a seconds 
clock and an ocular micrometer. 

Heights from 20c to 150c were measured be- 
tween crossed nicols, the thickness being given by 


A = (Ng — Np) 1, 


where A is path difference, (ig. _ Np) is birefring- 
ence, and e,; is thickness. These measurements 
could be made only on stepped crystals 3-10 yu 
thick, giving first- or second-order colors in the 
thin part. An increase of 150c in the second order 
gave a sharp change of color. The thicknesses on 
the two sides of a step gave the height of the step. 
Similarly, the mean thickness of several layers can 
be deduced in the same way, as in Fig.2. 

Smaller thicknesses were measured in reflec- 
tion by Perrin's method [8]: A = 2ne,, where A is 
path difference, n is refractive index for the given 
direction, and e, is thickness. 

Combination of these two formulas gives e;/e, 
= 2n/ (ng — np); the ratio is 15 for p-toluidine, so a 


given path difference is produced by a thickness 15 
times smaller. A crystal only 1 » thick gives no 
color between crossed polarizers, while the height 
of a step that produces a sharp change of color in 
the second order is only 10c. There is an appre- 
ciable color change from the passage of very thin 
identical layers. Measurements were made in this 
way on layers of thickness up to 10 units. 

The results were verified and extended by us- 
ing Senarmont's method with a A/4 plate. (The 
Senarmont compensator is suitable for very small 
path differences.) The path difference A is 


20 


where @ is the angular displacement of the analyzer 
from the crossed position, and A is wavelength. The 
mercury green line (A = 546 my) is equivalent to 

30 A per 1°, which corresponds to 136 A thickness 
(about 6c) for p-toluidine. 

Growth occurs via spiral layers under the con- 
ditions we used [6], but we took measurements only 
at points where the radius of curvature of the growth 
front was not less than 100 yw, in order to eliminate 
any effects from the curvature. The usual formula 
[2] indicates that the tangential rate for a front of 
radius 100 pw differs by only 0.2% from that for a 
straight front if the height is 15c and the distance 
between steps is 10 u. The mean rates were meas- 
ured to 4-5%, so the correction for curvature is 
quite negligible. 

The thickness was reduced to the number of 
elementary layers n = h/c, in which h is height; the 
actual heights were always close to integers on this 
basis, and the rate always showed a distribution re- 
lated to the height, with few-intermediate values. 

There was good agreement between growth 
rates measured in different ways for layers of the 
same thickness (Table 1). The third line is the 
volume deposited on unit length of a step in 1 sec. 

The amount is least for the layer of unit height; 
it is rather larger for 2c but then remains roughly 
constant up to 25c, after which it rises again. This 
difference arises from variation in the spacing be- 
tween successive steps in spiral-layer growth; the 
steps on thin layers are close together, whereas 
those on thick ones are far apart. The influx of ma- 
terial per unit area of face is presumably fixed 
under specific conditions, so the distance between 
steps should be determined by the thickness and 
tangential growth rate. This feature requires more 
detailed examination. 

It is clear that the crystal grows only via prop+ 
agation of layers under these conditions. The 
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TABLE 1 


Height, in units of c 
Growth rate, p /sec 


Volume deposited 
per unit length, 
u2+ 104 


Note, The asterisks denote Kowarski's values for crystals in vapor under conditions similar to ours. The thick- 


nesses have been converted to c units, 


Fig, 3. Propagation of successive layers 
from growth centers. 


sources of the layers are indicated by the current 
theory as lying at the points where dislocations 
emerge, these giving rise to steps that are the 
growth centers for the layers [2]. 

These dislocations (Fig. 3) do act as layer 
courses; we have previously considered the origin 
of the dislocations [6], which arise where parts of 
a skeletal crystal meet. Paratoluidine usually 
starts to grow from the corners, and the points 
where layers meet may be ones where the two lat- 
tices do not fit exactly, so a dislocation arises. 
This provides a highly active center for further 
growth and a source of spirals, as in Fig. 4. Turns 
far from the center are of large radius and produce 
nearly straight growth fronts (Figs. 2 and 38). 

Figure 5 shows the growth rate as a function of 
thickness for straight fronts. 

It is clear that the rate for an elementary layer 
is very different from that for a thick one at the 
same temperature and supersaturation. The differ- 
ence may be due to growth via a mobile adsorbed 


Fig. 4, Dislocations arising at the junction 
of parts of a skeletal structure. 


layer instead of direct from the vapor [9]. Frenkel 
[10] discussed surface fluctuations in relation to the 
observed roughness of crystal faces; he showed that 
there should be such an adsorbed layer, in which 
the atoms or molecules can migrate freely over the 
surface. Barton et al. [2] also assigned an import- 
ant part to surface diffusion in the growth of ele- 
mentary steps. 

The growth rates for straight fronts show that 
surface diffusion does play a very large part in the 
growth of the elementary layers. The diffusion flux 
reaches a limit for a height of 3c, and the curve of 
Fig. 5 is a hyperbola, as Table 1 indicates. The 
proportion entering by surface diffusion becomes 
relatively small for large thicknesses, for then 
much material is received directly from the vapor. 

These processes cause the growth rate to be- 
come independent of the thickness; the curve of Fig. 
5 becomes nearly horizontal above 25c, where it 
tends asymptotically to a constant value. The growth 
rate for an infinitely thick layer (by extrapolation) is 
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Fig. 5. Tangential growth rate as a function of thickness: filled circles, 
Senarmont compensator; open circles, Perrin's method; triangles, meas- 
urements between crossed polarizers. 
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A study is made of the kinetics of the tangential growth of elementary layers of 


various thicknesses on the surface of a growing crystal. 


On the basis of the 


concept of "creeping" molecules adsorbed by the surface, and on the assumption 
of slight supersaturation, a formula is derived connecting the rate of movement 


of the front of the elementary layer with the thickness of the latter. 


son is made with experimental data. 


According to modern concepts, first developed 
by Frenkel’ [1], the crystal face of any substance, 
in equilibrium with another phase of the same sub- 
stance, comprises a multitude of microscopic steps 
of different heights. The "staircase" formed by 
these steps has rises and falls, so that the macro- 
scopic face is plane, and is characterized by simple 
indices. The front of each step is also not astraight 
line, but has "Frenkel' discontinuities" of molecu- 
lar order of magnitude. The molecules adsorbed on 
the face of the crystal form an adsorbed layer in 
equilibrium with the "Frenkel' discontinuities" and 
the medium surrounding the crystal. If the equlib- 
brium is disturbed, a transfer of molecules to the 
solid phase occurs at the positions of these discon- 
tinuities, leading to the formation of plane (super- 
ficial) diffusion streams. The consumption of sub- 
stance in the adsorbed layer is made good by the 
deposition of fresh molecules from the medium sur- 
rounding the crystal. It is thus significant that di- 
rect condensation of the substance to the solid 
phase occurs only on the lateral surfaces of the 
steps, which we shall call the absorbent or active 
surfaces. 

Applying these concepts to steps of the height 
of one molecule, Burton, Cabrera, and Frank [2] 
developed a theory of the growth of crystals from 
the gas phase. Direct microscopic observation of 
the face of a growing crystal, however, shows that 
steps of a thickness of up to 1000 molecules or 
more play a very important part in the growth pro- 
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A compari- 


cess. Furthermore, experiments by Lemmlein 
and Dukova [3] have revealed that the thinnest of 
the propagating layers have a thickness of one 
unit cell of the structure, and this, in the case, 
for example, of the organic substance paratoluid- 
ine, forms four molecules. From this follows the 
necessity to study the kinetics of the movement of 
the fronts of layers of polymolecular thickness. The 
present work examines the possibility of sucha 
movement. 

We shall examine the mechanism macroscopi- 
cally and characterize the absorbent surface of a 
step, the mean profile ab of which is shown in the 
figure, by the coefficient of absorption 8B, the coef- 
ficient of diffusion D, = a2ve"U/KT Where U is the 
activation energy of plane movement, a~ 5 A, the 
size of a molecule, v ~ 1012 sec}, the frequency 
of vibration of a molecule as a whole), and some 
effective length / = 9h. We shall characterize 
the points of the profile by the distance s to the 
center of the profile O. We denote the mean path 
of the adsorbed molecules on this surface by §;, 
and their chemical potential by 4,. The same 
quantities, but referring to the nonabsorbent sur- 
face, will be provided with the subscript 2. The 
chemical potentials of the gas and solid will be wg 
and Ug 

As will be seen from Fig. 1, we are dealing 
with diffusion along the bent surface of an aniso- 
tropic body. Therefore, the equilibrium concentra- 
tions of the adsorbed molecules depend, generally 
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speaking, already on the coordinate of s. Due to 
this, it is more convenient to operate with the 
chemical potential of these molecules, which in 
equilibrium is contant along the entire surface. 

The experiments of Lemmlein and Dukova [3] 
indicate that the front of a step moves at a velocity 
of not more than about 10 u/ sec, which is much 
less than the mean velocity of the molecules diffus- 
ing over the surface (of the order of hundreds of 
meters per second). When considering the diffusion 
of such molecules to the front of the step, their 
velocity may therefore be ignored. 

With regard to the dependence of the velocity 
of movement of the front of the layer on the thick- 
ness of the latter, which is of interest to us, itmay 
be noted that asymptotically for large values of h 
we should have 


v(h) = Vo + a3} /h, (L) 


where Va is determined by the quantity of sub- 
stance arriving on the absorbent surface directly 
from the gas phase, while j is the density of flow of 
plane diffusion. Subsequently, however, we shall 
be mainly concerned with the range of small values 
of h. 

We consider the case where D,; = %, i.e., py 
= const = U4. We shall consider the diffusion pro- 
cess to be stationary. The density of diffusion flow 
of adsorbed molecules to the nonabsorbent surface 
is V9 (dpp/ ds). The number of molecules deposited 
from the vapor per unit of this surface of the crys- 
tal, for small deviations from thermodynamic 
equilibrium, is evidently a! (Ug — Up); here, y and a 
are coefficients of proportionality, connected, in the 
case of low densities n of molecules adsorbed on 
the surface, with the usual coefficient of plane dif- 
fusion D and the mean lifetime of a molecule in the 
adsorbed state t by the relationships 
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Generally speaking, they depend on »,, but for small 
values of 4, this dependence may be ignored. From 
the foregoing, it is clear that the law of conserva- 
tion of matter leads to the following equation for the 
chemical potential on the nonabsorbent surface (s 

> 1/2): 


dis 
Y2 ds? 


— A(t, — pg) = 0. (2) 


We now formulate the boundary condition for s =//2. 
In unit time on the part aO, B1/2(u4,, — us) mole- 
cules are absorbed, (1/2) g — [4.) molecules are 
deposited from the gas, and ay,(du./ds) molecules 
diffuse through the point a. Therefore, 


diy l l 
Ye od oF + (Hg = 100) = Gs me {'g)- (3) 


By virtue of the continuity of the chemical potential, 


Lint alec LenCer 
du. v l 
eget +S (1g — pe) ae (paps): (4) 


By solving (2) together with (4), it is easy to see 
that 


v = Bla (4105 — ps)/h 
= Vo + a86(tg— ps/1 + BU/2V toe. 


We shall now take into account the finitenature 
of the quantity D,. For this purpose, j44.. in (4) 
must be replaced by some mean value 14; the differ- 
ence between these values y44,, — Wy may be repre- 
sented by a series of powers of the ratio of the 
length of the step to the value of the mean path of a 
molecule in the adsorbed state 1/8, commencing 
with a term of first order (forD,;=~, Ss, = ~ and 
[L400 — 44 Ought to disappear). The finite nature of 
D, may therefore be allowed for by adding to the 
right-hand part of (4) terms with successive powers 
of 1/S;. It is clear, however, that in the first 
place, for calculating the coefficients of such anex- 
pansion, some additional considerations of a physi- 
cal nature are required and, in the second place, it 
will not be true for large values of h. It is there- 
fore more expedient to start from a natural model 
of the uniformly absorbing surface of a step, and 
solve the following boundary problem exactly: 


(5) 


d2u, 
1 + 4 (4 g— th) + B (Hs 7 Pi) = 9, Cee 
a 

so +o (Hg— te) =O, s>H 


d l 

foro —10) — =0, for eer y eat Oy recat or rate (6) 
Here, the term 6 (ug — 1) takes into account the 
transfer of substance to the solid phase. 


From this, considering that B= a, we have 


Ye 


l 
v(t) = $1 B(m —n,) ds 


0 
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where op = MoUs: 
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TABLE 1 
h/4a | { | 2 3 4 | 5 | 6 7 | 10 | 45 | 50 | 100 
Yexp, | 8 ee 5 3.8 | a | 25 2.4} 4.5 1.0] 04] 0.3 
u/ sec | a 
Veale? GM Meret a Su ear5 2.9 | 201 2,21) 1.6) fA Oe Ows 
u/sec | 


For Di >”, Y1 ~®, and (7) becomes (5). Here, 
the quantity 6 = 1/ h, determining the mean devia- 
tion of the absorbent surface from the vertical ob- 
viously varies from values of the order of several 
units for small values of h (in the comparison with 
experiment, we considered @ = 4), and approaches 
unity as h increases. At the same time, for low 
values of h, the first term of the sum in (7) plays 
altogether a relatively minor part, and in it we may 
therefore put 6 =1. In the second term, 9 forms 
part of an expression which rapidly tends to a con- 
stant limit as h increases, so that the variation of 
h may here be ignored and this quantity may be re- 
garded as constant. The problem of the strict de- 
termination of the dependence of @ on h is one of 
surface structure and currently has not been 
solved. 

Proceeding to a brief discussion of the expres- 
sion obtained, we note first of all that it is propor- 
tional to 6u. Just as in expression (1), the first 
term represents the contribution of the flow of sub- 
stance from the gas phase directly onto the absorb- 
ing surface, while the second term represents the 
plane diffusion flow. The coefficients in formula (7) 
vary according to the properties of the substance. 
With increase in the energy of sublimation and the 
bond energy per molecule proportional to it, the 
mean length of the plane path of the adsorbed mole- 
cule increases and the coefficient a~ ve~W/2kT 
falls (W is the bond energy per molecule). In view 
of the fact that the first term of (7) is proportional 
to a, while the second is proportional to Va, plane 
diffusion flow plays a relatively larger part in the 
growth of difficultly volatile substances than in the 
growth of volatile substances. A decrease in the 
frequency of vibration v of the adsorbed molecules 
as a whole acts in the same way as an increase in 
the energy of sublimation. 


The second term of (7) shows that those values 
of h, commencing with which the dependence of in- 
terest to us assumes the character (1), are deter- 
mined by the quantity 8 or, more exactly, by the 
ratio between the mean paths s, and Ss, of the mole- 
cules in the adsorbed state on the absorbent and 
nonabsorbent surfaces, respectively. 

The presence of an impurity, i.e., the pres- 
ence of foreign molecules in the adsorbed layer, 


must alter the value of the diffusion coefficients D, 
and D,, and hence the rate of growth itself. 

Since the same kinetic mechanism evidently 
operates in the growth of crystals from solution, 
the results obtained are also applicable to that case. 

A number of circumstances make comparison 
with experiment difficult: first, the value of the 
constants @ and y; ~ vy, are expressed merely in 
order of magnitude through sublimation energy and 
frequency V; second, the only measurements of this 
kind (of G.G. Lemmlein and E.D. Dukova) have 
been made on one substance only,i.e. , paratoluid- 
ine. We considered 96 = 4, W/kT ~ 20, correspond- 
ing to a sublimation energy of ~15 kcal/ mole, and 
S, = 200 a, and hence found @ and yz. It was then 
possible to select in (7) a value of B>a (8B ~ 40 @), 
such that there was sufficiently close approxima- 
tion of the experimental data in the limits of the ap- 
proximate examination. The selectedvalue! of 6u 
determines the additive component in (7), while B 
determines the form of the curve. Table 1 shows 
the calculated and experimental values of the velo- 
cities of movement of the fronts of layers of differ- 
ent thicknesses as calculated and found experimentally. 

It is not difficult to see from Table 1 that there 
is an appreciable departure from law (1) at h <~12 
a, i.e., three heights of the unit cell of the struc- 
ture. Diffusion flow plays an important part up to 
h ~ 300 - 4a. 

We are happy to take this opportunity of thank- 
ing I.M. Lifshits and G.G. Lemmlein for their in- 
terest in the work, and valuable discussions. 
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In investigation of molecular statistical phenomena in solid bodies close to 
phase transition points, measurements were made of the intensity of the 
molecular scattering of light in quartz, over the temperature range 15-600 
°C. The phenomenon of opalescence was detected in the crystal at 573°C, 
indicating intense and reversible formation of small optical nonuniformities 
in the quartz in the vicinity of the phase transition temperature. 


Introduction and Experimental Method 


In the detailed investigation of the molecular 
mechanism of phase transformations in solids, the 
topics to be studied may conveniently include phe- 
nomena bearing on the statistical side of molecular 
processes. The phenomena in this category include 
the molecular scattering of light, the existence of 
which, in crystals, was first established by Lands- 
berg [1]. 

Molecular scattering of light occurs at small 
nonuniformities in the optically pure medium, these 
being continuously created and reabsorbed as a re- 
sult of thermal movements of the molecules. 

The intensity of the light scattered in the me- 
dium is I ~ I,Ae?/d*4, where I, is the intensity of 
the primary beam and A is the wavelength of the 
light [2]. The quantity Ae? = (€ — €,)?, which sup- 
plies a quantitative expression for the factor of 
most importance in considering the phenomenon, is 
the mean-square fluctuation of the optical dielectric 
permeability from its average value €). If we use 
the variables p and T, in the case of an isotropic 
liquid the quantity Ae? can be written in the form 


ae de\2 > OSEAN are 
ss (tar + tar 


lo (1) 
where p and T are the density and temperature of 
the medium, and Ap2 and AT are the correspond- 
ing mean-square fluctuations of these quantities in 


the microvolume under consideration. Experience 
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has shown that in all known cases of scattering of 
light, (&/ 8p)%pAp? >>(8€/AT)AT*. Then, Ac? 
= p*(de/dp )ipkTB | sinceA p2 = p2kTB/v, where k 
is Boltzmann's constant and 8 is the compressibil- 
ity of the medium [2]. The value of Ae? in an aniso- 
tropic medium will not be considered for the mo- 
ment. It is important to note that the quantity Ae? 

is proportional to the compressibility of the medi- 
um, which shows that a relationship exists between 
the intensity of the light scattered and the elastic 
properties of the medium. The quantity AT? 

= kT’/Cy, where C, is the specific heat of the 

bulk material within the fluctuation range. It willbe 
shown below that under solid-body phase-transition 
conditions, the dependence of the scattered-light in- 
tensity on the temperature term cannot be neglected, 
which lends a special character to the present light- 
scattering study in comparison with earlier investi- 
gations. The characteristic relationship between I 
and 4, holding only in the case of optical nonuni- 
formities of size 1 < A, provides an important ex- 
perimental criterion in considering the average size 
of the light-scattering nonuniformities in the medi- 
um. 

From the foregoing it will be apparent that the 
intensity of the scattered light will be intricately 
bound up with statistical phenomena occurring in 
the medium, and can be used in investigating these 
phenomena. It is only necessary to select effective 
conditions in applying the scattering of light to the 
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study of the fluctuation processes taking place dur- 
ing solid-body phase transitions. 

As the first case for investigation, we selected 
the a = £ transformation in quartz (T = 573°).Our 
choice was determined by the high quality of quartz 
single crystals, which have been thoroughly studied 
under the most diverse conditions [38,4]. The 
quartz phase transformation must be considered, 
in the thermodynamic sense, as a transition of the 
second kind, since it is virtually unaccompanied by 
stepwise changes in specific volume and entropy 
[5,6,7]. Supporting the view that the quartz phase 
transition is a A-point, is the increase in sym- 
metry of the crystal on attaining the high-tempera- 
ture form: A33A, >A,6A, [8]. According to re- 
sults given by Perrier and Mandrot [9], the elastic 
constants of quartz pass through a minimum close 
to the transition point, and the refractive index of 
quartz alters extremely rapidly under the same 
conditions [10]. From these facts we may expect, 
according to formula (1), that the intensity of the 
scattered light will increase markedly during the 
phase transformation, and may serve to indicate 


statistical phenomena taking place within the quartz. 


An important consideration in favor of the choice of 
quartz as the experimental material was the availa- 
bility of data on scattering of light in it. This in- 
formation appears in the series of articles [1], 
while an exhaustive experimental and theoretical 
study of the molecular scattering of light in quartz 
was carried out by Motulevich [11], who investi- 
gated the phenomenon up to a temperature of 250°C. 
Thus, the normal course of a projected experiment 
could be reliably checked by comparing it with the 
results of work carried out previously. 


Description of Experiments 


To pursue the objectives stated above, the ap- 
paratus shown diagrammatically in Figs. 1 and 2 
was assembled. A light beam from a mercury lamp 
passes through a quartz crystal held in a furnace. 
The light scattered by the quartz at right angles to 
the primary light beam is registered by a photo- 
electronic multiplier. 

Oriented 20 x 20 x 40-mm parallelepipeds of 
polished quartz were prepared from best-quality 
specimens of optical quartz, not containing twins. ! 
However, the ultimate criterion for the suitability 
of a specimen for further work was the temperature 
—intensity behavior [1] shown by scattered light up 
to 250°, and according to experimental and theore- 
tical data the relationship should be linear; during 
our experiments we were able to compare the be- 


havior with the painstaking measurements made by 
Motulevich [11]. 

In order that the results of the projected ex- 
periments would be of interest with regard to the 
practical study of phase transitions in crystals, the 
experimental temperature conditions needed par- 
ticular attention. It was necessary to keep to amini- 
mum the temperature gradients in the crystals, be- 
cause of the extreme variations with temperature of 
the properties of the crystal close to the transition 
point. Analysis of the appropriate data showed that 
it would be inadvisable to work with a furnace tem- 
perature gradient over 1 deg/cm. This was a very 
harsh requirement for a furnace working at a tem- 
perature of about 600°C, and having three windows 
on the level of the working zone to boot.2 On the 
other hand, it was necessary to bring about a tem- 
perature distribution in the crystal which would sa- 
tisfy the optical demands of the experimental lay- 
out (see below). The furnace had a square 120 mm 
horizontal cross section, was 250 mm high, and 
was welded from 8-mm thick iron (See Fig. 2). Win- 
dows of diameter about 20 mim were arranged at a 
height of 160 mm above the base of the furnace. To 
reduce heat losses through radiation, at 600°, the 
windows were occupied by heatproof glass light 
filters. Inside the furnace, the windows were pro- 
vided with thin lightproof metal sleeves. The outer 
surface, under the Nichrome heating coil and above 
this, was smeared with a mixture of alumina and 
fireclay. 

The removable lid to the furnace was thermal- 
ly insulated by a ceramic fitting containing a heat- 
ing coil. The quartz crystal Q was placed on a 
metal platform P attached to a thin porcelain tube 
PT. The furnace stood on a massive ceramic slab 
CS. A ceramic fitting containing a heating coilwas 
placed inside the furnace on this ceramic slab. 
After the crystal had been placed in position, the 
removable furnace lid and the lower perimeter of 
the furnace were hermetically sealed with fireclay; 
the furnace thus contained a closed thermal zone. 
A metal jacket MJ was placed over the furnace, 
forming an air pocket 50 mm thick between each of 
the outer walls of the furnace and the metal jacket. 
The three heating coils mentioned above were inde- 


2 


Through the kind cooperation of Academician A. V. Shubnikov, 
some of the quartz specimens used were prepared by A. B, Gil'varg 
at the Institute of Crystallography of the Academy of Sciences; to 
these persons the authors are deeply grateful. 


The authors are deeply grateful to V. A. Sokolov for valuable con- 
sultations on the problems involved in achieving the necessary ex- 
perimental temperature conditions, 
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Fig. 1, General layout of the apparatus (horizontal projection). L) Ultra high- 
pressure mercury lamp; Q,, Og, O3) object lenses; F,) heatproof light filter; F,) 
colored light filter; Ny, No) Nicol prisms; PC) photocell; PEM) photoelectronic 
multiplier; M) mirror gate; T) telephoto lens; F) furnace. 


Fig. 2. Vertical section through 
the furnace. MJ) Metal jacket; 
Q) quartz crystal; CS) ceramic 

slab. 


pendently fed from a battery supply. The currents 
through the coils were maintained to an accuracy of 
0.05 A, using an ELA-type ammeter. Tempera- 
tures were measured with a Pt/ PtRh thermocouple 
and a potentiometer; one of the thermocouple junc- 
tions was placed in melting ice. The thermocouple 
and the crystal could be adjusted by microscopic 
amounts in vertical, horizontal, and rotational 
senses, while in the sealed furnace. Temperature 
gradients were measured by altering the thermo- 
couple position under rigidly constant furnace con- 
ditions. Gradients detected in the working zone 
were compensated for by altering the currents in 
the base and lid heating coils. At temperatures near 
573°, it was found possible to keep the vertical tem- 
perature gradient in the crystal zone down to 0.01 
deg/mm. The horizontal gradient in the direction 
of propagation of the light beam was kept down to 


0.03 deg/mm. The highest temperature in the crys- 
tal was that near the entry window, brought about 
by the use of unequal-sized entry and exit windows 
in the furnace. From the above, it will be seen 
that the isothermal layers in the crystal lay prac- 
tically perpendicular to the light beam passing 
through the crystal. This meant that as the furnace 
was heated, the phase transformation region, pass- 
ing along the light beam because of the temperature 
gradient, automatically had to pass by the observa- 
tion window in the furnace and the photoelectronic 
multiplier or telephoto lens and camera placed at 
this window. Finally, some special experiments 
were carried out to determine the temperature — 
time curves for which the thermocouple tempera- 
ture corresponded to the temperature of the crystal 
alongside it. The furnace was heated extremely 
slowly; a complete experimental cycle (15° ~—600° 
15°) occupied more than 72 hours. The rate of 
temperature change near the phase transition point 
was about 0.3 deg/h. 


Experimental Results: The 
Opalescence Phenomenon 


The final measurements on the relationship be- 
tween the intensity of light scattered by quartz and 
the temperature of the quartz, were made on three 
specimens. Two of the quartz specimens came from 
the same deposit, while the third was obtainedfrom 
another source. Several experiments were carried 
out on each specimen. Investigations were made of 
various combinations of directions of primary and 
scattered beams in the quartz, using different po- 
larization conditions; these are indicated in the 
table. 
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TABLE 1 


Direction 
of propagation 


Separate measurements of the individual inten- 
sity components were necessary to determine the 
depolarization of the scattered light, this being an 
essential parameter in considering the anisotropy 
of the scattering medium. The sum of the quantities 
Ix and ly (in the second case) is given because, when 
the scattered light is propagated along the z axis, 
the optical activity of the quartz leads to a differ- 
ence appearing between the values of Ix and ly. 
Measurements were made of each of the values I,, 
ly, l, for the total light of the mercury lamp, for 
the 4360, 4078, and 4047 A group of mercury lines, 
and for the 5460 A line. From the relationship be- 
tween the intensity of the scattered light and its 
wavelength, the dimensions of the optical nonuni- 
formities could be estimated, this being one of the 
salient points of our experiment. It was important 
to use light filters, so that the percentage of para- 
Sitic light in the complete beam of light being re- 
gistered would be lower for the shortwave part of 
the spectrum than for white light. The results ob- 
tained from measurements using a violet light fil- 
ter were therefore the most reliable. These are the 
ones shown in the figures which follow. However, in 
our apparatus, because of the constancy of the weak 
parasitic light, the relationships found between the 
scattered light intensity and the temperature were 
exactly the same for blue, white, or green light. 

In Figs. 3, 4, and 5 are shown the results of 
our measurements, in the order in which they are 
listed in Table 1. 

The x axis of each graph shows the tempera- 
ture, and the y axis the ratio IT/Inq°, i.e., the in- 
tensity of the light scattered by the quartz at tem- 
perature T over the intensity at room temperature. 
A correction for the temperature-independent para- 
sitic light,in our experiments not accounting for more 
than 0.3I,9°, was introduced by the usual method of 
extrapolation of the linear part of the curve to ab- 
solute zero.* In Fig. 5, the maximum of I'p/In9°, 
corresponding to the phase transition point, was in- 
serted separately, since the main experiment had 
been accidentally stopped at the peak of the maxi- 
mum, and we did not feel it was expedient to linkup 
results of subsequent tests on a single point of the 
curves with the results of the other experiments. 
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Fig. 3. Plot of intensity of scattered light against 
temperature (case 1 in Table 1). @ — Heating; O — 
cooling; x —repeat run; dashed line — theoretical 
curve (for the 573-600° range, this lay exactly on 
the experimental curve). 
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Fig. 4. Plot of intensity of scattered light against 
temperature (case 2 in Table 1), @— Heating; O— 
cooling; x — repeat run, 


3No correction was introduced into our results for the Raman effect 


because of its uncertain value under these conditions, The thermal 
radiation of the furnace was, of course, allowed for by making adi- 
rect measurement of the photoelectronic multiplier current when 
the mercury lamp was switched off. Thanks to the use of heatproof 
filters, even at 600° the thermal radiation of the furnace had only 
5% of the intensity of the light scattered by the quartz at this tem- 
perature, 
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Fig. 5. Plot of intensity of scattered light against 
temperature (case 3 in Table 1). x, @ — Heating. 


The slopes of the linear parts of our curves were 
all within 4% of Motulevich's results [11], obtained 
by the photographic photometry method in the tem- 
perature range 20-250°. 

Figure 6 shows typical results for the depolari- 
zation—temperature relationship of the scattered 
light, with the temperature plotted along the x axis 
of the graph and the depolarization A = Ix /I, along 
the y axis. The results refer to case 1 in Table 1, 


but similar results were obtained in the other cases. 


In Figs. 3, 4, and 5, the maxima on the curves 
are marked with an arrow and the figure 1.4 - 104. 
This indicates that over a 0.1-deg temperature 
range, close to 573°, there is a particularly rapid 
increase in the intensity of the scattered light, this 
reaching 10,000 times its room-temperature value 
I,92. This intensity jump is due to a specific phe- 
nomenon, quartz opalescence, occurring near the 
phase transition temperature. Visual and photo- 
graphic observations of the crystal, carried out in 
parallel with the scattered-light intensity measure- 
ments, showed that the "boundary" of the a and £ 
phases is an area of marked optical nonuniformity 
in the crystal. Ata temperature of 573°, it ap- 
peared as if a strongly light-scattering band of mist 
(I ~ 1.4 + 10*Iy90) began to move slowly from the 
hotter end of the quartz‘ to the colder end. On one 
side of the mist band lay a-phase quartz, and on 
the other side lay G-phase. As the temperature of 
the furnace was raised further, the mist band tra- 
veled to the opposite end of the crystal, and the 
whole crystal became transparent and weakly light 
scattering again. 

On cooling the furnace, the mist band again ap- 
peared, at the cooler end of the crystal, and then 
moved away from this point toward the hotter end 
of the crystal. The band extended across the whole 


THE a=fB TRANSFORMATION 95 


A% 
30 


20 


10 


J 


100 200 300 400 500 600 T°C 


Fig. 6. Plot of depolarization of scattered 
light against temperature, 


width of the crystal, and was 0.5 to 3 mm thick. We 
were able to vary its thickness by altering the hori- 
zontal temperature gradient in the quartz. The 
temperature range in the band was ~01°. By sta- 
bilizing the furnace temperature, we were able to 
hold the mist band in the center of the crystal for 
3h. In Figs. 7,8, and 9. are shown a series of 
photographs, showing the primary beam of light in 
the quartz before the appearance of the mist band 
at the phase boundary (exposure 1 h), the intersec- 
tion of the mist band and the primary light beam 
for a horizontal temperature gradient of 0.3 degper 
cm (exposure 1 sec), and another photograph of the 
latter for a gradient of 1 deg/mm. The photo- 
graphically very unfavorable relationship between 
the intensity of scattered light in the mist band and 
that in the transparent parts of the crystal, on both 
sides of the band, meant that it was not possible to 
photograph the primary light beam satisfactorily in 
both the mist band and the adjacent transparent 
crystal phases simultaneously. In the last two 
photographs, therefore, the outlines of the primary 
light beams have been picked out with dashed lines. 
The phenomenon in question occurred in all the 
quartz Specimens investigated, and could be repro- 
duced several times with any specimen. X-ray 
tests on specimens before and after heating did not 
show up any changes in their structures as a result 
of the treatment. There are thus no irreversible 
changes occurring in the crystal which are asso- 
ciated with the phenomenon we observed. The de- 
polarization of the light scattered by the mist band 
was ~6%, while that of light scattered by quartz 
under normal conditions is ~12%, which indicates 


‘Tt will be recalled that the crystal virtually contained only a hori- 
zontal temperature gradient, directed along the primary light beam. 
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Fig. 7, Photograph of the prim- 
ary light beam in quartz (expo- 
sure 1 h), 


Fig. 8. Photograph of the inter- 
section of the primary lightbeam 
and the mist band in quartz (ex- 
posure 1 sec; horizontal tempera- 
ture gradient ~0,03 deg/mm). 


Fig. 9. Photograph of the inter- 

section of the primary lightbeam 
and the mist band in quartz (ex- 

posure 1 sec; horizontal tempera- 
ture gradient ~1 deg/mm), 


that the structure of the mist band is more weakly 
anisotropic than that of a normal quartz crystal 
(see Fig. 6). 

Visual observations of the mist band in white 
light showed that when the band was viewed from 
the beam entry side, scattered blue light predomi- 
nated, but when viewed from the exit side, more 
yellow light was scattered. Such behavior is typical 
for the passage of light through a strongly scatter- 
ing medium. In order to estimate the size of the op- 
tical nonuniformities in the mist band, it was neces- 
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sary to verify that the i~4 law, for the intensity of 
scattering of light of different wavelengths, was 
obeyed. Over the whole experimental temperature 
range, beginning at room temperature, we moni- 
tored the scattering intensity ratio for blue to green 
light. The value of this ratio was ~5. Under our 
conditions, the value of this ratio for the direct 
light of the mercury lamp was 1.6. For the light 
scattered by the mist band, this quantity was in- 
creased to a value of 10. So that the deductions to 
be drawn from this value should be unambiguous, 
we measured in our apparatus the ratio of intensi- 
ties of scattered blue and green light for pure benz- 
ene, which has an undisputed molecular light- 
scattering character. Benzene, which is a liquid 
which scatters light extremely strongly, was 
chosen because its scattering was of roughly equal 
intensity to that observed in the mist band under 
consideration. The value for benzene of the ratio 
we were interested in was 4.2, in our apparatus. 
This provided grounds for the assumption that in 
the mist band we were dealing with optical nonuni- 
formities smaller in size than the wavelength of the 
light. The increase in the ratio of intensities of 
blue to green light scattered by the mist band, above 
the corresponding value for scattering by benzene, 
was due to the fact that the relative proportion of 
parasitic direct light from the lamp was less in the 
experiments with the strongly light-scattering mist 
band than in the experiments with benzene. 


Discussion of Experimental Results 
General Conclusions 


The experimental results described above can 
be most simply discussed by following the natural 
division of the whole temperature range investigated 
into a series of intervals, characteristic of the 
phenomena described. 

Over the temperature interval 20-400°, alinear 
relationship was established between the intensity 
of scattered light I, and the temperature T. In 
Motulevich's paper [11], this linear relationship 
was followed up to temperatures of 250°. In view 
of the fact that theoretical calculations also lead to 
a linear function for the explicit relationship be- 
tween I and T, our results indicate that up to 400° 
the elastic and elasto-optical parameters of quartz 
which, with the temperature, determine the inten- 
sity of scattered light, are weakly temperature-de- 
pendent. 

The temperature intervals 400-573° and 573- 
600° are two regions in which the relationship be- 
tween I and T is not linear. Both the sharp increase 
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and the rapid fall in intensity of scattered light 
within these two intervals are in accord with other 
experimental data, indicating that quartz takes on 
quite different properties in the vicinity of its phase 
transition point. Thus, the sharp drop in the elas- 
tic modulus E of quartz, observed in [9] (to 1.7 of 
the original for E]|z and to 2.5 times the original 
for E,z), points to an increase in the compressi- 
bility of quartz at about 573°, leading to greater 
fluctuations of density in the medium. In Fig. 3, 
the dashed curve shows the theoretical relationship 
between I and T for quartz, allowing only for the 
change in compressibility which is due to the varia- 
tion with temperature of Young's modulus. The path 
of the dashed curve is in good agreement with the 
experimental results. In constructing the theoreti- 
cal curve, it would be necessary to take into ac- 
count not only the dependence of the compressibility 
on temperature, but also the relationship between 
refractive index and temperature, i.e., the tem- 
perature term in formula (1). The measurements 
made by Rinne and Kolb [10] established a relation- 
ship between n and T, but under the natural condi- 
tions of their experiments they measured this value 
at constant pressure, while formula (1) has the ex- 
pression (9€/9 T)p. The transformation from 

(de /8 T)p to (de/ 8 T)p does not present any mathe- 
matical difficulties, since (8&/9T)p») = (d€/98 T)p 

+ (9€/ dp)q(a/ Bry), where o is the coefficient of 
volumetric expansion for the substance.> However, 
the value of (0& /Op) p was not measured for the 
temperature range of interest to us, which forms 
an obstacle to the conversion.® Qualitatively, how- 
ever, the relationship between € and T will undoubt- 
edly also have an effect in the light-scattering be- 
havior under study here, where, because of the 
large values of (d€/9 T)p » fluctuations in tempera- 
ture AT will develop [see (1)]. 

Let us return now to discussion of the opales- 
cence phenomenon arising in quartz in the immedi- 
ate vicinity of the phase transformation point. Our 
experiments were the first to detect the phenome- 
non of reversible opalescence in a solid body. It 
turns out here, when we compare our measure- 
ments with opalescence intensity measurements 
made at the critical points of various materials [12, 
13], that opalescence is stronger in quartz than in 
certain liquids in their critical state. The very fact 
that opalescence appears is an indication of the 
abrupt formation of small optical nonuniformities 
in the quartz. The explanation of the origin of these 
nonuniformities is directly derived from the quite 
unambiguous interpretation given above of the non- 
linear portions of the I—T curves. In a case such 


as this, we must recognize that opalescence occurs 
as a result of a special intensification of fluctuations 
in dielectric permeability at the phase transition 
point. But, since we have no grounds for assuming 
that, in the solid state the compressibility may at- 
tain values similar to those at the critical point of 
the substance, we must look for the cause of forma- 
tion of large optical nonuniformities, above all, in 
large values of (0&/ 8T)p, which provide for the 
possibility of temperature fluctuations appearing 
over small volumes. Calculation shows that the 
value of (AT?) 1/2 is of the same order of magnitude 
as the mist-band temperature range in our experi- 
ments, and further calculations lead to a value for 
the intensity of the scattered light which is accept- 
able in explanation of the effect observed. This 
point of view corresponds, in our opinion, to treat- 
ing the phase transition in quartz as a A point. 

In this phenomenon, it is undoubtedly signifi- 
cant that the w= f transition in quartz is close to 
the critical Curie point on the line of phase transi- 
tions of the second kind [8]. 

We will observe that the behavior shown in our 
experiments, and the visual picture of a single in- 
termediate zone between the two phases, can beas- 
sumed to solidly confirm the existence of a single 
phase transformation in quartz near 573°. Wecon- 
sider it necessary to make this observation, since 


-Steinwehr [14] has asserted that other intermediate 


phases exist between the a and # phases of quartz. 
Steinwehr's statement was based on interferometer 
measurements of the difference in path of two light 
beams, one of which passed through a quartz plate. 
The step-type relationship between the temperature 
and this path difference was.interpreted by the 
author as reflecting a series of phase transforma- 
‘ions. It is perfectly obvious that measurements of 
»ptical path differences in quartz close to the phase 
transformation point, without preliminary detec- 
tion and careful localization of the mist band, will, 
in general, be meaningless and cannot lead to any 
well-founded conclusions. 

In conclusion, we will list the following basic 
results of the work reported. 

1. An experimental method has been devised 
which permits quantitative study of the intensity of 
molecular Rayleigh scattering of light in quartz in 
the 15-600° temperature range, and also allows lo- 


*For greater physical clarity, we will again use the relationships for 
an isotropic medium. 


SThere is a further difficulty in that the variables p and T are not 
statistically independent like p and T. 
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cation and observation of the boundary region be- 
tween the two crystal structures of quartz. 


2. For the first time, a A-type intensity maxi- 


mum has been shown to exist in the scattering of 
light near the phase transition point of a solid body. 
From the existence of this maximum, it was de- 
duced that fluctuations in the dielectric permeabil- 
ity of quartz were intensified as its temperature 
approached the phase transition point. 

3. For the first time in a solid body, the phe- 
nomenon of molecular opalescence was detected 
and quantitatively investigated, in a crystal lying 
near its polymorphic transformation temperature. 
On the basis of the experimental results, an 
examination was made of a possible molecular me- 
chanism for the abrupt formation of optical non- 
uniformities in a crystal at its phase transition 
temperature. ! 
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The effective speed of a computing machine is 
presently only 2-3 times that of the manual method 
in two-dimensional Fourier syntheses, as experi- 
ence here and elsewhere has shown.! One reason 
for this is that at least four workers are needed in 
the machine method, in the later stages at least, 
whereas one person can do all the work in the man- 
ual method. Computing machines are also excep- 
tionally costly, and we have very little access to 
foreign machines, so it would seem that the strips 
method is likely to remain with us as the basic me- 
thod for a long time yet. This makes it extremely 
desirable to rationalize the work; the number of 
operations is not reduced, but the strain on the 
worker is reduced. 

One of the tiring operations in summing the 
numbers in the columns arises from the need to 
add the black figures but subtract the red. The ne- 
gative values on the strips may be replaced by the 
difference from 100, whereupon all the numbers 
are simply added; it merely remains to subtract 
from each column the appropriate number of 100s 
(the number used in the column). 

Red ink (or a bar over the number) is used al- 
so for these differences, but the color is ignored 
in the summation, the red entries merely being 
counted at the end to give the number of 100s to be 
subtracted. 

This treatment is of considerable value if the 
numbers are added by means of a manual or elec- 
trical adding machine such as the SASL- Pz by 
Rheinmetall (and then it is convenient to add up 
three columns at once); but it does not work so well 
with those who prefer oral counting, because the 
sums increase continuously in this method and give 
a large final number, whereas the sequence of plus 
and minus values usually gives a relatively small 


final sum. 


Oe 


The following is a new cosine strip for h= 1 
and amplitude 100 (sixtyfold division of the edge of 
the elementary parallelepiped): 


100 C 1 = 100 99.98 95 91 87 81 14 67.59 50 41 51 21 10 0 


O06 L696 OS HE IWEHME L210 FT 7 oO 
The first 16 numbers (X, X,,X,..., X45) correspond 
to the positive side of the strip, and the last 16 
(X15, X4g,---+»X39) to the negative side; it is assumed 
that the long strip will be folded along the dotted 
line, the second half being glued to the first in that 
position (the numbers in the second set of 16 are 
best written upside down if this is actually done). 


100 99 98 95 91 87 81 74 67 59 50 41 31 21 10 O 
99 90 
98 79 
95 69 
94 59 
87 50 
81 44 
74 33 
67 26 
59 19 
50 13 
44 9 
34 5 
24 B 
10 T 
0 90 79 69 59 50 44 33 26 19 13 9 5 23 4 0 


Scheme No, 1 


1 Machines give a very great saving in time in more routine pro- 
cesses, €.g., refinement of coordinates. 


100 


TABLE 1. Cosines 
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The strip for h = 2 (positive side) will contain 
the zero, second, fourth, etc., entries up to the 
30th on the previous long strip provided that the 
final zero on the first half and the first zero on the 
second half will read identically as number 15, and 
both be black. The strip for h = 3 similarly con- 
sists of the zero, third, sixth, and so on; position 
10 will contain No. 30 of the basic strip. No. 11, 
and so on, may be obtained by adding to this a 
similar strip of 31 entries (Kg), X31,...,Xg9) with a 
mirror-image array, namely red figures first and 
black ones last; here, again, identical red figures 
at position 30 on the first strip and the zero on the 
reversed one will be given the same number, name- 
ly 30. 

Position 11 on the strip for h= 3 has entry 
33; position 15 has entry 45, which finishes h = 3. 

The double strip with 60 entries is the cosine 
part of the fundamental, whose length equals aj, 


the side of the unit cell. The h = 2 and h = 3 strips 
give harmonics corresponding to A= a;/2 and A 

= a,/3. All the strips for the harmonics are de- 
rived from the double strip if this is folded to 
bring entry 60 into coincidence with No. 1 (scheme 
No. 1), the strip for h = 4 being obtained by taking 
the entries zero, fourth, and so on up to the 60th, 
which occupies position 15. No. 64 reproduces (in 
sequence) No. 4, and so on. 

The formula for the strip for a given h may be 
written as follows on the basis of this very ele- 
mentary method. Let xh, xb, a tae ,xh be the num- 
bers appearing at the positions on such a strip;then 


x= Xp XP = ap oe = ye ke a ee ee 
ceeds 60, we take the remainder from division by 


: : OEE 
EBL may be written as > Xn > Xm? 
m = nh (mod 60). 

The cosine is symmetrical about argument T 
[cos (7 +@) = cos (tT — Q)], so we may take as m 


where 
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TABLE 3. Cosines 
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TABLE 3. Sines 
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the least in absolute value of the remainders from 
division by 60, no matter whether this is positive 
or negative, i.e. , we always have one of the first 
30 entries on the strip for h=1. 

Table 1 gives the cosine strips for amplitude 
100 and h = 0,1,2,...30;the values will be altered 
proportionately for other amplitudes. The positive 
sides are shown, the negative sides bearing (in red) 
the differences from 100 of the black numbers on 
the positive side. Experience over more than 15 
years has shown that it is unnecessary to use more 
than two significant figures in the entries. 

The sine strips are derived from the cosine 
ones by 90° shift; Table 2 gives these for h from 1 
to 30. 

The strips are drawn up as follows. We com- 
pile two tables for h = 1 with 100 lines having am- 
plitudes F from 1 to 100. On one side we insert 
poe = F cos mb for the positive amplitudes, and on 
the other (rear side), the same for the negative 
ones. All those on the rear side are negative, so 
the differences are written; but the red 100 is most 
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simply replaced by a black zero. One completed 
sheet is cut up horizontally and gives 100 strips for 
h = 1; the other is cut up into 16 narrow bands along 
vertical lines. The positive sides of the latter are 
numbered 0 to 15; the negative ones, 15 to 30. These 
are laid out in accordance with these numbers and 
xh = ea to form tables for h = 2,3,4, and so on, 
which are then cut up along horizontal lines. 

Further economy is attained by restricting the 
strips to steps of unity for amplitudes 1-10, etc., 
for multiples of ten; any amplitude with an inter- 
mediate value will correspond to two strips. Fifteen 
strips for each h suffice if two significant figures 
are used for differences as well as sums. 

Table 3 gives the complete set of 15 strips 
for h=1 (divisible once horizontally and twice 
vertically) . 
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1. N.V. Belov, Vestn. AN SSSR, 11,75-80 (1954). 
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Convolution is becoming increasingly used in 
structure analysis and scattering theory (see [1] 
for review); this is capable of giving new results, 
and also of giving known results more simply. This 
note deals with the latter aspect. 

The f curves are defined by the Fourier inte- 
gral over the distribution p (r) of the scattering 
matter in the atom at rest: 


f(s) = \o(n) et") doz = F{p(r)} (1) 


(F is the Fourier operator). Consider the distribsu- 
tion p7-(r) in an atom subject to thermal motion 
whose effect on the average distribution of the cen- 
ter of gravity about the equilibrium is a known 
probability function w(r); here we envisage the ran- 
dom part of the thermal motion, which is not de- 
pendent on the bonds between the atoms. We nor- 
malize w(r) via 


| w(x) dv, = f (2) 


An atom displaced from its equilibrium position to 
point r' has a distribution p(r—r'). The probability 
of finding it in that position is w(r'). So py(m is to 
be found by integrating p(r—r')w(r’), over all r': 


or (r) = \e(r—r’) w(t’) doy. (3a) 
But this is the convolution of p(r) and w(r): 
er (r) = pu (r) (3b) 


in complete accordance with the physical meaning 
of the convolution operation. The convolution theo- 
rem states that the Fourier integral of a convolu- 
tion of two functions is the product of the Fourier 
integrals of each of the functions: 


F {ow (r)} = F {o(r)}-F {w(v)}. (4) 


Then we have the result that the atomic tem- 
perature factor (the scattering power of an atom in 
thermal motion) is f(s) - f(s), namely the prod- 
uct of the f(s) of (1) by the temperature factor 
f p(s) as defined by the Fourier integral of w(r): 


fr (8) = F fw (r)} = \ w(x) ec) doy. , (6) 


Here the f (s) of (1) and the f 7.(s) of (5) have been 
put in the most general form, which imposes no re- 
strictions on p (r) and w(r). It is usually assumed 
that both of these functions have spherical sym- 
metry, in which case (1) and (5) are expressed in 
spherical coordinates; an ellipsoidal form has been 
used for w(r) in some structure studies. A spheri- 
cal Gaussian form for w(r), namely w(r) = (272) 

- exp (2? /2u}, gives us, from (5), the Debye tem- 
perature factor f p(s) = eM, w= '4u?s?, where w2 
is the mean-square displacement from the equilib- 
rium position. 
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This communication deals with one of the spe- 
cial questions of the propagation of elastic waves 
in anisotropic media — the case of the propagation 
of ultrasonic waves in the direction of the two-fold 
(Y) axis in a crystal of the tetragonal system. 

It is known from theory [1,2] that three ultra- 
sonic plane waves may be propagated in this direc- 


tion —a longitudinal wave and two transverse waves. 


The velocities of these waves are determined, re- 
spectively, by the elastic constants C41, C44, and Ce¢. 
We shall here confine ourselves only to the trans- 
verse waves. The directions of the displacements 
in transverse waves are perpendicular to eachother 
and parallel to the crystallographic axes Z and X, 
where Z is the four-fold axis and X the two-fold axis. 
In the case of ammonium dihydrogen phosphate, C4, 
= 8.30'- 10!" and eg, = 6.04 * 10'°dyn/cm?,i.e., a 
wave with a direction of displacement along Z is 
propagated at a lower speed than a wave with a di- 
rection of displacement along X. 

By means of a Y-cut quartz plate, a short 
ultrasonic pulse of transverse vibrations of a fre- 
quency of 1.67 Mc is sent into a single crystal spe- 
cimen of ammonium dihydrogen phosphate, 60 mm 
long. The ultrasonic pulses passing through the 
specimen are received by a second Y-cut plate, are 
amplified, and are fed to the plates of a cathode- 
ray tube. The Y-cut quartz plate irradiates and re- 
ceives the transverse waves with a fixed direction 
of displacement, and, in this sense, it is analogous 
to a polarizer in optics. In what follows we shall 
use the terminology adopted in optics, both for 
brevity of the discussion and to emphasize the very 
close analogy between the phenomena studied and 
similar phenomena in optics. 

This analogy has already been noted in the 
literature [3]. 
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Figure 1 shows oscillograms of pulses, which 
have passed through the specimen with change in 
the direction of displacement of the excited oscilla- 
tions relative to the Y and Z axes. It will be seen 
that if the wave excited in the crystal is polarized 
parallel to the X or Z axis (la and 1c), one pulse is 
received, the time of propagation of which depends 
on the length of the specimen and the value of the 
elastic constant. The direction of displacement in 
each of the pulses may be determined by rotation of 
the receiving plate. In the crossed position of "po- 
larizer" and “analyzer,” the pulses disappear. Ifthe 
excited wave is polarized in the direction of the bi- 
sector between the axes, the screen shows two 
pulses, propagating at different speeds (Fig. 1b). 
This is direct confirmation of the existence of two 
"permitted" directions of displacement of trans- 
verse waves fora given direction of propagation. 

As in the case of the propagation of light per- 
pendicularly to the optical axis of a uniaxial crystal, 
two plane transverse waves are propagated in one 
direction at different speeds. 

Figure 2 shows the variation in the oscillo- 
grams received when the receiver (analyzer) is ro- 
tated for constant diagonal position of the radiator 
(polarizer). Figure 2a repeats Fig. 1b, but with 
lower amplification. If the direction of the oscilla- 
tions of the receiver coincide with one of the per- 
mitted directions of displacement, sometimes one, 
sometimes the other wave, is received. It will be 
seen from Figs. 2a and 2c that the patterns are 
practically the same for parallel and crossed 
Nicols. Similar to this, in optics for diagonal po- 
sition of the crystal, the patterns with crossed and 
parallel Nicols in monochromatic light are similar, 
while in white light the color changes to the com- 
plementary color. 
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Fig. 1. Oscillograms of ultrasonic pulses, which have passed through 
a specimen of ammonium dihydrogen phosphate,as a function of the 
direction of displacement in the excited wave relative to the per- 
mitted directions X and Z, "Polarizer" and analyzer" parallel. 
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Fig, 2. Variation in oscillograms on rotation of the receiver n around 
the direction of propagation Y. The radiator u excites an ultrasonic 
wave with a direction of displacement along the bisector of the angle 


between the X and Z axes. 


These phenomena may be observed not only in 
crystals of the tetragonal system, but in any crys- 
tals with directions in which two purely transverse 
waves are propagated at different speeds. This re- 
quirement is met by any of the three orthogonal 
axes in crystals of the rhombic system, any direc- 
tion perpendicular to the six-foldaxisin crystals of 
the hexagonal system, and the [110] and similar di- 
rections in crystals of the cubic system, etc. For 
cubic crystals, in the [110] direction, the differ- 
ence in the speeds of transverse waves increases 
with increase in anisotropy of the crystal (c4, — 

— Cig — 2Cy4). 


We thank 1I.M. Sil'vestrov for discussing the 
questions examined. 
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It is well known that when crystals dissolve, 
multifaceted depressions, called etch figures, are 
formed on their faces. The formation of etch fig- 
ures is due to the different solubility of crystals in 
different directions and is connected with the ini- 
tial stage of solution of the crystals. In manycases, 
a knowledge of the etch figures enables a correct 
conclusion to be made regarding the symmetry ofa 
crystal. 

For the production of etch figures, unsaturated 
mother liquor is used. However, it is not always~ 
possible to obtain etch figures in an unsaturated so- 
lution. It is often necessary to select special sol- 
vents, which alter considerably the form of the etch 
figures. 

In the new method of producing etch figures in 
an ultrasonic field, the crystal is placed in satu- 
rated mother liquor, and is exposed to short-time 
irradiation (of the order of 2-3 min) in a weak 
ultrasonic field. In this way, etch figures have been 
produced on the octahedral face of a crystal of pot- 
ash alum (Fig. 1), on the cleavage plane parallel to 
the pinacoid of a crystal of benzophenone, and on a 
piece of a potassium bichromate crystal of random 
form. 
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For producing the ultrasonic field, a magneto- 
striction generator with a power of 1 kW and a fre- 
quency of 40 ke was used. 

The idea of producing etch figures by means of 
ultrasonic vibrations was suggested by A. V.: Shub- 
nikov. 


Fig. 1. Etch figures on the octahedral face of a potash alum crystal. 
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It is known that the optical spectra of crystals 
depend on the electron polarizability of the ions of 
which they are composed. The polarizabilities 
(refractions) of the ions depend on their radii [1]. 
The melting points of crystals are functions of 
their lattice energies, and are hence also connected 
with the ionic radii [2]. Since the refractive in- 
dices and melting points of ionic crystals dependon 
a single common source, these two characteristics 
ought to be interlinked with each other. The pres- 
ent article is devoted to establishing this link. 

The table gives the refractive indices (upper 
line) and melting points (lower line) of alkali hal- 
ides. The refractive indices in all cases were ex- 
trapolated to » = © by the Wulff method [3] in order 
to avoid the effects of dispersion. 

These data are shown graphically in Fig. 1 
(points enclosed by circles indicating the limits of 
experimental error). 

We see from the graph that the dependence of 
the No of the alkali halides on their melting point is 
linear and inverse. 

This type of relationship is quite understand- 
able. On passing from the fluorides to the iodides, 
and maintaining the type of structure, there is an 
increase in the interatomic distances and conse- 


TABLE 1 
Anions 
_ Cations 
I Cl br I 
Na 1.320 [4] 1.526 [°| 1.613 Fy 4725) (F| 
980 [>| 804 [>] 760 [5] 663 [°] 
K chests 0) 1.474 [7] {036 [7 1.628 [7] 
856 [>] TIA YT 730 [8] 681 [>] 
Rb 1.390 [7| ETD pat 1.529 [7] 1.610 [7] 
760 [?| 726 |} 683 [°] 642 [5| 
Cs 1.468 ["] down fF 1.538 [4] ga i] 
680 [5] 646 |°} 636 [>] 621 [>] 


quent fall in the crystal-lattice energy and melting 
point. In the same direction the polarizability of the 
anion increases: 


Rp sr '298 om? (Rey 2 534 
Roy = 11.32,andRy- = 16.63 [7] 


and hence so does the refractive index. 

The rectilinear nature of the relationship is due 
to the approximately additive character of the op- 
tical and energy properties of ionic crystals. In the 
already-cited paper by G.B. Bokii and S.S. 
Batsanov, it was shown that the refraction (and 
hence the refractive indices) of alkali halides is 
connected with their composition in an additive 
fashion. On the other hand, the crystal-lattice ener- 
gy may also be represented as a sum of terms cor- 
responding to the individual ions [9]. Hence, it is 
clear that the melting point of ionic crystals may 
also be connected with their composition by an ad- 
ditive function. 

Thus, both the refractive indices and melting 
points of ionic crystals are determined by functions 
having an additive character and depending finally 
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on the same thin, the interatomic distances. This 
also produces the simple proportional relation be- 
tween n and the mp of the crystals. 

In analytical form the relationship may be de- 
scribed as follows: the refractive indices of alkali 
salts are 


Ma = 2.420 — 0.001039 - 1° 
nx = 2.500 — 0.004250 - 1°, , 
ry = 2.650 — 0.001553 - 2° 
Nos = 3.017 — 0.002208 - 1... 


It is interesting to note that an equation of 
exactly the same type connects the boiling point of 
organic compounds with their refractive indices 
[10]. Hence, we may conclude that the refractive 
index of chemical substances depends linearly on 
their melting point or boiling point. 
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The generally accepted mechanism for the de- 
velopment of solid solutions of nonmetals (nitrogen, 
oxygen, etc.) in metals, in particular those formed 
by ds elements (titanium, zirconium,etc.), pro- 
vides for the splitting of the nonmetal molecule in- 
to atoms on interaction with the metal. 

At the same time, in regard to the hydrides of 
the same ds elements, since the time of Hagg's 
work [1], opinion has strengthened (see, for 
example, [2]), that unsplit hydrogen molecules 
might penetrate directly into the lattice of a metal, 
with the formation of a structure analogous to that 
of ThC,. We indicated in [3] that this model (pro- 
posed by Hage, for example, for the "dihydride" of 
zirconium ZrH,) was improbable, since the effect 
of the considerable expansion of the metal lattice 
as the hydrogen penetrated could scarcely be com- 
pensated from the energy side by forces arising be- 
tween the H, molecules and the Zr crystal. We in- 
dicated that, from general theoretical considera- 
tions, we should expect the formation of structures 
with the pairs of atoms (H,) splitup, i.e., with 
isolated hydrogen atoms. Otherwise, if Hagg's 
model for the structure of ZrH, were confirmed, 
quantum chemistry would have the difficult theo- 
retical problem of handling the mechanism and na- 
ture of this linkage. This question is of fundament- 
al interest. 

In the light of what has been said, we cannot 
fail to turn attention to the fact that in contemporary 
literature the description of the structure of zirco- 
nium hydride on Hagg's model is still continuing 
(see [4], 1955). 

It therefore seems appropriate, for elucidating 
the true state of affairs, to consider the important 
study of the dihydrides (dideuterides) of thorium 


and zirconium by x-ray, neutron-diffraction, and 
magnetic methods made by Rundle, Shull, and 
Wollan [5]. These authors found a tetragonalized 
cell belonging to the fluorspar type. Unfortunately, 
the authors, after obtaining extremely important 
experimental data regarding the Me—Me and Me—D 
interatomic distances, did not publish the D—D dis- 
tances or consider the question of the relation be- 
tween these data and the Hagg model. The absence 
of a sketch of the unit cell in this paper may be re- 
sponsible for the fact that subsequent discussion of 
this structure is usually based on the "old" unit 
cell which, after [5,6], ought to be considered as 
unfounded. 

We therefore think it appropriate to demon- 
strate the connection between the fec unit cell of 
fluorspar (Fig. 1 shows two such cells) and the be 
tetragonal unit cell of the dideuterides (dihydrides) 
of zirconium and thorium [6], and to derive the cor- 
responding interatomic distances D—D. 

The positions of the metal are 000 and 4,%,%. 

Analysis of intensities led the authors of [5] to 
the conclusion that of the five possible variants of 
hydrogen-atom distribution, the correct one was 
000; 4%, %, %A+"%, “A: 4, 4. This description 
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(unusual in our own literature) is equivalent to the 
following: os whe Wey Othe Joy ee, Loe. /. 

The unit cell of the dihydride is shown in final 
form in Fig. 2. Without deformation, the c/a 
ratio following from the fluorspar structure is V2 
= 1.414. In zirconium dihydride it equals 4.443; 
3.520 = 1.26, in thorium dihydride 5.03:4.1=1.23. 
Thus, some contraction takes place along the c 
axis. As seen from Fig. 2 and from the coordi- 
nates of the D atoms, the D—D distance denoted by 
the broken line equals c/ 2 = 2.224 A for zirconium 
dideuteride and 2.515 A for thorium dideuteride. 
The Me—D distance in the first case equals 2.09 
and in the second 2.41 ne as in [5]. 

Thus, the D—D distance in the dideuterides is 
three times larger than in the D, molecule (where 
it equals 0.74 A), i.e., there are no hydrogen 
molecules or groups of hydrogens in the lattice. 

It should be noted that, as shown in [6], the 
structure of ThC, also differs from that given by 
all authors prior to [6], and, in particular, the 
C-—C distance reaches 1.5 A as against 1.2 A in 
acetylene. 
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